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ABSTRACT 


MacRootLocus,  a  computer-aided  design  program,  was  developed  as  an 
analysis  and  design  tool  for  linear  feedback  control  systems.  A  variety  of  programs 
are  presently  available  for  the  IBM-PC  and  IBM  mainframe.  The  Apple  Macintosh 
offers  just  a  one— parameter  root  locus  method  capability.  MacRootLocus  supports 
both  one—  and  two— parameter  root  locus  methods  on  the  Apple  Macintosh.  It  is 
written  in  the  computer  language  Turbo  Pascal  which  is  the  native  language  of  the 
Apple  Macintosh  and  is  designed  with  the  same  user-friendliness  and  standard 
interface  philosophy  the  Macintosh  was  designed  for. 
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I.  MACROOTLOCUS  SYSTEM 


A.  INTRODUCTION 

MacRootLocus  system  is  a  program  written  for  the  Apple  Macintosh 
computer.  It  was  designed  to  allow  a  user  to  analyze  and  design  linear  feedback 
control  systems.  It  is  written  in  the  computer  language  Turbo  Pascal  which  is  the 
native  language  of  the  Apple  Macintosh  and  designed  with  the  same 
user-friendliness  and  standard  interface  philosophy  the  Macintosh  was  designed  for. 

The  Macintosh  gets  most  oi‘  its  input  from  the  user  through  the  'mouse'  which 
is  a  small  cigarette  pack— sized  control  that  is  moved  across  the  table  much  like  a 
pencil  across  paper  to  move  the  cursor  or  pointer  on  the  computer  screen.  Rather 
than  typing  in  commands  like  the  IBM,  the  Macintosh  lets  you  select  the  function 
you  want  performed  with  the  mouse.  Since  the  commands  are  not  typed  in.  the 
commands  do  not  have  to  be  remembered  as  with  other  computers.  It  is  the 
user-friendliness  that  sets  MacRootLocus  apart  from  other  system  analysis 
programs.  Prior  computer  experience  is  not  required  to  use  MacRootLocus.  A  few 
minutes  to  learn  how  to  use  the  mouse  and  pull  down  menus  is  all  that  is  needed. 
All  dialogs  are  easy  to  use  and  there  are  on-line  Help  menus.  So  the  first-time  user 
can  get  desired  results  without  using  trial  and  error. 

Apple  Macintosh  models  supported  by  MacRootLocus  are  the  Macintosh  SE 
and  the  Macintosh  II  system. 

MacRootLocus  was  tested  with  several  examples  and  is  now  available  to  any 
user  on  the  Naval  Postgraduate  School  Controls  Laboratory  computers  under  the 
icon  of  MacRootLocus  system. 
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B.  MATHEMATICAL  CONCEPT 

The  root  locus  method  is  a  graphical  technique  for  determining  the  roots  of  the 
closed— loop  characteristic  equation  of  a  system  as  a  function  of  the  static  gain. 
Consider  the  general  feedback  control  system,  as  shown  in  Figure  1.1. 


G(s) 

1  ^ 

k 

- 

H(s) 

Figure  1.1  Closed-loop  for  Control  System 

In  order  to  find  the  roots  of  the  characteristic  equation,  it  is  required  that 

1  +  G(s)H(s)  =  0.  (l.I) 

Of  course,  Eq.  (1.1)  may  be  rewritten  as 

G(s)H(s)  =  -1.  (1.2) 

Since  s  is  a  complex  variable,  Eq.  (1.2)  may  be  rewritten  in  parametric  form 

|G(5)H(s)|  =  1  (1.3) 
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L  G(s)H(s)  =  (2k-l)x 


(1.4) 


where  k  is  an  integer. 

For  a  specific  value  of  s  to  be  a  root  of  the  characteristic  equation,  it  must 
satisfy  both  Eqs.  (1.3)  and  (1.4).  Since  the  roots  are  those  values  of  s  that  satisfy 
both  equations,  then  the  root  points  are  points  where  these  curves  intersect. 

The  original  development  of  the  root  locus  method  was  concerned  with  the 
determination  of  the  locus  of  roots  of  the  characteristic  equation  as  the  system  gain, 
K,  is  varied  from  zero  to  infinity.  It  appears  that  the  root  locus  method  is  a  single 
parameter  method;  fortunately  it  can  be  readily  extended  to  the  investigation  of  two 
or  more  parameters. 

The  characteristic  equation  of  a  dynamic  system  may  be  written  as 

a  s"  4-  a  S"'*  +  ....  +  aS  +  a  =0.  (1.5) 

n  n-l  10  ^ 

Clearly,  the  effect  of  the  coefficient  ai  may  be  ascertained  from  the  root  locus 
equation 


a  S 

1+ - * - .  (1.6) 

a  S"  +  a  S"-‘  +  .  .  .  .  +  a  S2  +  a 

n  n-l  2  0 

The  parameter  of  interest,  A  (in  MacRootLous),  can  be  isolated  as 


a  S"  -t-  a  S"-* 

n  n-l 


-h _ +  (a  -  A)S"-'’  +  AS"-'^-> 

n-q 

-l-aS-t"aS  =  0. 

1  0 


+  .  . 


(1.7) 
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For  example,  a  third-order  equation  of  interest  might  be 


4-  (3  +  A)S2  +  3S  +  6  =  0.  (1.8) 

In  order  to  ascertain  the  effect  of  the  parameter  A,  we  isolate  the  parameter  and 
rewrite  the  equation  in  root  locus  form  as  shown  in  the  following  steps; 

S3  +  3S2  +  AS2  +  3S  +  6  =  0  (1.9) 

1+ - — - =  0.  (1.10) 

S3  +  3S2  +  3S  +  6 

Then,  to  determine  the  effect  of  two  parameters,  we  must  repeat  the  root  locus 
approach  twice.  Thus  for  a  characteristic  equation  with  two  variable  parameters,  A 
and  B  (in  Mac  RootLocus), 

a  S"  +  a  S"-'  +  ....  +  (a  -  A)S"-‘»  -f  AS"-^*^  +  .  . 

n  n-1  n-q 

+  (a  -B)S"-' +  BS"-^-^  +  . .  +  aS  +  a  =  0.  (1.11) 

n-r  I  0 

The  two  variable  parameters  have  been  isolated  and  the  effect  of  A  will  be 
determined,  followed  by  the  determination  of  the  effect  of  B.  For  example,  for  a 
certain  third— order  characteristic  equation  with  A  and  B  as  parameters,  we  obtain 

S3  + S2+ BS+ A  =  0.  (1.12) 
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In  this  particular  case,  the  parameters  appear  as  the  coefficients  of  the  characteristic 
equation.  The  effect  of  varying  B  from  zero  to  infinity  is  determined  from  the  root 
locus  equation 


1  + 


BS 

S3  +S2  +  A 


=  0. 


(1.13) 


One  notes  that  the  denominator  of  Eq  (1.13)  is  the  characteristic  equation  of  the 
system  with  B  =  0.  Therefore,  one  first  evaluates  the  effect  of  varying  A  from  zero 
to  infinity  by  utilizing  the  equation 


$3  +  S2  -f  A  =  0 


(1.14) 


rewritten  as 


1+ - - - =0  (1.15) 

S2(S  +  1) 

where  B  has  been  set  equal  to  zero  in  Eq.  (1.12).  Then,  upon  evaluating  the  effect 
of  A,  a  value  of  A  is  selected  and  used  with  Eq.  (1.13)  to  evaluate  the  effect  of  B. 
This  two-step  method  of  evaluating  the  effect  of  A  and  then  B  may  be  carried  out 
as  a  two-parameter  root  locus  procedure.  First,  we  obtain  a  locus  of  roots  as  A 
varies,  and  we  select  a  suitable  value  of  A;  the  results  are  satisfactory  root  locations. 
Then  we  obtain  the  root  locus  for  B  by  noting  that  the  poles  of  Eq.  (1.13)  are  the 
roots  evaluated  by  the  root  locus  of  Eq.  (1.15). 

In  Mac  RootLocus,  this  design  approach  is  used  to  calculate  the  roots  for  the 

plot. 
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II.  OPERATING  THE  MACROOTLOCUS 


A.  BASIC  MACINTOSH  USE 

This  document  must  serve  as  a  user's  manual  as  well  as  a  technical  explanation 
of  the  program  and  its  capabilities.  For  this  reason,  the  following  brief  explanation 
of  Macintosh  use  is  included.  It  is  by  no  means  a  substitute  for  the  Apple 
Macintosh  Users'  Manual  but  it  will  contain  enough  information  for  the  beginner  to 
be  able  to  use  Mac  Root  Locus. 

1.  Basic  Operation. 

The  Macintosh  has  a  finder,  a  special  application  you  use  to  organize  and 
manage  your  document  and  to  start  other  applications.  You  use  the  finder  every 
time  you  start  your  Macintosh,  or  whenever  you  move  from  one  application  to 
another. 

The  Macintosh  screen  looks  like  a  light  gray  desktop,  rather  than  a 
textual  list  of  commands  and  responses.  The  desktop  simulates  the  working 
environment.  It  is  initially  clean,  displaying  small  graphic  images,  called  icons, 
with  short  titles  directly  under  them  for  each  disk  presently  being  used  and  a  trash 
can  in  the  lower  right  corner.  An  icon  is  an  image  representation  of  an  application 
document  or  a  control  to  a  usable  function.  They  offer  quick  recognition  as  to  the 
type  of  item  they  describe  and  are  easier  to  identify  than  lists  or  directories  of  file 
names  with  extensions. 

The  main  interface  between  the  user  and  Macintosh  is  the  mouse.  The 
Macintosh  responds  instantly  to  every  movement  you  make  with  the  mouse.  You 
can  start  applications  and  get  documents,  work  on  them,  and  pul  them  away  again 
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just  by  moving  the  mouse  and  pressing  the  mouse  button. 

There  are  three  techniques  for  the  mouse;  pointing,  clicking  and  dragging. 
Moving  the  mouse  moves  the  cursor  or  pointer  on  the  screen  in  the  same  direction. 
Positioning  the  pointer  on  an  item  is  called  pointing  to  it.  The  mouse  is  used  to 
select  various  items  or  icons  on  the  screen.  You  select  any  item  or  icon  to  let 
Macintosh  know  this  is  what  you  want  to  work  on  next.  You  select  icons  by  using  a 
technique  called  clicking.  As  you  click  the  icon,  it  becomes  highlighted.  This 
highlighting  shows  that  you  select  it.  This  mouse  can  also  be  used  to  drag  across 
something.  This  means  the  button  is  pressed  and  held  while  the  mouse  is  moved. 
This  action  is  called  dragging.  When  dragging  the  mouse  across  the  screen,  a 
rectangle  is  outlined.  When  the  button  is  released,  everything  within  the  rectangle 
is  now  highlighted  and  selected.  Dragging  also  refers  to  moving  items  on  the  screen. 
This  is  done  by  pointing  to  an  item,  pressing  and  holding  the  button  and  moving 
the  mouse.  This  will  also  move  an  outline  of  the  icon  selected  and  when  the  button 
is  released,  the  icon  is  moved  to  the  new  location. 

Whenever  you  work  with  Macintosh,  you  tell  it  two  things:  what  you 
want  to  work  on  and  what  you  want  to  do.  First,  you  tell  the  Macintosh  what  you 
want  to  work  on  by  selecting  it  as  you  have  been  doing  with  icons  on  the  desktop. 
Then  you  tell  the  Macintosh  what  you  want  to  do  with  the  selection.  You  usually 
do  this  by  choosing  a  command  from  a  menu. 

Along  the  top  of  the  screen,  in  the  menu  bar,  are  titles  of  the  menus. 
Pressing  the  mouse  button  while  you  are  pointing  to  a  menu  title  causes  the  title  to 
be  highlighted  and  a  menu  to  appear,  much  like  a  window  blind  being  pulled  down. 
The  menu  contains  commands  you  can  carry  out  on  what  you  have  selected. 
Commands  that  you  cannot  use  right  now  appear  dimmed  in  the  menu.  When  you 
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release  the  mouse  button,  the  menu  disappears. 


To  choose  a  command  from  a  menu,  you  use  the  same  dragging  technique 
you  used  to  move  icons.  As  you  drag  through  a  menu,  each  usable  command  is 
highlighted  in  turn.  If  you  change  your  mind  about  choosing  a  command,  move  the 
pointer  off  the  menu  and  release  the  mouse  button.  Nothing  is  chosen  unless  you 
release  the  mouse  button  while  one  for  the  commands  is  highlighted.  You'll  follow 
this  same  pattern  whenever  you  work  with  the  Macintosh;  'select'  some  information, 
the  'choose'  an  action  for  it.  For  example  in  MacRootLocus,  if  you  double  click  the 
mouse  on  MacRootLocus  icon,  the  icon  will  be  selected  and  MacRootLocus 
application  will  start  running.  You  can  see  the  menu  bar  along  the  top  of  the  screen 
with  the  greeting  message  in  the  center  of  the  screen.  Now,  you  can  select  and 
choose  the  items  to  use  MacRootLocus  for  your  design. 

2.  Manipulating  Window 

The  window  is  the  area  that  displays  information  on  the  desk  top.  You 
view  application  documents  and  folders.  Folders  are  a  way  of  storing  and  organizing 
things,  much  like  in  a  file  drawer.  Double  clicking  on  a  folder  will  open  a  window 
that  displays  its  contents.  Folders  can  be  located  in  disks  or  within  other  folders. 
Usually  8  to  12  windows  are  the  maximum  that  can  be  open  at  any  one  time.  When 
several  windows  are  displayed  at  once,  they  will  usually  overlap.  If  the  window  you 
want  to  view  is  partially  covered  by  another  window,  pointing  anywhere  in  your 
window  and  clicking  will  select  it  and  make  it  the  active  window.  The  active 
window  is  always  in  front  of  all  the  others.  It  is  the  place  you  want  the  next  action 
to  happen,  such  as  move  or  select  icons  or  open  other  folders. 

The  active  window  can  be  identified  by  its  highlighted  title  bar  with 
narrow  horizontal  lines  on  either  side  of  the  title.  The  active  window  also  usually 
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has  a  close  box  (to  close  the  window)  at  the  top— left  corner;  at  the  top^-right  corner 
is  a  zoom  box  that  expands  the  window  until  it  nearly  covers  the  screen.  On  the 
bottom— right  corner  is  the  size  box  you  use  to  change  the  size  of  a  window.  As  you 
click  or  drag  these  boxes,  you  can  get  the  function  that  you  want. 

3.  Handling  Input 

Most  information  in  the  form  of  data  or  text  is  entered  through  the 
keyboard.  The  keyboard  includes  character  keys,  numerical  keys,  direction  keys  and 
other  special  keys.  The  return  key  tells  the  computer  that  the  data  just  typed  in 
should  be  accepted  now.  In  MacRootLocus,  the  apple  key  in  combination  with 
another  key  is  often  a  shortcut  to  choose  a  command  from  a  menu.  The  tap  key  is 
used  to  move  between  data  input  points  in  the  dialog  for  entering  data.  The 
direction  key  is  important  in  MacRootLocus.  It  will  be  used  for  some  word 
processing  applications  in  a  small  box  of  the  plot  window.  This  will  be  further 
explained  in  the  next  section. 

When  the  Macintosh  requires  information  from  you,  it  will  display  a 
dialog  box  like  the  one  shown  in  Fig  2.1.  It  will  tell  you  exactly  what  data  is  needed 
and  shows  you  where  to  enter  it.  Data  insertion  points  are  small  boxes  in  the  dialog 
box  that  let  you  type  in  numbers  or  text.  There  are  usually  several  such  insertion 
points  in  each  dialog  box.  The  tab  key  lets  you  move  from  one  point  to  another  to 
enter  data.  Usually  there  will  already  be  data  in  an  insertion  point  box.  This  is  the 
default  data.  You  can  change  it  if  you  want  but  you  do  not  have  to.  When  you 
enter  new  values  in  an  insertion  point  box,  they  will  become  the  new  default  values. 
Data  insertion  boxes  can  also  be  selected  by  pointing  and  clicking  with  the  mouse. 
Clicking  the  mouse  between  two  characters  in  the  box  will  let  you  insert  characters 
between  them.  Double  clicking  a  box  will  highlight  the  entire  number  or  an  entire 
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Figure  2.1  Sample  Dialog 


word  if  text  is  entered.  You  can  also  select  all  or  portions  of  numbers  or  text  l)y 
dragging  the  mouse  across  the  text  you  want  to  select.  When  all  or  part  of  a  text  is 
selected,  it  will  be  highlighted.  It  can  be  removed  by  using  the  delete  or  backspace 
key,  or  it  can  be  replaced  by  typing  in  whatever  you  want.  After  all  data  in  the 
insertion  box  is  correct,  you  can  enter  the  data  by  hitting  the  return  key  or  by 
clicking  the  OK  button.  If  you  click  on  Cancel,  any  changes  to  the  data  in  the 
insertion  boxes  will  not  be  saved  and  the  operation  which  called  the  dialog  box  will 
be  canceled.  If  you  type  in  data  that  does  not  apply  to  the  insertion  box,  such  as 
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typing  in  letters  in  the  AMin  Gain  insertion  box  of  one— parameter  plot  data  dialog, 
the  Macintosh  beeps.  It  informs  you  that  you  inserted  the  wrong  input  and  you 
should  correct  the  input  that  you  just  typed. 

4.  Printing  Out 

In  order  to  print  out  your  work,  there  are  a  few  ways  available  in 
MacRootLocus.  One  way  is  to  select  the  print  command  in  the  File  menu.  This 
allows  the  user  to  get  a  hard  copy  of  any  plot  displayed  by  MacRootLocus. 

Another  way  is  to  do  a  screen  dump  which  will  print  the  contents  of  the 
active  window  immediately.  This  is  done  by  holding  the  command  key  and  then 
typing  the  number  '4'.  This  is  a  fast  way  to  get  print  out  of  a  plot  in 
MacRootLocus.  You  can  also  create  a  MacPaint  document  by  pressing  the 
command  and  shift  key  and  typing  the  number  '3'.  You  can  take  up  to  10  of  these 
'snapshots'  and  can  then  alter  them  with  MacPaint  for  transferring  to  a  word 
processor  for  lab  writeup. 

B.  MENUS  AND  DIALOGS 

There  are  five  menus  for  MacRootLocus:  Apple,  File,  Edit,  Plot  and  Help. 
Actually,  the  File,  Plot  and  Help  menus  are  used  to  get  the  plot.  These  have  to  be 
used  in  order  which  will  be  explained  later. 

The  Apple  and  Edit  menu  are  not  directly  used  by  MacRootLocus,  but  they  are 
used  in  order  to  follow  the  standard  Macintosh  programming  philosophy  for 
user-friendliness.  This  allows  the  experienced  user  to  easily  adapt  to  a  new 
program  since  many  of  the  operations  are  already  familiar.  Also,  this  is  to  allow  for 
easy  interaction  with  various  other  programs  and  desk  accessories  [Ref.  1]. 

As  mentioned  earlier,  MacRootLocus  presents  commands  in  menus  you  pull 
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down  from  the  menu  bax.  As  soon  as  you  choose  the  command  you  want,  the  dialog 
box  appears  for  each  command.  This  allows  you  to  insert  input  data.  Now  more 
details  for  menus  and  dialogs  in  MacRootLocus  will  be  explained. 

1.  Apple  menu 

The  Apple  menu  is  identified  by  a  small  apple  in  the  top  left  corner.  It  is 
used  primarily  for  desk  accessories.  It  is  also  used  by  most  programs  as  a  way  of 
offering  program  information.  This  is  usually  the  first  item  of  the  Apple  menu. 

There  are  several  accessories  for  some  applications.  For  example,  the 
Mac's  clipboard  and  scrapbook  are  used  to  interact  with  the  word  processing 
program.  It  is  highly  recommended  that  the  user  read  the  Macintosh  Users  Guide 
as  it  explains  the  use  of  the  Mac's  Accessories  which  will  be  of  use  but  will  not  be 
discussed  here. 

2.  File  Menu 

MacRootLocus  starts  with  the  file  menu  first.  It  offers  EQ  Parameter, 
Get  Coeff,  Print  Screen,  Print  Window  and  Quit  commands.  These  selections  also 
have  a  keyboard  shortcut  by  holding  down  the  command  key,  which  has  a  clover 
leaf  symbol  on  it,  and  hitting  the  first  letter  of  the  menu  item  at  the  same  time. 

a.  EQ  Parameter 

The  EQ  Parameter  stands  for  Equation  Parameter.  As  mentioned 
the  previous  chapter,  the  MacRootLocus  program  calculates  the  roots  of  the 
characteristic  equation  for  plotting  points.  It  is  neccessary  to  get  the  degree  of 
polynomials  and  some  equation  parameters  shown  in  Figure  2.3  since  the  Laguerre 
algorithm  [Ref.  2]  is  used  in  order  to  calculate  the  roots. 

The  degree  of  the  polynomial  should  be  between  1  and  10  in 
.MacRootLocus.  Then  there  are  the  default  values  for  the  other  parameters.  These 
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values  avoid  the  convergence  error  for  almost  all  polynomials.  But  if  convergence 
error  messages  appear  on  the  screen,  you  can  change  these  parameter  values.  These 
parameters  must  satisfy  the  following  conditions: 

(1)  Initial  guess  >  0 

(2)  Maximum  Iteration  >  100 

(3)  Tolerance  >  0 


Characteristic  Equation  Parameter 

Degree  of  the  polynomial 

O 

InitGuess 

I 

MaHiter 

Tolerance 

100 

lE-6 

OK 

V  -- 

••  s 

Cancel 

J 

Figure  2.2  Equation  Parameter  Dialog 
b.  Get  Coeff 

After  getting  the  degree  of  the  polynomial,  you  should  insert  the 
coefficients  of  the  polynomial.  The  'Get  Coeff  stands  for  Get  Coefficients.  Ilic 
dialog  for  the  coefficients  i.s  shown  in  Figure  2.3.  This  dialog  is  varied  depending  on 
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the  degree  of  the  polynomial  as  shown  in  Figures  2.3  and  2.4. 

The  algebraic  expression  for  the  coefficients  of  the  characteristic 
equation  of  the  system  may  have  up  to  two  undetermined  parameters  (A  and  B).  In 
the  case  of  the  one— parameter  root  locus  method,  you  use  only  one  undetermined 
parameter  (A).  The  routine  uses  standard  algebraic,  or  infix,  notation  with 
parenthesis  allowed.  Operators  can  include  +,  — ,  *,  /,  and  “  (exponentiation).  The 
unary  minus  sign  is  allowed.  For  example,  the  characteristic  equation  is 
+  (10  +  A)S2  +  (10  *  A  +  5000  *  B)S  +  5000  *  A  =  0 
Figure  2.3  shows  how  you  insert  these  coefficients. 


Figure  2.3  Characteristic  Equation  Coefficient  Data  Dialog 
Box  for  Third  Degree  Polynomial 

If  you  choose  the  Get  Coeff  command  without  the  degree  of  the  polynomial. 
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the  message  shown  in  Figure  2.5  appears  on  the  screen.  This  message  tells  you  that 
degree  of  the  polynomial  has  not  yet  been  entered. 


Figure  2.4  Characteristic  Equation  Coefficient  Data  Dialog 
Box  for  Seven  Degree  of  Polynomial 


Figure  2.5  Message  (1) 
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c.  Print  Screen  and  Print  Window 

These  comands  allow  the  user  to  get  a  hard  copy  of  any  plot 
displayed  by  MacRootLocus.  When  you  choose  the  'Print  Screen'  command,  the 
whole  contents  of  the  screen  is  printed.  If  you  choose  the  'Print  Window'  command, 
the  plot  that  is  to  be  printed  should  be  on  the  active  window  of  the  display. 

Before  attempting  this  method,  ensure  the  printer  is  properly  set  up 
for  friction  feed. 

d.  Quit 

The  last  item  under  the  File  menu  is  Quit.  Selection  of  this  will 
cause  you  to  leave  MacRootLocus  and  return  to  the  desk  top. 

3.  Edit  Menu 

Since  these  operations  are  not  actually  used  in  MacRootLocus,  no  further 
explanation  of  the  Edit  menu  will  be  presented  here.  Any  additional  information 
regarding  the  Edit  menu  can  be  found  in  the  Macintosh  Users  Manual. 

4.  Plot  Menu 

There  are  two  commands.  One  Parameter  and  Two  Parameter.  These  are 
called  to  display  the  dialog  boxes  of  Figures  2.6  and  2.8  to  insert  the  plot  data. 

These  commands  are  followed  after  selecting  the  EQ  Parameter  and  Get 
Coeff  commands.  If  not,  the  message  shown  in  Figure  2.7  appears  to  tell  you  that 
the  degree  of  polynomial  and  polynomial  coefficients  should  be  entered  before  you 
choose  the  plot  menu.  After  reading  the  message,  just  click  once  and  this  box  will 
disappear. 

a.  One  Parameter 

When  you  choose  the  'One  Parameter'  command,  the  dialog  box 
shown  in  Figure  2.6  for  plotting  data  of  the  one— parameter  root  locus  method 
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One  Parameter  Root  Locus  Plot  Data 
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Points  To  Plot  50 


Figure  2.6  One  Parameter  Root  Locus  Plot  Data  Dialog 


appears.  The  plot  default  values  are  shown  in  Figure  2.6.  They  can  be  changed  as 
desired. 

First,  the  user  enters  the  minimum  and  the  maximum  gain  values 
into  the  'Min  Gain'  and  the  'Max  Gain'  insertion  box.  Next,  the  user  selects  one  of 
two  types  of  interval.  Linear  and  Logarithmic.  When  you  click  the  radio  button, 
the  desired  type  of  interval  is  chosen.  For  the  'Linear'  interval,  the  gain  step  size  is 
calculated  by  subtracting  the  minimum  gain  from  the  maximum  gain  entered  in  the 
dialog  box,  and  then  dividing  by  the  number  of  points  to  plot.  Using  'Logarithmic' 
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interval  can  best  be  described  as  giving  equally  spaced  intervals  on  a  logarithmic 
scale. 

Most  MacRootLocus  programs  calculate  the  gain  intervals  using  the 
'Linear'  interval.  This  emphasizes  gain  values  that  are  closer  to  the  max  gain.  This 
becomes  more  evident  as  the  maximum  gain  to  the  minimum  gain  ratio  increases. 
Using  the  'Logarithmic'  interval  gives  more  emphasis  to  the  lower  gains  so  more 
continuous  loci  can  be  drawn.  As  a  basic  rule  of  thumb,  if  the  maximum  gain  to 
minimum  gain  ratio  is  greater  than  100,  selecting  'Logarithmic'  interval  will  give  a 
more  continuous  plot. 


There  is  no  Initial  Degree  or 

Characteristic  Equation  Coefficient. 


Figure  2.7  Message  Box  (2) 

Next,  the  scale  for  the  axis  will  be  chosen.  For  the  'Auto  Scale'  the 
system  calculates  the  minmum  and  maximum  value  of  each  axis.  When  the 
'Manual  Scale'  is  chosen,  the  user  should  insert  the  minimum  and  maximum  values 
for  each  axis. 

The  last  item,  'Points  to  Plot'  sets  the  plot  resolution.  The  bigger 
the  value  you  choose,  the  better  resolution  plot  you  get,  but  the  calculation  time 
will  be  longer. 
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Tluo  Parameter  Root  Locus  Plot  Data 
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Figure  2.8  Two  Parameter  Root  Locus  Plot  Data  Dialog 


6.  Two  Parameter 

The  'Two  Parameter'  command  calls  the  two— parameter  plot  data 
dialog.  It  is  shown  in  Figure  2.8.  The  items  shown  in  Fig^’re  2.8  are  similar  to 
those  shown  in  Figure  2.6,  but  several  items  are  different. 

There  exists  one  more  undetermined  parameter  'B'  to  be  inserted. 
The  'How  many  loci'  item  lets  you  decide  how  many  loci  are  to  be  drawn  for  each 
parameter.  The  number  of  loci  will  be  from  1  up  to  10  for  each  parameter.  In 
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Figure  2.8,  this  value  is  5.  There  is  no  auto-scale  for  axes.  Only  the  manual  scale 
is  available.  You  focus  on  the  interesting  area  for  your  design.  The  last  item  is  the 
marking  and  justification  in  order  to  draw  the  selected  'A'  and  'B'  values  on  tlic 
plot.  There  are  four  radio  buttons.  Two  buttons  are  chosen  each  time  for  each 
parameter,  one  for  position,  the  other  justification  to  draw.  There  are  sixteen 
combinations  available  for  this  work  as  shown  in  Table  2.1.  Figure  2.10  shows  you 
the  ninth  case  in  Table  2.1. 


Table  2.1  The  Combination  for  Marking  and  Justification 
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5.  Help  menu 

MacRootLocus  supports  an  on-line  help  menu  so  that  the  first— time  user 
can  get  desired  results  without  using  trial  and  error. 

Help  is  the  last  item  in  the  menu  bar.  There  are  the  same  item  names  in 
the  menu  bar.  It  makes  it  easy  to  look  for  the  item  for  which  the  user  wants 
information.  The  contents  of  each  item  are  the  subject  of  section  2.B. 
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6.  Information  Box 

Finally,  MacRootLocus  gives  you  a  convenient  way  to  identify  your  plot. 
It  is  a  small  box  in  which  you  can  type  the  information  you  want  to  memorize.  As 
soon  as  the  plot  has  been  completed,  the  box  appears  at  the  bottom  of  the  plot 
window  automatically. 

When  the  up  direction  key  is  pressed  twice,  the  cursor  comes  out  on  the 
box.  Then  you  can  use  the  keyboard  just  like  a  typewriter.  The  capacity  of  the  box 
is  160  to  200  letters.  If  you  do  not  want  that  box,  just  click  once.  It  will  disappear. 
It  is  shown  in  Figure  2.10. 
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Tuio  Parameter  RootLocus 


Figure  2.10  Sample  Plot  (1) 
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III.  DETAILED  PROCEDURE  MODULE  DESCRIPTION 


This  chapter  will  basically  be  programmer's  notes  covering  the  significiaut 
procedures,  functions  and  libraries  to  be  used  in  MacRootLocus.  MacRootLocus 
follows  the  Macintosh  programming  technique  to  make  the  program  user  friendly. 
It  is  a  fairly  simple  Macintosh  application  that  uses  menus,  windows,  dialog  boxes, 
and  graphics.  Most  procedures  and  functions  called  in  MacRootLocus  were 
developed  in  Turbo  Pascal  verson  1.0  for  Macintosh. 

MacRootLocus  consists  of  one  main  program,  one  resource  file  and  six  units. 
The  main  program  integrates  the  resource  file  and  units,  then  it  shows  you  the 
menu  bar  and  tiie  greeting  message  to  start  the  work.  The  details  for  these  will  now 
be  explained. 


A.  MAIN  PROGRAM 

The  Main  Program  is  named  MacRootLocus. Pas.  This  program  consists  of  a 
main  body  and  some  precedures  to  handle  the  system. 

1.  Main  Body 

The  structure  of  the  main  body  uses  the  concept  of  event  driven 
programming.  It  looks  something  like  this: 


Initialize; 

repeat 

SystemTask; 

if  GetNextEvent (everyEvent, theEvent) 
HandleEvent (theEvent) ; 
until  Finished; 

Cleanup; 


{  set  everything  up  ) 

(  keep  doing  the  following  } 

i  update  desk  accessories  ) 

then  {  if  there's  an  event...  } 
{  . . .then  handle  it  ) 

{  until  user  is  done  1 


The  program  is  set  up  once  with  the  user-defined  routine  'Initialize.'  It 
then  enters  a  loop  that  continues  until  some  condition  (such  as  the  user  selecting 
Quit  in  a  menu)  causes  it  to  set  the  boolean  flag  'Finished'  to  true.  Within  that 
loop,  it  performs  two  major  tasks. 

First,  it  calls  'System  Task'  (a  Toolbox  routine),  which  allows  the  Mac 
operating  system  to  update  any  desk  accessories  that  might  be  in  use.  Second,  it 
calls  'GetNextEvent'  (another  Toolbox  routine)  to  see  if  any  events  have  occurred. 
If  any  have,  the  highest  priority  event  is  returned  in  the  data  structure  'the  Event.' 
The  program  then  passes  the  event  to  'Handle  Event,'  which  is  a  user-defined 
routine  that  handles  all  the  different  events  that  might  occur.  Such  events  include 
key— presses,  selection  of  menu  items,  mouse  clicks,  and  windows  being  opened, 
closed,  or  resized.  When  the  program  is  ready  to  terminate,  it  calls  the 
user-defined  routine  clean  up. 

2.  Handle  Event  Procedure 

When  an  event  occurs,  the  operating  system  creates  an  event  record  and 
puis  it  in  a  queue,  ready  for  you  to  handle.  To  see  if  there  is  one  waiting,  you  call 
'GetNextEvent,'  a  boolean  function  that  returns  true  if  there  is  an  event  there  for 
you.  You  give  it  a  mask  of  the  events  you  are  interested  in;  you  can  use  the 
predefined  mask  'Every  Event'  to  look  at  all  events.  This  event  is  passed  to  'Handle 
Event,'  which  takes  care  of  it  .  'HandleEvent'  is  just  a  case  statement  using  the 
'what'  field  in  'the  Event'  to  determine  which  of  the  procedures  to  call. 

3.  DoMouseDown  Procedure 

The  routine  'DoMouseDown'  determines  which  window  the  mouse  was  in  when 
the  clicking  took  place  and  where  exactly  it  happened.  Like  'HandleEvent,' 
'DoMouseDown'  is  mostly  a  case  statement. 
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4.  DoUpdate  Procedure 

The  Macintosh  keeps  track  of  a  lot  of  things  for  you.  For  one,  it  tells  you 
when  some  portion  of  a  window  needs  to  be  drawn,  because  of  resizing  or  removing  a 
covering  window.  This  is  known  as  an  update  event,  and  it  requires  special 
handling.  To  handle  an  update  event,  the  routine  'DoUpdate'  saves  the  current 
'grafport'  into  'SavePort'  and  makes  'the  Window'  the  current  port  so  that  you  can 
write  to  it.  'ReginUpdate'  limits  all  output  to  the  section  of  'theWindow'  that 
needs  updating.  You  then  do  whatever  redrawing  is  needed.  When  you  are  done, 
'EndUpdate'  lifts  those  limits,  and  'SetPort(SavePort)'  restores  the  old  'grafport'. 

5.  Dokeypress  Procedure 

This  routine  handles  the  'Key  down'  and  the  'Autokey'  events.  It  is  a 
check  to  see  if  a  command-key  combination  was  pressed;  if  so,  it  checks  if  the  key  is 
a  menu  command  and  takes  appropriate  action. 

6.  Handle  Menu 

The  procedure  'Handle  Menu'  decodes  the  mouse  position  and  figures  out 
which  menu  and  which  item  in  that  menu  were  selected.  It  uses  a  case  statement  to 
select  the  action  for  the  appropriate  menu;  the  menu  value  is  the  ID  assigned  when 
the  menu  is  created.  The  commands  in  a  menu  are  numbered  from  the  top  down, 
with  the  first  command  having  a  value  of  one.  The  action  itself  is  usually  a  second 
case  statement,  based  on  the  menu  item. 

When  an  item  in  a  menu  is  selected,  the  name  of  that  menu  (in  the  menu 
bar  at  the  top  of  the  screen)  is  highlighted,  that  is,  inverted  to  white-on-black. 
When  you  are  done  processing  the  menu  command,  the  menu  bar  is  restored  to 
normal  by  calling  'HiliteMenu(O),'  which  is  at  the  bottom  of  'HandleMenu.' 
Another  procedure  is  called  before  you  can  leave  'HandleMenu,'  'UpdateMenu.'  a 
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local  procedure  that  tests  to  see  if  certain  items  are  to  be  enabled  or  disabled.  To 
enable  and  disable  menu  items,  the  standard  Macintosh  procedures  'Enableltem' 
and  'Disableltem'  are  called. 

7.  Initialization 

The  initialization  procedure  for  the  MacRootLocus  program  includes  the 
following  structure:  Call  'Init'  routines,  set  up  menus,  set  up  windows,  do  other 
graphic  initialization  and  do  program— specific  initialization. 

There  is  an  'Init'  routine  for  most  of  the  major  managers.  The  first,  and 
most  important,  is  'InitGraf  (thePort).  That  sets  up  'QuickDraw'  (which  is  used 
by  just  about  everything  else)  and  sets  up  a  'grafport'  for  the  screen.  Other  'Init' 
routines  are:  InitFonts,  InitWindow,  InitMenus,  TEInit,  and  InitDialogs  (NIL). 

Setting  up  menus  involves  four  steps.  First,  it  defines  the  menus 
themselves.  If  a  resource  file  is  used,  just  do  a  call  to  'GetMenu'  for  each  menu 
handle,  or  even  a  single  call  to  'GetNewMBar.'  Otherwise,  it  has  to  build  each 
menu  using  an  initial  call  to  'NewMenu',  followed  by  a  call  or  calls  to  Append 
Menu.  Second,  if  it  is  handling  desk  accessories,  call  'AddResMenu'.  Third,  add  all 
the  menus  to  the  menu  bar  by  marking  successive  calls  to  'InsertMenu.'  'Finally, 
call  'DrawMenuBar'  to  display  the  menu  titles  and  make  them  active. 

As  with  menus,  the  window  initialization  takes  several  steps.  If 
MacRootLocus  needs  a  window  at  start  up,  create  it  using  either'GetNewwindow' 
(reading  in  from  resources)  or  Newwindow  (building  it  in  place).  Having  created 
the  window,  make  it  the  current  'grafport'  by  calling  'SetPort',  then  make  it  the 
active  window  by  calling  'Select Window.' 

The  program-specific  initialization  should  probably  come  here.  All  of 
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the  default  values  for  the  dialogs  are  defined  and  the  array  vectors  are  initialized 
here. 

B.  GLOBALVAR  UNIT 

This  unit  declares  all  of  the  whole  variables  to  be  used  in  MacRootLocus 
except  a  few  of  local  variables.  This  unit  is  called  by  the  main  program  and  all 
units.  It  defines  constant,  data  type,  and  variables. 

C.  MAKEROOT  UNIT 

The  unit  MakeRoot  is  very  important  in  MacRootLocus.  This  unit  provides 
several  procedures  and  functions  to  find  the  roots  and  the  array  vectors  for  plot. 
Since  the  characteristic  equation  is  derived,  the  program  must  be  able  to  parse  the 
user's  characteristic  polynomial  coefficient  equations  in  order  to  understand  the 
relations  and  be  able  to  iteratively  substitute  in  values  for  undetermined 
parameters;  A  (for  one-parameter  root  locus  method)  or  A  and  B  (for 
two— parameter  root  locus  method). 

This  unit  has  two  main  procedures,  the  Get  Root  1  procedure  for  the 
one— parameter  root  locus  method  and  the  Get  Root  2  procedure  for  two-parameter 
root  locus  method.  The  simplified  algorithm  for  the  two-parameter  root  locus 
method  is  outlined  in  Figure  3.1. 

There  are  several  procedures  to  perform  this  algorithm.  These  will  be 
explained  in  the  following  section  except  for  the  Rootfinder  unit.  The  InfixtoPoIish 
and  the  ComputePolish  are  especially  interesting. 

1.  InfixtoPoIish  and  ComputePolish  Procedure 

The  coefficients  of  the  polynomial  equation  are  written  in  algebraic,  or 
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Given  a  system's  characteristic  polynomial  : 


C.E  —  a  S°  +  a  S'*  ^  +  ....  +  aS  +  a  —  0 

n  n-l  1  0 

where  an,  an-i,  etc.  are  algebraic  expressions  in  A  and  B.  Also  A 
is  to  be  stepped  from  the  minimum  A  to  the  maximum  A  and  B  varied 
from  the  minimum  B  to  the  maximum  B  and  the  reverse  is  processed. 

Then 

SET  A  =  Min  A  .  SET  B  =  Min  B. 

SET  A  DeltaStep  =  abs((A  Max  -  A  Min)  /  Step) 

SET  B  DeltaStep  =  abs((B  Max  -  B  Min)  /  Step) 

FOR  1  =  1  to  2  do  {  case  of  A  and  B  parameter) 

FOR  j  =  1  to  Step  do  {Step  =  quantity  of  losi) 

IF  Step  A  then  A  =  A  Min  +  A  DeltaStep  *  (j  -  1) 

ELSE  B  =  B  Min  +  B  DeltaStep  *  (j  -  1) 

WHILE  point  no  <=  (points  -  1)  do 
IF  Step  A  then  B  =  NextGain2 
ELSE  A  =  NextGain2 
FOR  term  =  Initial  Degree  downto  0  do 
CONVERT  ai  from  Infix  to  Polish 
SUBSTITUTE  values  for  A  and  B 
COMPUTE  ai 
END  (FOR) 

CALL  RootFinder 

CALL  Results  (make  plot  array  and  numerical  data) 
ENDfWHILE) 

CALL  Plot  Root  LOcus2 
END{FOR} 

Step  A  =  False  (Next  case) 

ENDfFOR) 


Figure  3.1  Two  Parameter  Root  Locus  Algorithm 


'infix'  notation.  The  available  operators  include  (+)  addition,  (-)  subtraction,  (*) 
multiplication,  (/)  division,  and  (“)  exponentiation.  These  operators  follow  a 
hierarchical  precedence  with  exponentiation  operations  being  done  first,  followed  by 
multiplication  and  division,  and  finally  addition  and  subtraction.  Operations  like 
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multiplication  and  division  which  have  the  same  precedence  are  performed  from  left 
to  right  when  conflicts  arise.  To  change  the  order  of  precedence,  parentheses  may 
be  used  around  any  set  of  operations.  These  parenthetical  expressions  have  the 
highest  priority  and,  when  nested,  the  innermost  operations  within  parentheses  are 
done  first.  This  scheme  follows  closely  the  protocol  used  in  most  calculators  and 
high  level  programming  languages. 

Infix  notation,  while  convenient  for  the  program  user,  does  not  lend  itself 
well  to  computer  manipulation.  A  better  way  to  represent  equations  for  the 
computer  is  the  so  called  'reverse  Polish  notation.'  In  reverse  Polish  notation,  the 
operands  of  an  equation  are  entered  first,  followed  by  the  operator.  For  example, 
the  infix  expression 

3’''4  +  5 


would  be  represented  as 

34*5  + 

in  reverse  Polish  notation.  The  numbers  3  and  4  are  entered  and  multiplied,  then  5 
is  entered  and  added  to  the  previous  result.  Using  the  concept  of  a  'stack'  the 
reverse  Polish  expression  is  easy  to  evaluate. 

Recall  that  a  stack  is  a  last— in— first— out  queue  whose  operation  is 
analogous  to  a  stack  of  trays.  To  operate  the  stack,  the  program  calls  a  'push' 
procedure  to  place  an  item  on  the  stack,  and  a  'pop'  procedure  to  remove  the  top 
item.  Now,  using  the  example  given  above,  an  arithmetic  evaluation  procedure  can 
be  illustrated.  Figure  3.2  demonstrates  such  an  implementation. 

The  basic  equation  evaluation  algorithm  can  be  outlined  in  three  steps: 

( 1 )  Scan  the  reverse  Polish  equation  term— by— term. 

(2)  If  the  term  is  a  constant  then  push  it  onto  the  stack. 
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3 

4 

(4.  3  *)=  12 

5 

(5,  12  +)-  17 

17 

<-stack-> 

3 

{pop4} 

{pop3} 

{multiply} 

12 

{pop5} 

{pop12} 

{add} 

{push  3}  {push  4}  {push  121  (push  17} 

^  {push  5} 


Figure  3.2  Example  of  the  Stack  Operation 

(3)  If  the  term  is  an  operator  then  pop  the  first  two  items  off  the 

stack,  apply  the  operator,  and  push  the  result  back  onto  tlie 
top  of  the  stack. 

When  the  algorithm  is  completed,  the  answer  to  the  expression  will  be  on  the  top  of 
the  stack. 

To  get  the  infix  equation  into  reverse  Polish  form  is  a  bit  more  difficult 
than  simply  evaluating  the  Polish  expression.  Of  special  consideration  when 
building  the  Polish  form  of  the  equation  are  the  operator  priorities  and  the  use  of 
parenthesis  to  change  those  priorities.  A  set  of  rules  can  be  written  which  outline 
the  conversion  procedure  (Ref.  2].  These  rules  are  discussed  below  with  an 
illustrative  eaxmple.  The  infix  expression 

8 +  (7-6/ 3)  *2 

will  be  parsed  using  the  following  operator  priority  table. 

Operator _ Priority 

4 

*,/  3 

+,-  2 

operand  1 

(,), space  0 
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Using  a  result  string  called  RPN,  and  an  operator  stack  the  rules  for  infix  to  reverse 
Polish  conversion  follow: 

(1)  If  an  operand  is  encountered,  move  it  to  RPN. 

(2)  If  an  operator  is  encountered,  move  all  higher  priority 
operators  on  the  stack  to  RPN  and  push  the  new  operator 
onto  the  stack. 

(3)  If  a  left  parenthesis  is  encountered  push  it  onto  the  stack. 

(4)  If  a  right  parenthesis  is  encountered,  pop  all  operators  off  the 
stack  and  append  them  to  RPN  until  a  left  parenthesis  is 
encountered.  Discard  both  parenthesis. 

(5)  When  finished  with  the  infix  expression,  pop  all  remaining 
operators  from  the  stack  and  append  them  onto  RPN. 

Applying  these  rules  to  the  example  problem  above  is  shown  in  Figure  3.3. 


PRN 

STACK 

INFIX 

RULE 

8 

8+(7-6/3)*2 

+(7-6/3)*2 

1 

8 

+ 

(7-6/3)*2 

2 

8 

+  ( 

7-6/3)‘2 

3 

8  7 

+  ( 

-  6  /  3  )  *  2 

1 

8  7 

+  (  - 

6  /  3  )  *  2 

2 

8  7  6 

+  (  - 

/  3  )  *  2 

1 

8  7  6 

+  (  -  / 

3  )  *  2 

2 

8  7  6  3 

+  (  -  / 

)*2 

1 

8763/- 

+ 

*  2 

4 

8  7  6  3  /- 

+  * 

2 

2 

8  7  6  3  /  -  2 

+  * 

1 

8  7  6  3  /  -  2  *  + 

Figure  3.3  Conversion  from  Infix  to  Reverse  Polish 
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2.  NextGainl  Function 

This  function  makes  two  kind  of  intervals,  linear  interval  and 
logarithmic  interval.  This  function  is  called  by  the  'GetRootT  procedure. 

3.  NextGiun2  Function 

This  is  like  the  'NextGainl'  function  except  it  has  two  parameters.  It 
supports  two  kind  of  intervals  for  each  parameter.  This  function  is  called  by  the 
'GetRoot2'  procedure. 

4.  Results  Procedure 

This  procedure  outputs  the  calculated  roots  to  the  device  'Out  File'. 
Then  it  makes  it  possible  for  the  user  to  access  the  numerical  data  from  the  hard 
disk.  Also,  it  supports  the  plot  array  named  'GraphArray'  for  the  plotting  data. 

5.  PlotRootLocusl  and  PlotRootLocus2  Procedure 

The  structure  of  these  procedures  is  exactly  the  same  except  for  the 
auto-scale  function  of  the  coordinate  to  be  supported  by  the  'PlotRootLocusl.' 
These  procedures  draw  the  roots  and  some  information  in  the  desired  coordinate. 
First,  the  'SelectWind'  is  called  with  the  integer  number  and  a  boolean  expression  in 
order  to  select  a  window  as  visible.  The  'Define  Header'  is  called  to  draw  the  title. 

Next  the  'OpenPic'  opens  a  picture  for  a  specific  window  and  only  shows 
the  drawing  if  the  boolean  is  set  True.  In  order  to  defind  the  coordinate  it  calls 
'FindWorld'  for  the  auto-scale  or  'FindWorldl'  for  the  manual— scale.  Then 
'DrawAxis'  draws  an  axis  with  Footers  and  optional  arrows  on  the  axis.  Finally, 
'DrawPolygon'  draws  a  polygon  defined  in  the  plot  array  with  the  predefined  shape. 

D.  ROOTFINDER  UNIT 

The  unit  which  solves  for  the  roots  of  a  polynomial  is  called  'RootFinder.' 
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This  unit  uses  Laguerre's  method  with  linear  deflation.  Since  this  method  is  a 
commercially  available  package  of  subroutines  in  the  'Turbo  Pascal  ToolBox 
(Numerical  Methods)',  a  brief  explanation  is  offered  here. 

1.  Laguerre's  Method 

Laguerre's  method  attempts  to  approximate  all  the  real  and  complex 
roots  of  a  real  or  complex  polynomial.  Laguerre's  method  is  very  reliable  and  quick, 
even  when  converging  to  a  multiple  root. 

To  motivate  (although  not  rigorously  derive)  the  Laguerre  formulas  we 
can  note  the  following  relations  between  the  polynomial  and  its  roots  and 
derivatives. 


Pn(x)  =  (x  -  Xt)(x  -  X2)  .  .  .  (x  -  Xn) 


(3.1) 


ln(Pn(x)(  =  ln|x  -  xi|  Injx  -  X2I  + 
+  ln|x-xn| 


(3.2) 


din 


|Pn(x)|  _  1  I  2 

OX  X  -  Xi  X  -  X2 

I  ^  _  P  n  -  p 

+  x  -  x„  -  -p;-  -  ^ 


(3.3) 


+ 


1 


d^ln|Pn(x)|  ^  1 

dx^  (X  -  Xi)^  (X  -  X2)^ 


+ 


1 


(x  -  Xn)^ 


Pn 

T7 


^=ti 


(3.4) 


Starting  from  these  relations,  the  Laguerre  formulas  make  what  Acton  nicely  calls 
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"a  rather  drastic  set  of  assumptions",  the  root  xl  that  we  seek  is  assumed  to  be 
located  some  distance  a  from  our  current  guess  x,  while  all  other  roots  are  assumed 
to  be  located  at  a  distance  b. 


a  =  x  — Xl  b  =  x  — xi  i  =  2,  3,.  .  n 


(3.5) 


Then  we  can  express  Eqs  (3.3),  and  (3.4)  as 


(3.6) 


=  H 


which  yields  as  the  solution  for  a 


(3.7) 


a  = - -  -  -  “  -  (3.8) 

G  ±  '^(n  -  l)(nH  -  G^) 

where  the  sign  should  be  taken  to  yield  the  largest  magnitude  for  the  denominator. 
Since  the  factor  inside  the  square  root  can  be  negative,  a  can  be  complex.  (A  more 
rigorous  justification  of  Eq.  (3.8)  is  in  [Ref.  3].) 

The  method  operates  iteratively;  for  a  trial  value  x,  a  is  calculated  by 
Eq.  (3.8).  Then  x  -  a  becomes  the  next  trial  value.  This  continues  until  a  is 
sufficiently  small.  In  the  next  section,  a  major  procedure  which  handles  this  method 
will  be  explained. 
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2.  InitAndTest  Procedure 

This  procedure  sets  the  initial  value  of  the  input  and  output  variables. 
This  procedure  also  tests  the  tolerance  (Tol),  maximum  number  of  iterations 
(Maxiter),  and  code.  Finally,  it  examines  the  coefficients  of  Poly.  If  the  constant 
term  is  zero,  then  zero  is  one  of  the  roots  and  the  polynomial  is  deflated  accordingly. 
Also  if  the  leading  coefficient  is  zero,  the  degree  is  reduced  until  the  leading 
coefficient  is  non-^reo. 

3.  FindOneRoot  Procedure 

This  procedure  approximates  a  single  root  of  the  polynomial  Poly.  The 
root  must  be  approximated  within  Maxiter  iterations  to  a  toleran^^’  of  Tol.  The 
root,  a  value  of  the  polynomial  at  the  root,  and  the  number  of  iterations(Iter)  are 
returned.  If  no  root  is  found,  the  appropriate  error  code  (Error)  is  returned. 

4.  EvaluatePoly  Procedure 

This  procedure  applies  the  technique  of  synthetic  division  to  determine 
the  value  (yValue),  first  derivative  (yPrime)  and  second  derivative  (yDoublePrime) 
of  the  polynomial,  Poly,  at  X.  The  0th  element  of  the  first  synthetic  division  is  the 
value  of  Poly  at  X,  the  1st  element  of  the  second  synthetic  division  is  the  first 
derivative  of  Poly  at  X,  and  twice  the  2nd  element  of  the  third  synthetic  division  is 
the  second  derivative  of  Poly  at  X. 

5.  ConstructDifference  Procedure 

This  procedure  computes  the  difference  between  approximations;  given 
information  about  the  function  and  its  first  two  derivatives. 

6.  TestForRoot  Function 

These  are  the  stopping  criteria.  Four  different  ones  are  provided.  If  you 
wish  to  change  the  active  criteria,  simply  comment  off  the  current  criteria 


35 


(including  the  appropriate  OR)  and  remove  the  comment  brackets  from  the  criteria 
(including  the  appropriate  OR)  you  wish  to  be  active. 

7.  ReducePoly  Procedure 

This  procedure  deflates  the  jwlynomial  Poly  by  factoring  out  the  Root. 
Degree  is  reduced  by  one. 

E.  MESSAGE  UNIT 

This  unit  provides  the  several  messages  which  inform  the  user  with  some 
warnings  and  help  informations.  There  are  11  procedures, which  have  the  same 
structure,  to  provide  the  content  of  a  message  and  one  procedure,  'SetupWindow,' 
to  define  the  window.  Here,  the  'SetupWindow'  procedure  and  one  sample 
procedure  to  make  message  will  be  explained. 

1.  SetupWindow  Procedure 

This  procedure  defines  all  of  the  windows  in  the  program.  It  calls  the 
'Define Window'  procedure  from  'TurboGraph'  Unit.  There  are  five  other  standard 
window  types:  documentPrc,  DocProc,  dBoxProc,  PlainDBox,  and  noGrowDocProc 
[Ref.  4:  P.  463].  We  can  insert  window  type,  window  ID,  and  size  of  window  in  this 
procedure. 

2.  MakeInfoScreen  Procedure 

This  is  one  of  nine  procedures  for  messages.  It  is  a  sample  to  explain 
this  kinds  of  procedure.  'MakeInfoScreen'  procedure  brings  up  a  window  with  a 
description  of  what  this  program  is.  First,  a  'GrafPort'  is  called.  It  is  simply  a 
Pascal  record  type  with  fields  that  control  QuickDraw's  behavior. [Ref.  4;  p.  405  — 
422]  'GrafPort'  allows  an  application  with  windows  to  invoke  drawing  operations 
that  are  appropriate  for  each  window.  The  'SelectWind'  chooses  the  window  to  be 
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defined  in  'SetupWindow'  using  just  window  ID  for  this  procedure.  The  selected 
window  moves  to  the  center  of  the  screen  using  the  'Movewindow'  and  the 
'SetVisibility'  sets  the  visibility  of  a  window.  The  'TextFont,'  'TextSize,'  and 
'TextStyle'  are  used  to  define  the  letter  to  be  drawn  in  the  window.  Finally  the 
'MoveTo'  assigns  the  location  to  be  drawn.  'DrawString'  draws  the  messages  in  the 
window. 

F.  MYDIALOG  UNIT 

This  unit  supports  4  dialog  boxes  for  input  data.  Each  dialog  procedure  has 
the  same  skeleton  to  make  a  program;  the  only  major  procedure  to  construct  the 
program  will  be  explained  in  this  section. 

A  dialog  box  communicates  with  the  background  text,  controls,  icons,  and 
'editText'  items  in  which  the  user  can  enter  and  edit  text;  the  user  talks  back  with 
the  mouse  and  keyboard. 

The  Dialog  Manager  [Ref.  5:  p.  53]  handles  this  communication.  The  Dialog 
manager  leans  heavily  on  the  resource  mechanism  [Ref.  5:  pp.  137  -  154]  to  provide 
it  with  data  structures.  Basically,  each  dialog  is  represented  on  disk  as  a  resource  of 
'resType'  DLOG-a  template  describing  the  dialog's  size,  window  type,  and 
title-and  a  resource  of  type  DITL  (Dialog  Item  List),  which  lists  the  content  of  the 
dialog  (controls  background  text,  and  so  on). 

A  program  that  intends  to  use  dialogs  sees  to  it  that  appropriate  resources  of 
type  DLOG  are  available  at  runtime.  Then,  the  program  loads  them  into  memory 
and  draws  them  on  the  screen  with  'GetNewDialog.' 

After  'GetNewDialog'  has  read  a  dialog  template  and  its  item  list  into  memory 
and  drawn  it  on  the  screen,  a  program  enters  a  loop  in  which  it  repeatedly 
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calls  the  lynchpin  of  dialog  processing,  'ModalDialog,'  and  acts  on  the  integer  value 
it  returns.  'ModalDialog'  allows  the  user  to  customize  its  operation  with  a  routine 
specified  by  the  user.  The  user  informs  'ModalDialog'  of  this  by  passing  a  specific 
parameter  that  points  to  the  routine,  which  then  gets  control  every  time  an  event  is 
generated  when  the  dialog  is  on  the  screen.  It  is  up  to  the  filter  routine  to  decide 
what  to  do  about  it. 

To  set  the  value  of  a  control,  you  need  to  call  'SetCtlValue'.  Since  the 
'SetCtlValue'  expects  a  handle  to  the  control  record,  you  will  have  to  get  a  handle 
to  the  control  that  needs  setting.  This  is  a  task  for  the  Dialog  Manager's 
'GetDItem'  routine.  'GetDItem'  takes  in  two  bits  of  information  and  returns  three. 
Given  the  indicated  dialog  and  item  number,  it  returns  the  information  about  that 
particular  item;  its  type  (such  as  'radioButton'  or  'staticText',  encoded  as  an 
integer),  a  handle  to  its  underlying  data  structure  (for  controls,  this  is  a 
'ControlHandle'),  and  finally,  its  bounding  box.  The  key  working  with  'editText' 
items  is  to  use  'GetIText'  procedures.  An  'editText'  item  begins  life  with  the 
starting  content  assigned  by  its  definition  in  an  RMaker  file  [Ref.  5:  p.  8].  After  the 
user  has  played  with  it  (as  determined  by  a  suitable  'itemHit'  value  returned  by 
'ModalDialog'),  'GetIText'  is  used  to  see  what  the  value  is  now. 

G.  TURBOGRAPH  UNIT 

This  unit  supports  many  procedures  for  graphics  under  the  Turbo  Pascal 
environment.  It  has  commercially  available  packages  in  the  Turbo  Pascal  Toolbox 
(Numerical  Method). 

Many  procedures  and  functions  in  this  unit  are  called  by  several  units.  Since 


38 


comments  in  source  code  of  this  unit  give  you  enough  information. 
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IV.  EXAMPLE  OF  DESIGN 


A.  OVERVIEW 

The  root  locus  method  is  very  valuable  in  the  analysis  of  dynamic  systems, 
and  is  also  used  for  design.  By  inspection  of  the  curve,  we  determine: 

(1)  Whether  any  loci  cross  the  jo;  axis  into  the  right  half  s— plane  ( If  such  a 
crossing  exists,  it  defines  a  stability  limit  for  the  system.), 

(2)  The  frequency  (value  of  oi)  at  such  a  stability  limit, 

(3)  Whether  the  loci  go  through  any  area  on  the  s— plane  where  we  want 
dominant  complex  roots  for  our  system  (i.e.,  can  we  get  what  we  want). 

We  can  also  determine  the  value  of  the  parameter  at  the  stability  limit,  the 
value  of  the  roots  for  any  specified  value  of  the  parameter,  and  the  value  of  the 
parameter  required  to  place  roots  at  any  selected  point  on  the  loci.  These  latter 
items,  however,  are  not  done  by  inspection  of  the  curves,  but  by  inspection  of  the 
computer  printout  or  by  separate  calculations.  On  the  MacRootLocus,  both  plot 
and  calculation  data  are  supported.  Since  two  parameters  of  the  system  are 
adjustable,  it  is  possible  to  calculate  and  plot  a  family  of  root  loci  for  the  pair  of 
parameters  The  technique  for  synthesizing  a  system  utilizing  the  rootlocus  method 
is  demonstrated  in  this  chapter. 

B.  GRAPHICAL  SOLUTION 

In  this  section,  four  different  kinds  of  control  problems  will  be  demonstrated. 
Such  problems  include  simple  cascade  compensators  for  systems,  subject  only  to 
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step  inputs  and/or  load  disturbances,  and  feedback  compensators  with  no  more  than 
two  adjustable  variables. 

1.  Example  1  (cascade  lead  compensation) 

The  block  diagram  of  Figure  4.1  shows  the  general  case  of  cascade 
compensation. 


Figure  4.1  Simple  Casade  Compensation  Block  Diagram 


Consider  a  plant  with 


Gp(S)  = - t  .1) - .  (4.1) 

S2(S  +  1)(S  +  100) 

The  uncompensated  system  is  unstable.  It  is  desired  to  stabilize  the  system  with  a 
low— pass  filter  to  reduce  bandwidth.  Specifications  are  now  in  the  form  of  a  desired 
location  for  the  roots  of  the  characteristic  equation.  A  cascade  lead  compensator  is 
to  be  used,  and  we  choose  the  simplest  possible  i.e., 
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Also  we  define  A  and  B  as 


Gc(S) 


A  = 


P  (S  +  Z) 

Z  (S  +  P)  ’ 


1 

z 


(4.2) 


(4.3) 


The  characteristic  equation  for  this  system  is 


+  (101  +B)S'‘  +  (100  +  101B)S3  +  (lOOB  +  400A)S2 

+  400A(1  +  B/A)S  +  400B  =  0.  (4.4) 


The  two  parameter  root  locus  family  for  this  system  is  given  on  Fig  4.2.  Since  the 
desired  roots  are 

S  =  1.684  ±j  3.579, 

an  appropriate  choice  from  this  plot  is  A  =  5.0,  B  =  4.8,  from  which  Z  =  0.96,  and 
P  =  5.0. 
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Two  Parameter  RootLocus 


Figure  4.2  Plot  of  Example  1 
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2.  Example  2  (cascade  lag  compensation) 

If  we  consider  a  type-one  third-order  plant  for  which 

Gp(S)  = - ^ - 

S(S  +  1)(S  +  10) 

and 

Gc(S)  =  ■ 

Z  (S  +  P) 

Then,  letting 

A  =  B  =  P, 

Z 

the  characteristic  equation  for  the  cascade  lag  compensated  system  is 

+  (11  +  B)S3  +  (10  +  11B)S2  +  (lOB  +  150A)S 

+  150B  =  0.  (4.8) 

The  two  parameter  root  Locus  for  the  lag  relocation  zone  is  shown  on 
Figure  4.v3.  Choosing  A  =  0.1  and  B  =  0.1  gives  P  =  0.01,  Z  =  0.1  and  provides 
complex  roots  at 

S  =  0.0372  ±  j  1.12, 
with  real  roots  at  S  =  —1.014,  —10.15. 

Stabilization  can  also  be  achieved  with  a  lead  filter,  but  the  resulting 
system  will  have  a  wider  bandwidth  and  faster  response. 


(4.5) 

(4.6) 

(4.7) 
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Figure  4.3  Plot  of  Example  2 
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3.  Example  3  (vdocity  feedback) 

The  basic  problem  is  given  by  the  block  diagram  of  Figure  4.4.  G(s) 
may  be  any  order. 


Figure  4.4  Velocity  Feedback  Compensation  Block  Diagram 


For  a  second-order  system,  let 


G(s)  =  100  /  S{S+2). 


(4.9) 


We  want  the  root  loci  for 


100  AS 

S(S  +  2)  +  100 


(4.10) 


so  the  characteristic  equation  is 
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-I-  (2  +  100A)S  +  100  =  0. 


(4.11) 


The  resulting  loci  are  shown  on  Fig  4.5.  Inspection  of  Fig  4.5  shows  that  increasing 
the  feedback  gain  increases  damping.  Any  desired  value  of  (  is  available  for  the 
complex  roots.  With  high  feedback  gain,  overdamping  (all  real  roots)  is  available, 
and  no  positive  value  of  gain  can  make  the  system  unstable.  To  design  the  system, 
pick  a  root  location  on  the  locus,  use  the  magnitude  rule  to  find  the  gain  as  in 
Figure  4.5,  or,  since  a  computer  program  was  used  to  generate  the  locus,  tlic 
tabulated  data  are  supplied  to  get  the  desired  information.  Part  of  the  tabulated 
data  are  shown  below. 


A  =  0.07971 

-4.98533003622524e+0  +  8 . 66870719484229e+0  j 

-4.98533003622524e+0  +-8 . 66870719484229e+0  j 

A  =  0.08633 

-5.31655371460862e+0  +  8 . 46960781852863e+0  j 

-5.31655371460862e+0  +-8 . 46960781852863e+0  j 

A  =  0.09351 

-5.67530563384649e+0  +  8 . 23352330186964e+0  j 

-5.67530563384649e+0  +-8.23352330186964e+0  j 

A  =  0.10128 

-6.06387368606130e+0  +  7 . 95169390252752e+0  j 

-6.06387368606130e+0  +-7 . 95169390252752e+0  j 

A  =  0.10969 

-6. 484735911754086+0  +  7 . 61237151975697e+0  j 

-6.484735911754086+0  +-7 . 61237151975697e+0  j 

A  =  0.11881 

-6.940576303174546+0  +  7. 19919443964476e+0  j 

-6.940576303174546+0  +-7.199194439644766+0  j 

A  =  0.12869 

-7.434301921120686+0  +  6. 68813538631070e+0  j 

-7.434301921120686+0  +-6.688135386310706+0  j 
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Figure  4.5  Plot  of  Example  3 
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4. 


Example  4  (velocity  and  acceleration  feedback) 

The  block  diagram  of  Figure  4.6  feedback  shows  the  general  case  of 
velocity  and  acceleration  feedback. 


Figure  4.6  Velocity  and  Acceleration  Feedback  Compensation  Block  Diagram 


If  we  consider  a  third  order  plant  for  which 


Gp(s) 


8 

(S  +  1)'  ’ 


(4.12) 


then  the  root  locus  form  is 


8(AS  +  BS^) 
(S  +  1)’  +8 


(4.13) 


and  the  characteristic  equation  is 
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S3  +  (3  +  8A)S2  +  (3  +  8B)S  +  9  =  0. 


(4.14) 


The  two— parameter  root  locus  family  is  given  on  Figure  4.7.  We  may  select 
desired  location  for  the  complex  roots  such  as  s  =  —0.73+  j  1.07,  for  which  Ai 
0.48,  B2  =0.8.  The  real  root  is  at  s  =  -  5.38  so  the  complex  roots  are  dominant. 
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Figure  4.7  Plot  of  Example  4 
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V.  CONCLUSION  AND  RECOMMENDATION 


MacRootLocus  is  a  useful  and  powerful  tool  for  the  designer  and  a  simple  and 
easy  to  use  package  for  the  user.  The  program  supports  the  two  parameter  root 
locus  method  intensively  as  well  as  the  one  parameter  root  locus  method.  Also,  it 
offers  high  resolution  plots,  the  tabulated  data  and  a  small  text  editor  for  recording 
some  information  of  the  design. 

The  program  does,  however,  have  some  areas  where  improvements  are 
possible.  These  are  discussed  below. 

(1)  According  to  the  engineer's  needs,  the  root  finder  can  be  replaced 
easily  with  a  more  reliable  and  faster  algorithm,  and 

(2)  the  capability  to  move  the  information  window  can  be  added. 

Since  this  program  was  developed  under  the  standard  interface  philosophy  the 
Macintosh  was  designed  for,  it  is  a  simple  matter  to  change  and  to  append  some 
subroutines. 

MacRootLocus  is  the  user  friendly  CAD  program  available  that  is  both  simple 
enough  for  the  beginning  student  to  easily  use  as  well  as  powerful  and  flexible 
enough  to  benefit  the  experienced  system  designer. 
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APPENDIX 


SOURCE  CODE 

This  appendix  contains  the  source  code  of  all  modules  in  MacRootLocus 
except  Standard  Apple  Macintosh  libraries.  Some  subroutines  in  Numerical  Method 
(Turbo  Pascal  Tool  Box)  are  also  included  since  they  are  slightly  modified.  A  disk 
is  available  from  Dr.  Thaler  that  contains  the  MacRootLocus  source  code  and  the 
MacRootLocus  resource  file. 
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MacRootLocus 
versioni  .0 


04  DEC  1989 


author 


KO.  SUNG  HOON 


r_ 

1= 

r_ 

r 


SYSTEM 

LANGUAGE 

LIBRARIES 

RESOURCE 


:  Apple  Macintosh  SE 
:  Turbo  Pascal  for  Macintosh  versioni  .0 
:  Numerical  Method  (Turbo  Pascal  Tool  Box) 
;  MacRootLocus.rsrc 


program  MacRootLocus; 

{=  This  main  program  handles  the  system  and  integrate  theresource 
{=  and  6  units. 


{$B+}  {  Set  the  bundle  bit } 

{$R  MacRootLocus.Rsrc}  {  Identify  resource  file  for  menu  and  icon  info  } 
{$T  APPLFFTD}  {  Set  the  application  type  and  creator } 

{$S-t-}  {  Generate  segmented  code} 

{$!-}  {  Turn  off  I/O  error  checking  } 

{$U  SpecVar) 

{$U  RootsFinder} 

{$U  MakeRoot) 

{$U  TurboGraph} 

{$U  Message) 

{$U  MyDialog} 


uses 

MemTypes,  QuickDraw,  OSIntf,  ToolIntf.Packlntf, 
PasPrinter,SANE,MacPrint,  RootsFinder,  SpecVar, 
{$S  Second  Segment) 

TurboGraph, 

Message, 

{$S  Third  Segment) 

MakeRoot, 

MyDialog; 
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function  L  ■)Case(Ch  :  char)  :  char; 


{=  Returns  the  upper  case  of  Ch  -} 


inline 

$301 F,  {  UpCase  MOVE.W  (SP)+.D0  ;  GetCh  } 

$0C40, 

$0061 ,  {  CMP.W  #'a’,D0  ;  skip  if  not  lower  case  } 
$6D0A.  {  BLT.S  @1  } 

$0C40, 

$007A,  {  CMP.W  #’z'.D0  } 

$6E04.  {  BGT.S  @1  } 

$0440. 

$0020.  {  SUB.W  #$20.D0  } 

$3E80:  {  @1  MOVE.W  DO.(SP)  } 

procedure  HideAllWindows; 

{=  Hide  the  three  main  windows  =} 

var 

WindNum  :  integer; 
begin 

for  WindNum  :»  RegendBox  downto  InfoScrnWind  do 
SetVisibility( WindNum.  FALSE); 
end;  {  HideAllWindows } 

procedure  ShowAIIWindows; 

{=  Make  the  three  main  windows  visible  =} 

var 

WindNum  :  integer; 
begin 

for  WindNum  :=  RootLocusWind  to  RegendBox  do 
SelectWind( WindNum.  TRUE); 
end;  {  ShowAIIWindows  } 

procedure  MakeWatchCursor; 
var 

CursorHandle  :  CursHandle; 
begin 

CursorHandle  GetCursor(WatchCursor); 
SetCursor(CursorHandle'''') ; 
end;  {  MakeWatchCursor } 

procedure  MakeArrowCursor; 
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begin 

InitCursor; 

end;  {  MakeArrowCursor } 
procedure  PrintScreen; 


{ . . ~) 

Dump  the  entire  screen  to  the  printer.  -} 


begin 

HideCursor; 
HardCopy(FALSE): 
ShowCursor; 
end;  {  PrintScreen  } 


procedure  DoDeskAcc(ltem  :  Integer); 

{=  start  up  desk  accessory  from  Apple  menu  =} 
var 

SavePort  r  GrafPtr; 

RefNum  ;  integer; 

DName  ;  string; 
begin 

GetPort(SavePort);  {  save  port  before  starting  it  } 

Getltem(MenuLlst(AM],  Item,  DName);  {  get  name  of  desk  accessory  } 
RefNum  OpenDeskAcc(DName);  {  and  start  that  sucker  upl  ) 
SetPort(SavePort);  { restore  grafport  and  continue  } 

end;  { DoDeskAcc } 

procedure  SetltemState(Mndx,lndx  :  Integer;  Flag  :  Boolean); 

{s  if  Flag  is  true,  enables  item  Indx  of  menu  Mndx;  else  disables  =} 
begin 

if  Flag  then 

Enableltem  (MenuList[Mndx],lndx) 
else 

Disabieltem(MenuList[Mndx],lndx); 
end;  {  SetItemState } 

procedure  HandleMenu(Menulnfo  :  Longint); 

{ . - . . . 

{=  Decode  Menuinfo  and  carry  out  command  -) 

var 

Menu  :  Integer;  {  menu  number  that  was  selected  } 
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Item  :  Integer;  { item  in  menu  that  was  selected  } 

WindNum  :  integer; 
begin 

if  Menuinfo  <>  0  then 
begin 

PenNormal;  { set  the  pen  back  to  normal  } 

Menu  :»  HiWord(Menulnfo);  { find  which  menu  the  command  is  in  } 
Item  >  LoWord(Menulnfo);  {  get  the  command  number  } 

case  Menu  of  { and  carry  it  out  } 

AppIMenu  :  if  Item  «  1  then 

DoAbout  {  bring  up  "About..."  window} 
else 

DoDeskAcc(ltem);  { start  desk  accessory  } 

FileMenu  :  case  Item  of  {  File  menu } 

1  :  GetEQParameter; 

2 ;  GetCoeff; 

3  :  PrintScreen;  {Print  screen} 

4  ;  begin 

HideCursor; 

HardCopy(TRUE);  {  Print  top  window } 

ShowCursor; 

end; 

5  ;  Finished  :=  TRUE;  {  Quit  command } 
end; 

EditMenu  :  case  Item  of  {  Edit  Menu } 

1..5  :  if  not  SystemEdit{ltem-1)  then 
{ Do  nothing } 

end; 

PlotMenu  :  case  Item  of  {  bring  up  the  dialog  } 

1  ;  PlotOneParameter;  {  for  one  parameter  plot  data} 

2  :  PlotTwoParameter;  {  for  two  parameter  plot  data} 
end; 

HelpMenu  : 

case  Item  of  {  Help  menu  } 

1  :  InfoGetEQParameter; 

2 ;  InfoGetCoeff; 

3  :  InfoPlotOneParameter; 

4  :  InfoPlotTwoParameter; 

5  :  InfoPrint 
end 

end;{  case } 

HiliteMenu(O);  { reset  menu  bar  } 
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end; 

end;  { HandleMenu } 


procedure  SelectOurWindow(WPtr  :  WindowPtr); 


{=  Select  the  window  pointed  to  by  WPtr  If  it's  one  of  oar's  «} 


begin 

if  (WPtr  -  Window[RootLocusWind].P)  then 
SelectWind(RootLocusWind,  TRUE); 
end;  { SelectOurWindow } 

procedure  HandleClick(WPtr  :  WindowPtr;  MLoc  :  Point); 

{=  Handle  a  mouse  click  within  a  window  =} 
begin 

if  WPtr  <>  FrontWindow  then  { if  it's  not  in  front...  } 
SelectOurWindow(WPtr); 
end;  {  HandleClick } 

procedure  HandleGoA  way  (WPtr  :  WindowPtr;  MLoc  :  Point); 

{=  Handle  a  mouse  click  in  the  go-away  box  of  a  window  =} 
var 

WPeek  :  WindowPeek;  { for  looking  at  windows  } 
begin 

if  WPtr  *  FrontWindow  then  {  if  it's  the  active  window  } 
begin 

WPeek  :=  WindowPeek(WPtr);  {  peek  at  the  window  } 
if  TrackGoAway(WPtr,  MLoc)  then  {  and  the  box  is  clicked  } 
begin 

if  WPeek''. WindowKind  =  userKind  then  { if  it’s  our  window  } 
HideAllWindows  { time  to  stop  } 

else 

CloseDeskAcc(WPeek''.WindowKind)  { close  DeskAcc  } 
end 
end 
else 

SelectOurWindow(WPtr); 
end;  { HandleGoAway } 

procedure  HandleZoom(WRr  :  WindowPtr;  MLoc  :  Point;  WLoc  ;  integer); 


{ . 

Handle  a  mouse  click  in  zoom  box  of  a  window  -} 
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var 

WPeek  :  WindowPeek;  { for  looking  at  windows  } 
begin 

if  WPtr  =  FrontWindow  then  {  if  it's  the  active  window  } 
begin 

WPeek  >  WindowPeek(WRr);  {  peek  at  the  window  } 
if  TrackBox(WPtr,  MLoc,  WLoc)  then  {  and  the  box  is  clicked  } 
begin 

if  WPeek^.WindowKind  -  userKind  then  { if  it's  our  window  } 
ZoomWindow(WPtr,  WLoc,  FALSE); 
end 
end 
else 

SelectOurWindow(WPtr) ; 
end;  { HandleZoom } 

procedure  HandleGrow(WPtr  :  WindowPtr;  MLoc  :  Point); 

{=  Handle  mouse  click  in  the  grow  box  of  a  window  =} 
type 

GrowRec  =  record 

case  integer  of 

0  :  (Result :  Longint); 

1  (Height,  Width  ;  integer); 

end: 

var 

Growinfo  :  GrowRec; 

WindNum  :  integer; 
begin 

if  WPtr  =  FrontWindow  then  {  if  it's  the  active  window  } 
with  Growinfo  do 
begin 

Result  :=  GrowWindow(WRr,  MLoc,  GrowArea);  {  get  amount  of  growth  } 
SizeWindow(WPtr,  Width,  Height,  TRUE);  { resize  window  } 
InvalRect(WRr^.portRect):  {  set  up  for  update  } 

WindNum  >  2  ; 

if  Window(WindNum].P  <>  FrontWindow  then 
OrawGrowlcon(Window[WindNum].P);  {  draw  grow  icons  } 
end 
else 

SelectOurWindow(WPtr): 
end:  {  HandleGrow } 

procedure  HandleDrag(WPtr  :  WindowPtr;  MLoc  :  Point): 

Handle  the  dragging  of  a  window  -} 
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var 

WindNum  :  integer; 
begin 

if  WRr  =  FrontWindow  then 
begin 

DragWindow(WRr,  MLoc,  DragArea);  { in  the  drag  bar } 

WindNum  :=  2  ; 

DrawGrowlcon(Window[WindNum].P);  {  draw  grow  icx>ns  } 
end 
else 

SelectOurWindow(WPtr) ; 
end;  {  HandleDrag } 

procedure  DoMouseDown(TheEvent :  EventRecord); 

{=  identify  where  the  mouse  was  clicked  and  handle  it  =} 
var 

theWindow  : 

WindowPtr; 

MLoc 

Point; 

WLoc 

integer; 

begin 

MLoc  :=  TheEvent.Where;  {  get  mouse  position  } 

WLoc  :=  FindWindow(MLoc,  theWindow) ;{  get  the  window  and  the  location  } 

case  WLoc  of 

InMenuBar  :  HandleMenu(MenuSelect(MLoc));  { in  the  menu  } 
InContent  :  HandleClick(theWindow,  MLoc);  { inside  the  window  } 
InZoomln  :  HandleZoom{theWindow,  MLoc,  WLoc);{  in  the  zoom  box  } 
InZoomOut  :  HandleZoom(the Window,  MLoc,  WLoc);{  in  the  zoom  box  } 
InGoAway  :  HandleGoAway(the Window,  MLoc);  { in  the  go  away  box } 
InGrow  :  HandleGrow(theWindow,  MLoc);  { in  the  grow  box  } 
InDrag 

HandteDrag(theWindow,  MLoc);  { in  the  drag  bar  } 

InSysWindow  :  SystemClick(TheEvent,  theWindow);  { in  a  DA  window  } 
end; 
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end;  { DoK^.ouseDown } 


procedure  DoUpdate(TheEvent :  EventRecord); 

{=  handles  window  update  event  -} 
var 

SavePort, 

theWindow  :  WindowPtr; 
begin 

theWindow  >  WindowPtr{TheEvent.Message);  { find  which  window  } 
if  (theWindow  -  Window[RootLocusWind].P)  or 
(theWindow  =  Window[RegendBox].P) 
then  { only  update  our  windows } 

begin 

MakeWatchCursor; 

GetPort(SavePort):  {  save  current  grafport  } 

SetPort(theWindow):  {  set  as  current  port  } 

BeginU^ate(theWindow);  {  signal  start  of  update } 

{  and  here's  the  update  stuffi } 

if  theWindow  =  Window[RootLocusWind).P  then 
begin 

ClearWindow(RootLocusWind) ; 

DrawGrowlcon(theWindow): 

DrawPic(RootLocusWind); 

end; 

{  now,  back  to  our  program...} 

EndUpdate(theWindow):  {  signal  end  of  update  } 

SetPort(SavePort):  { restore  grafport  } 

MakeArrowCursor;  { restore  cursor  } 

end 

end;  { DoUpdate } 


procedure  DoActivate(TheEvent  :  EventRecord); 


Handles  window  activation  event 


var 

AFlag  ;  boolean; 
theWindow  :  WindowPtr; 
WindNum  ;  integer; 
begin 

with  TheEvent  do 


begin 


theWindow  >  WindowPtr(Message);  { get  the  window  } 

AFlag  :=  Odd(Modifiers);  { get  activate/deactive  } 

if  AFlag  then 

begin  { if  it's  activated...  } 

SetPort(theWindow);  {  fnatf^e  it  the  port  } 

WindNum  >  RootLocusWind  ; 

DrawGrowlcon(Window[WindNum].P);  {  draw  grow  icons  } 
end; 
end 

end:  {  DoActivate  } 

procedure  DoKeypress(theEvent  :  EventRecord); 

{=  handles  keypress  (key[)own,  autoKey)  event  =} 

var 

KeyCh  ;  char; 
begin 

KeyCh  ;=  Chr(theEvent. Message  and  charCodeMask);  { decode  character } 
if  (the  Event. modifiers  and  cmdKey)  <>  0  then 
begin  {  menu  key  command  } 

HandleMenu(MenuKey(KeyCh))  { get  menu  and  item} 

end 
else 

if  TextInputEnabled  then 
begin 

TEKey (KeyCh,  textH); 

TEUpdate(thePort''.portRect,  textH); 
end 
else 

SysBeep(l);  { do 'something*  } 

end;  {  DoKeypress  } 

procedure  initialize; 

{=  Initialize  everything  for  the  program  =} 

var 

Indx 

integer; 

Result  :  real; 

File  Err  :  byte; 

begin  { initialize  all  managers  used  } 

lnitGraf(@thePort):  {  create  a  grafport  for  the  screen  ) 

InitFonts;  { start  up  the  font  manager  } 

InitWindows;  {  start  up  the  window  manager  } 
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InitMenus;  {  start  up  the  menu  manager  } 

TEInit; 

InitGraphic; 

SetUpWindows; 

TextFont(SystemFont):  {  initialize  the  font  } 

{  set  up  menus  } 

MenuUst[AM]  :=  GetMenu(ApplMenu);{  read  menus  in  from  resource  fork  } 
MenuList[FM]  >  GetMenu(FiteMenu); 

MenuList[EM]  :=  GetMenu(EditMenu); 

MenuList[PM]  :=  GetMenu(PlotMenu); 

MenuList[HM]  :=  GetMenu(HelpMenu); 

AddResMenu(MenuList[AM],'DRVR');  {  pull  in  all  desk  accessories  } 
for  Indx  :=  1  to  5  do  {  place  menus  in  menu  bar  } 
lnsertMenu(MenuList[lndx],0); 
for  Indx  :=  1  to  6  do 

SetltemState(EM, Indx, FALSE);  {  deactivate  items  in  Edit  menu  } 
DrawMenuBar;  { draw  updated  menu  bar  to  screen  } 

Finished  :=  False;  {  set  program  terminator  to  false  } 

{  set  drag  region  } 

SetRect(DragArea,  XMinGlb+1,  YMinGlb+38,  XMaxGlb-1,  YMaxGlb-1); 

{  set  grow  region  } 

SetRect(GrowArea.  XMinGlb+1,  YMinGlb+38,  XMaxGlb-1,  YMaxGlb-l); 
InitGuess.Re  :=  1.0; 

InitGuess.lm  0.0; 

Tolerance  :=  1e-6; 

Maxiter  :=  100; 

InitDegree  :=  0; 

AMinGain  ;=  0.1 ; 

AMaxGain  :=  10000; 

BMinGain  >  0.1; 

BMaxGain  :=  10000; 

XMn  >-10; 

XMx  >  5; 

YMn  >  -10; 

YMx  >  10; 

Step  :=  4; 

Points  >  50; 

FillChar(lnitPoly,  SizeOf(lnitPoly),  0); 

FillChar(xAnswer,  SizeOf(xAnswer),  0); 

FillChar(xlnitPoly,  SizeOf(xlnitPoly),  0); 

FillChar(lnfixArray,  SizeOf(lnfixArray),  0); 

MakeWatchCursor; 

MakeInfoScreen; 
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MakeArrowCursor; 

TextInputEnabled  False; 

InitDegreeStatus  False; 

GetCoeffStatus  False; 

end;  {  Initialize  } 

procedure  HandleEvent(TheEvent  :  EventRecord); 

{=  Decodes  an  event  and  handles  it 
begin 

case  The  Event.  What  of 
mouse  Down 

DoMouseDown(TheEvent);  {  mouse  button  pushed  } 
key  Down  :  DoKeyPress(TheEvent);  { key  pressed  down  } 

autoKey  :  DoKeyPress(TheEvent);  { key  held  down  } 
updateEvt  :  DoU|;^ate(TheEvent);  { window  needs  updating  } 
{activateEvt :  DoActivate(TheEvent);}  {  window  made  act/inact  } 
end 

end;  {  HandleEvent } 
procedure  Cleanup; 

{=  Do  any  last  minute  clean  up  work  =} 
begin 

end;  { Cleanup } 


begin  { MacRootLocus } 

Initialize;  { set  everything  up  } 

repeat  { keep  doing  the  following  } 

SystemTask;  {  update  desk  accessories  } 

if  GetNextEvent(everyEvent,theEvent)  then  {  if  there’s  an  event...  } 
HandleEvent(theEvent);  {  ...then  handle  it  } 

until  Finished;  {  until  user  is  done  } 

Cleanup; 

end.  { MacRootLocus } 


64 


unit  GlobalVar  (5000); 

{=  This  unit  declares  the  whole  variables  tobe  used  =} 
{=  in  MacRootLocus  except  a  part  of  local  variables.  -} 

{$U  RootsFinder} 
interface 


uses 

MemTypes,  QuickDraw,  OSIntf,  ToolIntf.PackIntf.PasPrinter, 
SANE.MacPrint,  RootsFinder; 


const 

InfoScrnWind  =  1; 
RootLocusWind  =  2; 
RegendBox  =  3; 
AboutBoxWind  =  4; 
AlertBox  =  5; 
HeIpWind  =  6; 


{  Window  number  of  info  window  } 

{  Window  number  of  root  locus  window } 

{  Window  number  of  informatiom  box  window } 
{  Window  number  of  about  box  window } 

{  Window  number  of  alert  box  window  } 

{  Window  number  of  help  window} 


MenuCnt  =  5;  { total  #  of  menus  } 

AppIMenu  =  1000;  { resource  ID  of  Apple  Menu  } 

FileMenu  =  1001;  { resource  ID  of  File  Menu  } 

EditMenu  >=  1002;  { resource  ID  of  Edit  Menu  } 

PlotMenu  =  1003;  { resource  ID  of  Plot  Menu  } 

HelpMenu  =  1004;  { resource  ID  of  Help  Menu  } 


AM 

FM 

EM 

PM 

HM 


1 ;  { index  into  MenuUst  for  Apple  Menu  } 


2;  { ditto  for  File  Menu  } 

3;  {  ditto  for  Edit  Menu  } 

4;  {  ditto  for  Plot  Menu  } 

5;  {  ditto  for  Help  Menu  } 


MaxPlotGIb  =  1024;{  The  maximum  number  of  points  in  a  plot  array  } 


type 

str255  =  string  [255]  ; 

StringArray  =  array[1..20]  of  str255;  {  Storage  for  plot  data  } 
NodePtr=  ''Node; 

Node  = 


record 

Data :  Extended; 

Next :  NodePtr; 
end; 

PlotArray  «  array[1.. MaxPlotGIb,  1..2I  of  Extended; 
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var 

Answer,  polish  :  str255; 

InfixArray  :  StringArray; 

A.B.ADeltaStep,  AIncrem  .BDeltaStep,  BIncrem:  Extended; 

InitGuess  :  TNComplex;  { Initial  approximation  } 

Tolerance  :  Extended;  { Tolerance  of  approximation } 

AMaxGain,  BMaxGain  ;  Extended;  {  Range  of  unknown  parameter  values } 
AMinGain,  BMinGain  :  Extended; 

Root,  Imag,  Value,  Deriv  ;  TNvector;  {  Resulting  roots  and  other  info  } 

Iter :  TNIntVector;  { Iterations  to  find  each  root } 

Maxiter :  integer;  {  Maximum  number  of  iterations  } 

InitDegree,  Degree  :  integer;  { Initial  and  final  degree  } 
xInitPoly  :  TNvector;  { of  polynomial } 

InitPoly,  Poly  :  TNCompVector;  { Initial  and  final  coefficients  } 

{ of  the  polynomial } 
xAnswer ;  TNCompVector; 
yAnswer :  TNCompVector; 

NumRoots,  Arrayindex,  pointno  ;  integer;  {  Number  of  roots  } 

Error :  byte;  {  Error  flag  } 

StepA,  StepB,  linear :  Boolean; 

Step, Points  :  Longint; 


OutFile  :  text; 
OutName  ;  string; 
lOerr  :  boolean; 


{  The  output  file  used  by  an  application} 
{ The  out  file  name} 

{  Flags  I/O  errors } 


Finished  ;  boolean;  {  used  to  terminate  the  program  } 

theEvent  :  EventRecord;  {  event  passed  from  operating  system  } 

MenuList  :  array[1..MenuCnt]  of  MenuHandle;  {  holds  menu  info  } 
DragArea  :  Rect  ;  {  Area  in  which  window  can  be  dragged } 

GrowArea  :  Rect;  {  Area  in  which  a  window's  size  can  change  } 


DataPicture  :  PicHandle;  {  Holds  a  picture  of  the  plotted  data  } 

theDialog  :  DialogPtr; 
itemHit  ;  Integer; 
theType  :  Integer; 
r  :  Rect; 

done,  AutoScale  :  Boolean; 
n  :  Integer; 
h,h3,h4,h5,h6  :  Handle; 

s  ;  Str255; 

XMn  :  Extended; 

XMx  :  Extended; 

YMn  ;  Extended; 

YMx  :  Extended; 


Graph  Array  ;  '^PIotArray; 
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Ds  :  Deci  r : 

AR Justification,  AMarkStatus.BRJustification,  BMarkStatus 
Da,  Db,  Dan,  Dax,  Dbn,  Dbx,  Dss  :  DecStr; 
txRect :  Rect; 
textH  :  TEHandle; 

TextInputEnabied :  Booiean; 

InitDegreeStatus  ;  Boolean; 

GetCoeffStatus  :  Boolean; 

implementation 

begin 

end.{GlobatVar} 


Boolean; 
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unit  MakeRoot(3000); 


{=  This  unit  provides  several  procedures  and  functions  to  == } 

{=  make  the  roots  and  the  array  vectors  for  ptot.  = } 

{$!-}  {  Disable  I/O  error  trapping  } 

{$8+}  {  Enable  segmentation  of  code } 

{$U  SpecVar} 

{$U  Message} 

{$U  RootsFinder} 

{$U  TurboGraphj 

interface 

uses 

MemTypes,  QuickDraw,  OSIntf,  Toolintf,  Packlntf.PasPrinter, 

SANE,  MacPrint, RootsFinder,  SpecVar, 

{$S  SecondSegment) 

TurboQraph, 

Message; 

var 

Error  :  byte; 

procedure  GetRootf ; 
procedure  GetRoot2; 

implementation 


procedure  lnfixToPolish(  Answer;str255;var  RPN:str255): 

{=  This  procedure  converts  an  algebraic  expression(infix)  = } 
{=  to  reverse  Polish  notation.  = } 


var 

Stack  ;  Array [1  ..50]  of  char; 
Top,l,N,p  :  Integer; 

Ch  ;  char; 

Ch1  :  string; 

Firstchar.PrevDigit  ;  Boolean; 

Function  Prlority(Ch:char):  Integer; 
Begin 
Case  Ch  of 
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:  Priority;*  4; 

:  Priority:*  3; 
r  :  Priority:*  3; 

:  Priority:*  2; 
'+'  :  Priority:*  2; 

'a'..'b'  :  Priority:*  1; 

:  Priority:-  1 : 
'0'..'9'  :  Priority:*  1; 


1  1 

:  Priority:*  1 ; 

'(■ 

:  Priority:*  0; 

■)■ 

:  Priority:*  0 

1  1 

end; 

end; 

begin 

RPN  :=  "; 

;  Priority;*  0; 

Top  :=0: 

FirstChar  :=True: 

PrevDigit  ;*  False; 

for  i:*  1  to  Length(Answer)  do 
begin 

chi  ;»  copy(Answer,i,1): 
ch:=  ch1[11: 

P:=  Priority(ch); 

if  FirstChar  and  (Ch  =  then  P:*  1 ; 

if  P  =  1  then 
begin 

if  PrevDigit  then 
RPN  ;*  concat(RPN,ch) 
else 

RPN  :*  concat(RPN,'  ’,ch): 
firstchar  ;*  False;  PrevDigit  :*  True; 
end; 

if  P>1  then 
begin 

while  (Top>0)  and  (Priority(stack[top])  >=  P)  do 
begin 

RPN  ;*  concat(RPN,'  '.Stackfrop]); 

Top  :»  Top  - 1 ; 
end; 

Top  :-  Top  +  1; 
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StackfTop]  :=  ch; 

Firstchar  :=  True;  PrevDigit  ;=  False; 
end; 
if  ch  -  ’(' 
begin 

Top  >  Top  +1 ; 

Stack[Top]  := '('; 

FirstChar  :=  True;  PrevDigit  :=False; 
end; 

if  ch  =  ')'  then 
begin 

whiie  Stack[Top]<> '('  do 
begin 

RPN  :=  ooncat(RPN,'  ',Stack[Top]); 

Top  ;=  Top  -1 ; 
end; 

Top  ;=  Top  -1 ; 

Firstchar  ;=  False;  PrevDigit  ;=  False; 
end; 
end; 

while  Top  >  0  do 
begin 

RPN  :=  concat(RPN,'  ',Stack[Topl); 

Top  :=  Top  -1 ; 
end; 
end; 

procedure  ComputePolish(polish:str255;a,b:Extended;  var  EvalArray;Extended); 


{=  This  procedure  uses  the  string  generated  in  Procedure  = } 

{=  InfixToPolish  to  evaluate  the  numeric  expression  =} 

var 

j.N  ;  integer; 
ch.ch2  ;  char; 
chi  :  string; 

temp  :  string(255]; 

Valuel  ,Value2, Values  :  Extended; 

StackPtr  :  NodeRr; 

StackEmpty :  Boolean; 

procedure  CreateStack; 

{«  Initialize  stack  -} 

begin 


70 


New(StackPtr): 
with  StackPtr''  do 
begin 

Next  :=  nil; 

Data  >  0.0; 
end; 

StackEmpty  :=  True; 
end; 

procedure  Pop(var  Val  :  Extended);{Push  a  number  onto  numeric  stack} 
{=  pop  a  number  off  numeric  stack  =} 
var 

NPtr  ;NodePtr; 
begin 

if  not  StackEmpty  then 
begin 

NPtr  :=  StackPtr'^.Next; 

StackPtr'^.Next  :=  NPtr''. Next; 

Val  :=  NPtr^'.Data; 

Dispose(NPtr); 

StackEmpty  :=  (StackPtr^'.Next  =  nil); 
end; 
end; 

procedure  Push(Val  :  Extended); 

=  pop  a  number  off  numeric  stack  =} 

=  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =} 
var 

NPtr ;  NodePtr; 
begin 

StackEmpty  ;=  False; 

New(NPtr); 

NRr''.Data  :=  Val; 

NRr^.Next  :=  StackPtr^'.Next; 

StackPtr^'.Next  :«  NPtr; 

end; 

procedure  DeleteStack;{ Delete  stack} 

{ssssssassss} 

{=  Delete  stack  s} 

{saasaasasas} 

var 

Temp :  Extended; 
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begin 

while  not  StackEmpty  do 
Pop(Temp); 
Dispose(StackPtr); 
end; 


function  Expon(y,x:Extended):Extended; 
{-  Computes  Y  raised  to  X  power 
begin 

Expon  exp(  x  *  {ln(y))): 
end; 


begin 

CreateStack;{lnitialize} 
temp  ;=' '; 

for  i:=  1  to  Length(polish)  do{do  one  char  at  a  time} 
begin 

chi  :=  copy(polish,i,1); 
ch:»  ch1[1l;{Get  a  char} 

case  ch  of  {and  evaluate  it} 

:  temp  :*concat(temp,ch);{Real  constant} 

:  temp  :=  concat(temp,ch); 

:  begin 

ch1  :=  copy(polish,i+1,1); 
ch2:=  ch1[1];{Get  a  char} 

if  (ch2  <>  ' ')  and  (i  <  length(polish))  then{Unary  minus} 
temp  :=  concat(temp,ch) 
else 

begin  {minus  operator} 

POP(Value1  ):POP(Value2); 

Values  Va'ue2  •  Valuel ; 

PUSH(Value3); 

end; 

end; 


'a' 

:  PUSH(a); 

•A’ 

:  PUSH(A); 

'b' 

:  PUSH(b); 

•B' 

:  PUSH(B); 

’+•  i)egin 

POP(Valuel);  POP(Value2); 
Values  Values  +  Valuel ; 

PUSH(Value3); 

end; 

’**  iboQin 

POP(Valuel):  POP(Value2); 
Values  Values  *  Valuel ; 
PUSH(  Values): 

end; 


T  :begin 

POP(Valuel):  POP(Value2): 
Values  :=  Values  /  Valuel ; 
PUSH(Value3): 
end; 


:begin 

POP(Valuel):  POP( Values); 
Values  :»Expon(  Values  ,  Valuel); 
PUSH(  Values): 
end; 

' '  :begin 

if  temp  <> '  'then 
begin 

Valuel  :=  Str2Num(temp); 
PUSH(Valuel): 
temp  :=  "; 
end; 
end; 


end; 

end; 

if  temp  <>  ' '  then 
begin 

Values  Str2Num(temp); 
PUSH{Value3); 
end; 

POP(Valuel): 
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EvalArray  :>  Valuel ; 

DeleteStack; 

end; 

function  Ten2  (power:  Extended)  :  Extended; 
begin 

Ten2  Exp(power  *  Ln(10)); 
end; 


function  NextGainI  :  Extended; 

{=  This  function  makes  two  kind  of  intervals  for  one  = } 

{=  parameter  root  locus  method.  = } 

var 

GainOut,  logmin,  logmax  :  Extended; 


begin 

if  linear  then 
begin 

AIncrem  :=  abs((AMaxGain  •  AMinGain)/(points  -  1)); 
GainOut  :=  AMinGain  +  AIncrem  *  pointno; 
end 
else 
begin 

if  (AMinGain  <  1E-5)  and  (AMaxGain  >  1E-1)  then 
logmin  :=  -5 
else 

logmin  :=  Ln(AMinGain)/Ln(10); 
logmax  :=  ln(AMaxGain)/Ln(10); 

AIncrem  :=  (logmax  -  logmin)/(Points  -1); 

GainOut  :=  Ten2(logmin  +  pointno  *  AIncrem); 
end; 

NextGainI  GainOut; 
end; 


function  NextGain2  :  Extended; 

{=  This  function  makes  two  kind  of  intervals  for  each  » } 

{=  parameter  of  two  parameter  root  locus  method.  =} 

var 

GainOut,  logmin,  logmax  :  Extended; 


begin 

if  StepA  then 
begin 
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if  linear  then 
begin 

BIncrem  :=  abs((BMaxGain  -  BMinGain)/(points  -  1)); 
GainOut  ;=  BMinGain  +  BIncrem  *  pointno; 
end 
else 
begin 

if  (BMinGain  <  IE-5)  and  (BMaxGain  >  1E-1)  then 
logmin  :=  -5 
else 

logmin  >  Ln(BMinGain)/Ln(10); 
logmax  ln(BMaxGain)/Ln(10); 

BIncrem  >  (logmax  -  logmin)/(Points  -1); 

GainOut  >  Ten2(logmin  +  pointno  *  BIncrem); 
end; 

NextGain2  :=  GainOut; 

end 

else 

begin 

if  linear  then 
begin 

AIncrem  :=  abs((AMaxGain  -  AMinGain)/(points  -  1)); 
GainOut  :=  AMinGain  +  AIncrem  *  pointno; 
end 
else 
begin 

if  (AMinGain  <  IE-5)  and  (AMaxGain  >  IE-1)  then 
logmin  :=  -5 
else 

logmin  ;=  Ln(AMinGain)/Ln(10); 
logmax  :=  ln(AMaxGain)/Ln(10); 

AIncrem  :=  (logmax  -  logmin)/(Points  -1); 

GainOut  :=  Ten2(logmin  +  pointno  *  AIncrem); 
end; 

NextGain2  >  GainOut; 
end; 
end; 

procedure  Results(A,B  :  Extended; 

xAnswer  :  TNCompVector; 

Error  :  byte); 


{=  This  procedure  outputs  the  results  to  the  device  OutFile 
{=  and  make  the  array  for  plot  data.  =} 
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var 

Term  ;  integer; 
begin 

if  Degree  >  0  then 
begin 

Writeln(OutFile,The  deflated  polynomial;'); 
for  Term  :=  Degree  downto  0  do 
Writeln(OutFile,  'Poly[',Term,'J:', 

Polyn‘0miJ.Re:23.  '  +',  PolyfTerm].lm:23,'i'); 

end; 

if  Error  <=  1  then 
begin 

Writeln(OutFile); 

Writeln(OutFile,  ’  A  =  ’  ,  A  :10:5.'  B  =  B  :10:5  ); 
for  Term  :=  1  to  NumRoots  do 
begin 

Writeln(OutFile,  xAnswer[Term].Re:23,  ’  +',  xAnswer[Term].lm:23,  'i'); 
GraphArray^[Arraylndex,1]  :=  xAnswer(Term].Re; 

GraphArray''[Arraylndex,2]  :=  xAnswer[Term].lm; 

Arrayindex  :=  Arrayindex  +  1 ; 
end; 
end; 

end:  {  procedure  Results  } 
procedure  PlotRootLocusI ; 

{=  This  procedure  draw  a  plot  of  one  parameter  root  locus  method  =) 

begin 

SelectWind(2,  TRUE); 

DefineHeader(2,  'One  Parameter  RootLocus'); 

DrawGrowlcon{Window[2].P): 

OpenPic(2,  TRUE);  {  Open  up  a  picture  for  window  #1  } 

{  Draw  the  axis  and  data  points  } 
if  AutoScale  then 

FindWorid(2,  GraphArray'',  Arrayindex) 
else 

FindWorldf  (2,XMn,YMn,XMx,YMx); 

DrawAxis(  Da  ,  ",  True): 

DrawPolygon(GraphArray'',  1,  -Arrayindex,  -3,  1,  0,  TRUE); 

ClosePicture;  {  Close  the  picture  for  window  #2} 
ValidRect{Window[2].P''.portRect): 
end:  {  PlotPowerExpLog  } 
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procedure  PlotRootLocus2; 


This  procedure  draw  a  plot  of  two  parameter  root  locus  method 

begin 

SelectWind{2,  TRUE): 

DefineHeader(2,  'Two  Parameter  RootLocus'); 
DrawGrowlcon(Window[2].P); 

OpenPic(2,  TRUE);  {  Open  up  a  picture  for  window  #2  } 

{  Draw  the  axis  and  data  points } 

FindWorldl  (2,XMn,YMn,XMx.YMx); 

DrawAxis(Da,  Db,  True); 

DrawPolygon(GraphArray'',  1,  -Arrayindex,  -3,  1,  0,  TRUE): 
ClosePicture;  {  Close  the  picture  for  window  #2  } 
ValidRect(Window(2].P''.portRect); 
end: { } 

procedure  GetRoott ; 

This  procedure  gets  the  root  of  the  characteristic  equation  for 
one  parameter  root  locus  method. 

var 

term  :  integer; 
f :  DecForm; 
begin 

Num2Str(f,AMinGain,Dan); 

Num2Str(f,AMaxGain,Dax) ; 

Da  ;=  'A  ; '+  Dan  +  '  To  ’+Dax; 

Ds  :=  ' 

Arrayindex  ;=  1; 

New  (GraphArray); 
pointno  :=  0; 

while  pointno  <=  (  Points  -  1)  do 
begin 

A  :=  NextGainI ; 

B  0; 

for  Term  :=  InitDegree  downto  0  do 
begin 

infixToPolish(lnfixArray[Term],polish); 

ComputePolish(polish,a,b,xlnitPoly[Term]); 

InitPolyfTermJ.Re  :=  xlnitPoly[TermJ: 
lnitPolyrrerm].lm  :=  0.0; 
end; 

Degree  :*  InitDegree; 

Poly  >  InitPoly; 


Laguerre(Degree,  Poly,  InitGuess,  Tolerance,  Maxiter, 

NumRoots,  xAnswer,  yAnswer,  Iter,  Error); 

Results(A,  B,  xAnswer,  Error); 
pointno  :=  polntno  +  1 ; 
end; 

Arrayindex  :*  Arrayindex  -  1; 

PlotRootLocusI  ; 
end; 

procedure  GetRoot2{(var  StepA  :  boolean;  A,B  :  Extended)}; 

{=  This  procedure  gets  the  root  of  the  characteristic  equation  for  =} 
{=  two  parameter  root  locus  method.  =} 

var 

i,  j,  k,  term  :  integer; 
f :  DecForm; 
begin 

ADeltaStep  ;=  abs((AMaxGain  -  AMinGain)/(Step)); 

BDeltaStep  :=  abs((BMaxGain  -  BMinGain)/(Step)); 
Num2Str(f,AMinGain,Dan); 

Num2Str(f,AMaxGain,Dax); 

Da  :=  'A  : '+  Dan  +  '  To  '+Dax; 

Num2Str(f,BMinGain,Dbn); 

Num2Str(f,BMaxGain,Dbx) ; 

Db  ;=  'B  :  ’+  Dbn  +  *  To  ’+Dbx; 

New  (Graph Array): 

StepA  :=  True; 
for  i  :=  1  to  2  do 
begin 

for  j  :=  1  to  Step  do 
begin 

Arrayindex  ;=  1; 

if  StepA  then  A  :=  AMinGain  +  ADeltaStep*{j-1) 
else  B  :=  BMinGain  +  BDeltaStep*(j-1); 
if  StepA  then 
begin 

Num2Str(f,A,Dss): 

Ds  :■  'A='+  Dss; 
end 
else 
begin 

Num2Str(f,B,Dss): 

Ds  >  'B>'-f  Dss; 
end; 

pointno  >  0; 

while  pointno  <»  (  Points  -  1 )  do 
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begi 

if  StepA  then  B  >  NextGain2 
else  A NextGain2; 
for  Term  :=  InitDegree  downto  0  do 
begin 

infixToPolish(lnfixArray[Term], polish); 
ComputePolish(polish,a,b,xlnltPoly[Term]); 
lnitPoly[Term].Re  :=  xlnitPoly[Term]; 
lnitPoly[Term].lm  :=  0.0; 
end; 

Degree  ;=  InitDegree; 

Poly  :«*  InitPoly; 

Laguerre(Degree,  Poly,  InitGuess,  Tolerance,  Maxiter, 
NumRoots,  xAnswer,  yAnswer,  Iter,  Error); 
Results(A,  B,  xAnswer,  Error); 
pointno  :=  pointno  +  1 ; 
end; 

Arrayindex  :=  Arrayindex  -  1 ; 

PlotRootLocus2; 

end; 

StepA  :=  False; 
end; 
end; 
begin 

end.  {  Unit  MakeRoot } 
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unit  RootsFinder(2000); 


Turtx)  Pascal  Numerical  Methods  Toolt)ox 
Copyright  (C)  1987  Borland  International 

This  unit  provides  procedures  for  finding  the  roots 
of  a  single  equation  in  one  real  variable. 


{$R+}  {  Enable  range  checking  } 

interface 

uses 

MemTypes; 

u)nst 

TNNearlyZero  =  IE-015:  {  Close  to  zero  } 

TNArraySize  =10;  {  Maximum  size  of  vectors  } 

type 

TNvector  =  array[O..TNArraySizel  of  Extended; 

TNIntVector  =  array[O..TNArraySize]  of  integer; 

TNcomplex  =  record 

Re,  Im  :  Extended: 
end; 

TNComp Vector  =  array(0.. TNArraySize]  of  TNcomplex; 


procedure  Laguerre(var  Degree  ;  integer; 
var  Poly  :  TNCompVector; 

InitGuess ;  TNcomplex; 

Tol  :  Extended; 

Maxiter  :  integer; 
var  NumRoots  :  integer; 
var  Roots  :  TNCompVector; 
var  yRoots  :  TNCompVector; 
var  Iter  :  TNIntVector; 
var  Error  :  byte): 


{- 

4 

- 

{= 

Input:  Degree,  Poly,  InitGuess,  Tol,  Maxiter 

«} 
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Output:  Degree,  Poly,  NumRoots,  Roots,  yRoots,  Iter,  Error 


1= 

(= 

[= 

[= 

[= 

[= 

(= 

[= 

[= 

[= 

[= 

[= 

(= 

(= 

(= 

[= 


[= 

[= 

[= 

(= 

{= 

(= 

{= 

(= 

(= 

(= 

{= 

(= 

(= 

(= 

(= 

[= 

(= 


Purpose:  This  unit  provides  a  procedure  for  finding  all  the 
roots  (real  and  complex)  to  a  polynomial. 

Laguerre's  method  with  deflation  is  used. 

The  user  must  input  the  coefficients  of  the  quadratic 
and  the  tolerance  in  the  answers  generated. 

Pre-defined  Types:  TNcomplex  -  record 
Re,  Im  :  Extended; 
end; 

TNIntVector  =  array[O..TNArraySize]  of  integer; 
TNCompVector  =  array[O..TNArraySize]  of  TNcomplex; 

Variables:  Degree  :  integer;  degree  of  deflated 

polynomial 

Poly  :  TNCompVector;  coefficients  of  deflated 
polynomial  where  PolyIn]  is 
the  coefficient  of  X'^n 

InitGuess  :  TNcomplex;  initial  guess  to  a  root 
(may  be  very  crude) 

Tol  :  Extended;  tolerance  in  the  answer 
Maxiter  :  integer;  number  of  iterations 
NumRoots  :  integer;  number  of  roots  calculated 
Roots  :  TNCompVector;  the  roots  calculated 
yRoots  :  TNCompVector;  the  value  of  the  function 
at  the  calculated  roots 

Iter  :  TNIntVector;  number  iteration  it  took  to 
find  each  root 

Error  :  byte;  flags  an  error 

Errors:  0:  No  error 
2:  Degree  <=  0 
3:  Tol  <=  0 
4:  Maxiter  <  0 


=} 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 

=  } 

=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=  } 
=} 
=  } 


implementation 


The  following  inline  procedure  and  functbn  are  used  to  call  the  user  » } 
defined  procedures  and  functions  pointed  to  by  the  ProcAddr  parameter.  >= } 


function  UserFunction(X  :  Extended;  ProcAddr :  ProcPtr)  :  Extended; 
inline 

$205F,  {  MOVE.L  (A7)+.  AO  } 

$4E90;  {JSR  (AO)  } 

procedure  UserProcedure(X  ;  TNcomplex;  var  Y  :  TNcomplex;  ProcAddr  :  ProcRr); 
inline 

$205F,  {  MOVE.L  (A7)+,  AO  } 

$4E90;  {JSR  (AO)  } 


procedure  Laguerre{(var  Degree  ;  integer; 
var  Poly  :  TNCompVector; 
InitGuess :  TNcomplex; 

Tol  :  Extended; 

Maxiter  :  integer; 
var  NumRoots  :  integer; 
var  Roots  :  TNCompVector; 
var  yRoots  :  TNCompVector; 
var  Iter  :  TNIntVector; 
var  Error  ;  byte)}; 


type 

TNquadratic  =  record 

A,  B,  C  :  Extended; 
end; 


var 

Additer  :  integer; 

InitDegree  ;  integer; 

InitPoly  :  TNCompVector; 

GuessRoot  :  TNcomplex; 

{ .  Here  are  a  few  complex  operations  . } 

procedure  Conjugate(var  Cl,  C2  :  TNcomplex); 
begin 

C2.Re  C1  .Re; 

C2.lm  :=  -1.0  *  Cl.lm; 
end;  {  procedure  Conjugate  } 

function  Modulus(var  Cl  :  TNcomplex)  :  Extended; 
begin 

Modulus  :«  Sqrt(Sqr(C1  .Re)  +  Sqr(Cl.lm)); 
end;  { function  Modulus  } 
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procedure  Add(var  C1,  C2,  C3  :  TNcomplex); 
begin 

C3.Re  ’=  C1  .Re  +  CS.Rei 
C3.lm  :=  Cl.lm  +  C2.lm: 
end;  {  procedure  Add } 

procedure  Sub(var  Cl,  C2,  C3  :  TNcomplex); 
begin 

C3.Re  :=  C1.Re-C2.Re; 

C3.lm  Cl.lm  -  C2.lm; 
end;  {  procedure  Sub } 

procedure  Mult(var  C1 ,  C2,  C3  :  TNcomplex); 
begin 

C3.Re  :=  Cl  .Re  *  C2.Re  -  Cl.lm  *  C2.lm; 

C3.lm  :=  Cl.lm  *  C2.Re  +  Cl.Re  *  C2.lm; 
end;  {  procedure  Mult  } 

procedure  Divide(var  C1,  C2,  C3  ;  TNcomplex); 
var 

Dum1,  Dum2  :  TNcomplex; 

E  :  Extended; 
begin 

Conjugate(C2,  Dum1); 

Mult(Cl,  Dum1,  Dum2); 

E  :=  Sqr(Modulus(C2)); 

C3.Re  :=  Dum2.Re  /  E; 

C3.lm  :=  Dum2.lm  /  E; 
end;  {  procedure  Divide  } 

procedure  Square Root(var  C1,  C2  :  TNcomplex); 
var 

R,  Theta :  Extended; 
begin 

R  >  Sqrt(Sqr(Cl.Re)  +  Sqr(Cl.lm)); 
if  ABS(CI.Re)  <  TNNearlyZero  then 
begin 

if  Cl.lm  <  0  then 
Theta  Pi !  2 
else 

Theta  (-1.0  *  Pi)  /  2; 
end 
else 

if  Cl.Re  <  0  then 

Theta  ArcTan(C1.lm  /  Cl.Re)  +  Pi 
else 

Theta  ArcTan{C1.lm  /  Cl.Re); 
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C2.Re  >  Sqrt(R)  *  Cos(Theta  /  2); 

C2.lm  ;=  Sqrt(R)  *  Sin(Theta  /  2); 
end;  {  procedure  SquareRoot } 

procedure  lnitAndTest(var  Degree  :  integer; 

var  Poly  :  TNCompVector; 

Tot  :  Extended; 

Maxiter  :  integer; 

InitGuess  iTNcomplex; 
var  NumRoots  :  integer; 
var  Roots  :  TNCompVector; 
var  yRoots  :  TNCompVector; 
var  Iter  :  TNlntVector; 
varGuessRoot  :TNcomplex; 
var  InitDegree  :  integer; 
var  InitPoly  ;  TNCompVector; 
var  Error  :  byte); 


[=  Input:  Degree,  Poly,  Tol,  Maxiter,  InitGuess 
(=  Output:  InitDegree,  InitPoly.  Degree.  Poly,  NumRoots, 
{=  Roots,  yRoots,  Iter,  GuessRoot,  Error 

(= 

{=  This  procedure  set*:  the  initial  value  of  the  above 
(=  variables.  This  procedure  also  tests  the  tolerance 
(=  (Tol),  maximum  number  of  iterations  (Maxiter),  and 
(=  code.  Finally,  it  examines  the  coefficients  of  Poly. 

[=  If  the  constant  term  is  zero,  then  zero  is  one  of  the 
[=  roots  and  the  polynomial  is  deflated  accordingly.  Also 
(=  if  the  leading  coefficient  is  zero,  then  Degree  is 
[=  reduced  until  the  leading  coefficient  is  non-zero. 


var 

Term  :  integer; 
begin 

Error  :-  0; 
if  Degree  <»  0  then 

Error :-  2;  {  degree  is  less  than  2  } 

if  Tol  <-  0  then 
Error  :-  3; 
if  Maxiter  <  0  then 
Error  :-  4; 

if  Error  ■  0  then 
begin 
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NumRc  'ts  :=  0; 

GuessRoot  >  InitGuess; 

InitDegree  :«  Degree; 

InitPoly  :*  Poly; 

{  Reduce  degree  until  leading  coefficient  <>  zero } 
while  (Degree  >  0)  and  (Modulus(Poly[Degreel)  <  TNNearlyZero)  do 
Degree  Pred(Degree); 

{  Deflate  polynomial  until  the  constant  term  <>  zero  } 
while  (Modulus(Poly[0])  =  0)  and  (Degree  >  0)  do 
begin 

{  Zero  is  a  root } 

NumRoots  Succ(NumRoots); 

Roots[NumRoots].Re  :=  0; 

Roots[NumRoots].lm  >  0; 
yRoots[NumRoots].Re  ;=  0; 
yRoots[NumRoots].lm  :=  0; 

IterfNumRoots]  :=  0; 

Degree  :=  Pred(Degree); 
for  Term  ;=  0  to  Degree  do 

Poly[Term]  ;=  Poly[Term  +  1]; 

end; 

end; 

end;  {  procedure  InitAndTest } 

procedure  FindOneRoot(Degree  :  integer; 

Poly  ;  TNCompVector; 

GuessRoot ;  TNcomplex; 

Tol  :  Extended; 

Maxiter  :  integer; 
var  Root  ;  TNcomp'ex; 
var  y  Value  ;  TNcomplex; 
var  Iter  :  integer; 
var  Error  :  byte); 


{=  Input:  Degree,  Poly,  GuessRoot,  Tol,  Maxiter  = } 

{=  Output;  Root,  yValue,  Iter,  Error  = } 

{=  =} 
This  procedure  approximates  a  single  root  of  the  polynomial  = } 

{>  Poly.  The  root  must  be  approximated  within  Maxiter  - } 

{«  iterations  to  a  tolerance  of  Tol.  The  root,  value  of  the  > } 

{«  polynomial  at  the  root  (yValue),  and  the  number  of  iterations  - } 

{>  (Iter)  are  returned.  If  no  root  is  found,  the  appropriate  error  ic } 

{-  code  (Error)  is  returned.  - } 


var 
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Found :  boolean; 

Dif  :  TNcomplex; 

yPrime,  yDoublePrime  :  TNcomplex; 

procedure  EvaluatePoly(Degree  :  integer; 
Poly  ;  TNCompVector; 

X  :  TNcomplex; 
var  y  Value  :  TNcomplex; 

var  yPrime  :  TNcomplex; 
var  yDoublePrime  :  TNcomplex); 


{=  Input;  Degree,  Poly,  X  = } 

{=  Output:  yValue,  yPrime,  yDoublePrime  = } 

{=  This  procedure  applies  the  technique  of  synthetic  division  to  = } 

{=  determine  value  (yValue),  first  derivative  (yPrime)  and  second  = } 

{=  derivative  (yDoublePrime)  of  the  polynomial.  Poly,  at  X.  = } 

{=  The  0th  element  of  the  first  synthetic  division  is  the  = } 

{=  value  of  Poly  at  X,  the  1st  element  of  the  second  synthetic  = } 

{=  division  is  the  first  derivative  of  Poly  at  X,  and  twice  the  = } 

{=  2nd  element  of  the  third  synthetic  division  is  the  second  = } 

{=  derivative  of  Poly  at  X.  = } 


var 

Loop  :  integer; 

Dummy,  yDPdummy  :  TNcomplex; 

Deriv,  Deriv2  :  TNCompVector; 

begin 

DerivfDegree]  Poly[Degree]; 
for  Loop  :=  Degree  - 1  downto  0  do 
begin 

Mult(Deriv[Loop  +  1],  X,  Dummy); 
Add(Dummy,  Poly[Loop],  Deriv[Loop]); 
end; 

yValue  :=  DerivlO];  {  Value  of  Poly  at  X  } 

Deriv2[Degree]  Deriv[Degree]; 
for  Loop  Degree  - 1  downto  1  do 
begin 

Mult(Deriv2[Loop  +  1],  X,  Dummy); 
Add(Dummy,  Deriv[Loop],  Oeriv2[Loop]); 
end; 

yPrime  :«  Deriv2[11;  {  1st  deriv.  of  Poly  at  X  } 
yDPdummy  Oeriv2[Degree]; 
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for  Loop  >  Degree  - 1  downto  2  do 
begin 

Mull(yDPdummy,  X,  Dummy); 

Add(Dummy,  Deriv2[Loop],  yDPdummy); 
end: 

yDoublePrime.Re  :=  2  *  yDPdummy.Re;  {  2nd  derivative  of  Poly  at  X  } 
yDoublePrime.lm  :=  2  *  yDPdummy.lm; 
end;  {  procedure  EvaluatePoly  } 

procedure  ConstructDifference(Degree  :  integer; 
yValue  :  TNcomplex; 

yPrime  :  TNcomplex; 

yDoublePrims  :  TNcomplex; 
var  Dif  ;  TNcomplex); 


{=  Input:  Degree,  yValue,  yPrime,  yDoublePrime  = } 

{=  Output;  Dif  = } 

{=  =} 

{=  This  procedure  computes  the  difference  between  approximations;  = } 

{=  given  information  about  the  function  and  its  first  two  = } 

{=  derivatives.  = } 


var 

yPrimeSQR,  yTimesyDPrime,  Sum,  SRoot, 

Numen ,  Numer2,  Numer,  Denom  :  TNcomplex: 

begin 

Mult(yPrime,  yPrime,  yPrimeSQR); 
yPrimeSQR.Re  ;=  Sqr(Degree  •  1)  *  yPrimeSQR. Re; 
yPrimeSQR.Im  ;=  Sqr(Degree  -  1)  *  yPrimeSQR.Im; 

Mult(yValue,  yDoublePrime,  yTimesyDPrime); 
yTimesyDPrime.Re  ;=  (Degree  -  1)  *  Degree  *  yTimesyDPrime. Re; 
yTimesyDPrime. Im  ;*  (Degree  -  1)  *  Degree  *  yTimesyDPrime. Im; 
Sub(yPrimeSQR,  yTimesyDPrime,  Sum); 

SquareRoot(Sum,  SRoot); 

Add(yPrime,  SRoot,  Numert): 

SubiyPrIme,  SRoot,  Numer2); 
if  Modulus(Numer1 )  >  Modulus(Numer2)  then 
Numer Numert 
else 

Numer :»  Numer2; 

Denom.Re  :»  Degree  *  yValue.Re; 

Denom.lm  ;«  Degree  *  yValue.lm; 
if  Modulus(Numer)  <  TNNearlyZero  then 
begin 
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Dif.Re  :=  0; 

Dif.lm  :=  0; 
end 
else 

Divide(Denom,  Numer,  DiO;  {  The  difference  is  the  } 

{ inverse  of  the  fraction } 
end;  {  procedure  ConstructDifference  } 

function  TestForRoot(X,  Dif,  Y,  Tol  :  Extended) :  boolean; 

{=  These  are  the  stopping  criteria.  Four  different  ones  are 
{=  provided.  If  you  wish  to  change  the  active  criteria,  simply 
{=  comment  off  the  current  criteria  (including  the  appropriate  OR) 
{=  and  remove  the  comment  brackets  from  the  criteria  (including 
{=  the  appropriate  OR)  you  wish  to  be  active. 


begin 

TestForRoot  ;= 

(ABS(Y)  <=  TNNearlyZero) 

or 

(ABS(Dif)  <  ABS(X  *  Tol)) 


(*  or  *) 

(*  *) 

(*  (ABS(Dif)  <  Tol)  *) 

(*  *) 

(*  or  *) 

(*  *) 
(*  (ABS(Y)  <=  Tol)  *) 


{ . } 

{-  Y=0  -} 

{-  -} 

{-  -} 

{-  -} 

{-  Relative  change  in  X  -} 

{-  -} 

{-  -} 

{-  -} 

{-  -} 

{-  Absolute  change  in  X  -} 

{-  -} 

{-  -} 

{-  -} 

{-  Absolute  change  in  Y  -} 

{ . } 


(n.  The  first  criteria  simply  checks  to  see  if  the  value  of  the 
{=  function  is  zero.  You  should  probably  always  keep  this  criteria 
{=  active. 

{= 

{=  The  second  criteria  checks  the  relative  error  in  X.  This  criteria 
{«  evaluates  the  fractional  change  in  X  between  iterations.  Note 
{»  that  X  has  been  multiplied  throught  the  inequality  to  avoid  divide 
by  zero  errors. 

{= 
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{=  The  third  criteria  checks  the  absolute  difference  in  X  between  = } 

{=  iterations.  = } 

{=  =} 

{=  The  fourth  criteria  checks  the  absolute  difference  between  = } 

{=  the  value  of  the  function  and  zero.  > } 


end;  {  procedure  TestForRoot } 

begin  {  procedure  FindOneRoot } 

Root  :=  GuessRoot; 

Found  :=  false; 

Iter  ;=  0; 

EvaluatePoly(Degree,  Poly,  Root,  y Value,  yPrime,  yDoublePrime); 

while  (Iter  <  Maxiter)  and  not(Found)  do 

begin 

Iter  :=  Succ(lter); 

ConstructDifference(Degree,  y Value,  yPrime,  yDoublePrime,  Dif); 

Sub(Root,  Dif,  Root); 

EvaluatePoly(Degree,  Poly,  Root,  yValue,  yPrime,  yDoublePrime); 

Found  :=  TestForRoot(Modulus(Root),  Modulus{Dif),  Modulus(yValue),  Tol): 
end; 

if  not(Found)  then  Error  :=  1 ;  {  Iterations  execeeded  Maxiter } 
end;  {  procedure  FindOneRoot } 

procedure  ReducePoly(var  Degree  :  integer; 
var  Poly  :  TNCompVector; 

Root  :  TNcomplex); 


{=  Input:  Degree,  Poly,  Root  = } 

{=  Output:  Degree,  Poly  = } 

{=  =} 

{=  This  procedure  deflates  the  polynomial  Poly  by  = } 

{=  factoring  out  the  Root.  Degree  is  reduced  by  one.  = } 


var 

Term  :  integer; 

NewPoly  :  TNCompVector; 

Dummy  :  TNcomplex; 

begin 

NewPoly[Degree  -  1]  :=  Poly[Degree]; 
for  Term  Degree  -  1  downto  1  do 
begin 

Mult(NewPoly[Term],  Root,  Dummy); 
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Add(Dummy,  Poly[Term],  NewPolylTerm  -  1]); 
end; 

Degree  >  Pred(Degree): 

Poly  >  NewPoly; 
end;  {  procedure  ReducePoly } 

begin  {  procedure  Laguerre  } 

I nitAndTest( Degree,  Poly.  Tol,  Maxiter,  InitGuess.  NumRoots,  Roots, 
yRoots,  Iter,  GuessRoot,  InitDegree,  initPoly,  Error); 
while  (Degree  >  0)  and  (Error  >  0)  do 
begin 

FindOneRoot(Degree,  Poly,  GuessRoot,  Tol.  Maxiter, 
Roots[NumRoots  +  1],  yRoots[NumRoots  +  1], 


lter[NumRoots  +  1],  Error); 
if  Error  =  0  then 
begin 

{=  The  next  statement  refines  the  approximate  root  by  =} 

{=  plugging  it  into  the  original  polynomial.  This  =} 

{=  eliminates  a  lot  of  the  round-off  error  =} 

{=  accumulated  through  many  iterations  =} 


FindOneRoot(lnitDegree,  InitPoly,  Roots(Num Roots  +  I), 

Tol,  Maxiter,  Roots[NumRoots  +  1], 
yRoots(NumRoots  +  1),  Additer,  Error); 
lter(NumRoots  +  1)  ;=  IterfNumRoots  +  1]  +  Additer; 

NumRoots  :=  Succ(NumRoots); 

ReducePoly(Degree,  Poly,  RootsfNumRootsJ);  {  Reduce  polynomial } 
end; 

GuessRoot  Roots[N  urn  Roots]; 
end; 

end;  {  procedure  Laguerre  } 
begin 

end.  { RootsOfEquat } 


90 


unit  Message  (6000); 


[=  This  unit  provides  the  several  messages  which  inform  the  « } 

[=  user  with  some  warnings  and  the  help  informations.  >} 

{$T  APPLFFTD}  {  Set  the  application  type  and  creator } 

{$S+} 

{$!-}  {  Turn  off  I/O  error  checking  } 

{$U  RootsFinder} 

{$U  SpecVar} 

{$U  TurboGraph} 

interface 

uses 

MemTypes.  QuickDraw,  OSIntf,  ToolIntf.Paddntf, 

PasPrinter.SANE.MacPrint,  RootsFinder,  SpecVar, 

{$S  Second  Segment} 

Turtx>Graph; 

procedure  SetUpWindows; 
procedure  MakeInfoScreen; 
procedure  OoAbout; 
procedure  AlertBoxI; 
procedure  AlertBox2; 
procedure  AlertBoxS; 
procedure  InfoGetEQParameter; 
procedure  infoGetCoeff; 
procedure  InfoPlotOneParameter; 
procedure  InfoPlotTwoParameter; 
procedure  InfoPrint; 

implementation 

procedure  SetUpWindows; 

[=  Define  all  of  the  windows  used  in  the  program  =} 
begin 

{  The  initial  information  screen  } 

DefineWindow(lnfoScrnWind,  100,  100,  362,  200,  dBoxProc); 

{  The  real  transform  window  } 

DefineWindow(RootLocusWind,  XMinGlb+5,  YMinGlb-t-43,  XMaxGlb-5,  YMaxGlb-5, 
documentProc); 

DefineHeader(RootLocusWind,  'Root  Lx>cus'); 
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DefineWindow(RegendBox,100,  100,362.150,  dBoxProc); 

{ The  about  box  window } 

DefineWindow(AboutBoxWind,  125,  60.  367,  272,  dBoxProc); 
DefineWindow(AlertBox,50,  100,360,170,  dBoxProc); 

DefineWindow(HelpWind, 20,60, 480.330,dBoxProc); 

end;  { SetUpWindows } 

procedure  MakeInfoScreen; 

{ . . . . . } 

{»  Bring  up  a  window  with  a  description  of  what  this  program  does  =} 

var 

SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  ;  integer;  {  Window  offset  from  top  } 

begin 

GetPort(SavePort);  {  save  the  current  grafport  } 

SelectWind(lnfoScrnWind,  FALSE); 
with  Window[lnfoScrnWindl.P^.portRect  do 
begin 

InLeft  :=  (XMaxGlfa  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  :=  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window(lnfoScrnWind].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility(lnfoScrnWind,  TRUE); 

T  extFont(SystemFont) ; 

TextSize(14); 

TextStyle((|); 

MoveTo(20,  20); 

DrawStringCWelcome  To'); 

TextSize(16); 

MoveTO(70,  50); 

DrawStringCMAC  RootLocus'); 

TextSize(12); 

TextStyle(n); 

MoveTo(90,  80); 

DrawStringCHave  a  fun  11'); 
repeat  until  Button; 

HideWindow(Window(lnfoScrnWlndJ.P); 

GetPort(SavePort);  {  save  the  current  grafport } 
end;  {  MakeInfoScreen } 
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procedure  DoAbout; 

{«  Bring  up  the  about  box  -} 

{SSSSSSKSESCSSSSSSS} 

var 

SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  :  integer;  {  Window  offset  from  top  } 

begin 

GetPort(SavePort);  {  save  the  current  grafport } 

Sel6CtWind(AboutBoxWind,  FALSE); 

with  Window[AboutBoxWind].P'^.portRect  do 

begin 

InLeft  :=  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  ;=  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MgveWindow(Window[AboutBoxWind].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility(AboutBoxWind,  TRUE); 

T  extFont(SystemFont) ; 

TextSize(18); 

TextStyle([]); 

MoveTo(60.  40); 

DrawStringCMAC  RootLocus'); 

TextStyle([l); 

TextSi2e(12): 

MoveTo(85,  60); 

DrawStringCversion  1.00'); 

MoveTo(35,  90); 

DrawStringCComputer  Aided  Desgin  Tool'); 

MoveTo(35,  110); 

DrawStringCFor  RootLocus  Analysis'); 

MoveTo(97,  160); 

DrawString('written  by’); 

TextSlze(13); 

MoveTo(85,  185); 

DrawStringCKO  SUNG  HOON'); 

T  extFont(SystemFont) ; 
repeat  until  Button; 

HideWindow(  Window[AboutBox  Wind]  .P) ; 

SetPort(SavePort); 
end;  { DoAbout } 

procedure  AlertBoxI; 


{ . . . } 

Bring  up  a  window  with  a  description  of  warning  message. 

{ . . . } 
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var 

SavePort  :  GrafPlr; 

InLeft  :  integer;  {  Window  offset  from  left } 
inTop  :  integer;  {  Window  offset  from  top  } 
begin 

GetPort(SavePort);  {  save  the  current  gra^rt } 

SetectWind(3,  FALSE); 

with  Window[3].P''.portRect  do 

begin 

InLeft  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 
InTop  >  (YMaxQlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window[3].P,  InLeft,  InTop.  TRUE);  {  Center  the  window  } 
SetVisibility(3.  TRUE); 

T  extFont(SystemFont) ; 

TextSize(12); 

TextStyledl); 

MoveTo(30,  30); 

DrawStringCThere  is  no  Initial  Degree.'); 
repeat  until  Button; 

HideWindow(Window[3].P); 

end; 


procedure  AlertBox2; 

[=  Bring  up  a  window  with  a  description  of  warning  message.  «) 
var 

SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  :  integer;  {  Window  offset  from  top  } 

begin 

GetPort(SavePort);  {  save  the  current  grafport } 

SelectWind(5,  FALSE); 

with  Window(5].P*.portRect  do 

begin 

InLeft  :«  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 
InTop  :«  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window[5].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility(5.  TRUE); 

T  extFont(SystemFont) ; 

TextSize(12); 

TextStyle{[]); 

MoveTo{30,  30); 
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DrawStringCThere  is  no  Initial  Degree  or*); 

MoveTo(50,50); 

DrawStringCCharacteristic  Equation  Coefficients.'); 
repeat  until  Button; 

HideWindow(Window[5].P) ; 
end; 

procedure  AlertBoxS; 

[a  Bring  up  a  window  with  a  description  of  warning  message.  »} 
var 

SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  :  integer;  {  Window  offset  from  top  } 
begin 

GetPort(SavePort);  {  save  the  current  gra^rt } 

SelectWind(5,  FALSE); 

with  Window(51.P''.portRect  do 

begin 

InLeft  :=  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 
InTop  :=  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window[5].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility{5,  TRUE); 

TextFont(SystemFont) ; 

TextSi2e(12); 

TextStyledJ): 

MoveTo(30,  30); 

DrawStringCYou  have  a  wrong  Input.'); 

MoveTo(50,50): 

DrawStringCCheck  the  Help  Menu  and  try  again.'); 
repeat  until  Button; 

HideWindow(  Window[5]  .P) ; 
end; 

procedure  InfoGetEQParameter; 


{=  Bring  up  the  Help  information  of  EQ  parameter  dialog. 

{= . . } 

var 


SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 
InTop  :  integer;  {  Window  offset  from  top  ) 


begin 

GetPort(SavePort);  {  save  the  current  grafport } 


SelectWind(HelpWind.  FALSE); 

with  Window[HelpWind].P^.portRect  do 

begin 

InLeft  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window(HelpWind].P,  InLeft.  InTop.  TRUE);  {  Center  the  window  } 
SetVisibility(HelpWind.  TRUE); 

TextFont(SystemFont) ; 

TextSize(14); 

TextStyle(n); 

MoveTo(150.  15); 

DrawStrIngfEQ  PARAMETER*); 

TextStyle(Il); 

TextSize(12); 

MoveTo(20.  40); 

DrawString(The  EQ  Parameter  stands  for  equation  parameter.  It  will  allow  you'); 
MoveTo(10.  55); 

DrawStringCto  input  the  degree  of  polynomial  and  some  equation  parameters. '); 
MoveTo(  10.70); 

DrawStringCThe  degree  of  the  polynomial  should  be  1  up  to  10.  Then  there  are  '); 
MoveTo(10,  85); 

DrawStringCthe  default  values  for  the  other  parameters.  These  values  avoid  '); 
MoveTo(10.  100); 

DrawStringCthe  convergence  error  for  almost  ail  polynomials.  But  if  '); 
MoveTo(10.  115); 

DrawStringCconvergence  error  messages  appears  on  the  screen,  you  can  change  '); 
MoveTo(10.  130); 

DrawStringCthese  parameter  values.  These  parameters  must  satisfy  the  '); 
MoveTo(10.  145); 

DrawStringCfollowing  conditions:'); 

MoveTo(30.  160); 

DrawString('(1 )  Initial  guess  >«  0  '); 

MoveTo(30.  175); 

DrawString('(2)  Maximum  Iteration  0  '); 

MoveTo(30,  190); 

DrawString('(3)  Tolerance  >  0  '); 

MoveTo(40.  260); 

DrawStringC"  CLICK  THE  MOUSE  ONCE  TO  RETURN  THE  MAIN  MENU "’); 
repeat  until  Button; 

HideWindow(Window[HelpWind].P); 

SetPort(SavePort); 

end; 

procedure  InfoGetCoeff; 


{-  Bring  up  the  Help  information  of  Get  Coeff  dialog.  -} 
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var 

SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  ;  integer;  {  Window  offset  from  top  ) 

begin 

GetPort(SavePort);  {  save  the  current  grafport } 

SelectWind(HelpWind.  FALSE); 

with  Window[HelpWind].P'^.portRect  do 

begin 

InLeft  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window(HelpWind].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility(HelpWind,  THUE); 

T  extFont(System  Font) ; 

TextSi2e(14); 

TextStyle([]); 

MoveTo(170,  15); 

DrawStringCGET  COEFP); 

TextStyledl); 

TextSize(12); 

MoveTo(20,  40); 

DrawStringCThe  Get  Coeff  stands  for  Get  Coefficients.  It  will  allow  you  to’); 
MoveTo(10,  55); 

OrawStringCinput  the  algebraic  expression  for  the  coefficients  of  the  '); 

MoveTo(  10,70); 

DrawStringCcharacteristic  equation.  It  may  have  up  to  two  undetermined  '); 
MoveTo(10,  85); 

DrawString(‘parameters(A  and  B).  In  case  of  the  one  parameter  root  locus  '); 
MoveTo(10,  100); 

DrawStringCmethod,  you  use  only  one  undetermined  paramerter(A).  The  routine'); 
MoveTo(10,  115); 

DrawStringCuses  standard  algebraic,  or  infix,  notation  with  parenthesis'); 
MoveTo(10,  130); 

DrawStringCallowed.  Operator  can  include  +,  -,  *,  /,  and  ''(expomentiation).'); 
MoveTo(10,  145); 

DrawStringCThe  unary  minus  sign  is  allowed.'); 

MoveTo(10,  160); 

DrawStringCIf  you  choose  the  Get  Coeff  command  without  the  degree  of  the'); 
MoveTo(10,  175); 

DrawString('polynomial,the  message  appears  on  the  screen.  It  tells  you  that'); 
MoveTo(10,  190); 

DrawStringCthe  degree  of  the  polynomial  has  not  been  entered  yet.'); 

MoveTo(40,  260); 

DrawStringC"  CLICK  THE  MOUSE  ONCE  TO  RETURN  THE  MAIN  MENU  "'); 
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repeat  until  Button; 
HideWindow(Window[HelpWind].P); 
SetPort{SavePort): 
end; 

procedure  InfoPlotOneParameter; 


{-  Bring  up  the  Help  information  of  One  parameter  dialog.  »} 


var 

SavePort  ;  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  :  integer;  {  Window  offset  from  top  } 

begin 

GetPort(SavePort);  {  save  the  current  grafport } 

SelectWind(HelpWind,  FALSE); 

with  WindowfHelpWindJ.P-'.portRect  do 

begin 

InLeft  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  :=  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window(HelpWind].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  ) 
SetVisibility(HelpWind,  TRUE); 

T  extFont(System  Font) ; 

TextSize(14): 

TextStyle(n): 

MoveTo(145,  15); 

DrawStringCONE  PARAMETER’); 

TextStyle([]); 

TextSize(12); 

MoveTo(20,  35); 

DrawStringCOne  parameter  command  will  allow  you  to  input  the  plot  data  of); 
MoveTo(10,  50); 

DrawStringCthe  one  parameter  root  locus  method.  The  default  values  are  shown'); 
MoveTo(10,65); 

DrawStringCin  the  dialog  box.  They  can  be  changed  as  desired.'); 

MoveTo(20,  80); 

OrawStringC-.  Enter  the  minimum  and  maximum  gain  values  into  the  Min  Gain'); 
MoveTo(10,  95); 

DrawStringCand  the  Max  Gain  insertion  box.'); 

MoveTo(20,  110); 

DrawStringC-.  Choose  the  type  of  interval.  There  are  two  types  of  interval:'); 
MoveTo(10,  125); 

DrawStringCLinear  and  Logarithmic  interval.  When  you  click  the  radio  button, '); 
MoveTo(10,  140); 

DrawStringCthe  desired  type  of  interval  is  choosen.'); 
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MoveTo(20,  155); 

DrawStringC-.Choose  the  type  of  scale  for  the  axis.  There  are  two  types  of  '); 
MoveTo(10.  170): 

DrawStringCscale:  Auto  and  Manual  scale.  When  you  click  the  radio  button,’); 
MoveTo(10.  185); 

DrawStringCthe  desired  type  of  scale  is  choosen.'); 

MoveTo(20,  200); 

DrawStringC-.Enter  some  number  into  Points  to  Plot  inserton  box  in  order  to  '); 
MoveTo(10,  215): 

Drawstring  ('decide  the  plot  resolution.'); 

MoveTo(40,  260): 

DrawStringC"  CLICK  THE  MOUSE  ONCE  TO  RETURN  THE  MAIN  MENU  -); 
repeat  until  Button; 

HideWindow{Window[HelpWind].P) : 

SetPort{SavePort); 

end: 

procedure  InfoPlotTwoParameter; 


{=  Bring  up  the  Help  information  of  the  two  parameter  dialog, 
var 

SavePort  :  GrafPtr; 

InLeft  ;  integer;  {  Window  offset  from  left } 

InTop  :  integer;  {  Window  offset  from  top  } 

begin 

GetPort(SavePort):  {  save  the  current  grafport  } 

SelectWind(HelpWind.  FALSE): 
with  Window[HelpWindl.P''.portRect  do 
begin 

InLeft  :=  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  :=  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window(HelpWind].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility(HelpWind,  TRUE); 

TextFont(SystemFont): 

TextSize(14); 

TextStyle(n): 

MoveTo(145,  15); 

DrawStringCTWO  PARAMETER'); 

TextStyle([]): 

TextSize(12); 

MoveTo(20,  35): 

DrawString('Two  parameter  command  will  allow  you  to  input  the  plot  data  of); 
MoveTo(10,  50); 

DrawStringCthe  two  parameter  root  locus  method.  The  default  values  are  shown'); 
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MoveTo(  10,65); 

DrawString(’in  the  dialog  box.  They  can  be  changed  as  desired.'); 

MoveTo{20,  80); 

DrawStringCThe  items  in  this  dialog  are  similar  to  those  in  the  one  parameter'); 
MoveTo{10.  95); 

DrawStringCpIot  data  dialog  box.  but  several  items  are  different.  '); 

MoveTo(20,  110); 

DrawString(There  exists  one  more  undetemined  parameter  B  to  be  inserted.The'); 
MoveTo{10.  125); 

DrawStringCHow  many  loci  item  lets  you  decide  how  many  loci  are  to  be  drawn  '); 
MoveTo(10,  140); 

DrawStringCfor  each  parameter.  It  will  be  1  up  to  10  for  each  parameter.'); 
MoveTo(10.  155); 

DrawString(There  is  no  auto  scale  for  axes.  Only  the  manual  scale  is  available'); 
MoveTo(10,  170); 

DrawStringCThe  last  item  is  the  marking  and  justification  in  order  to  draw'); 
MoveTo(10.  185); 

DrawStringCthe  selected  A  and  B  values  on  the  plot.  Two  buttons  are  chosen'); 
MoveTo(10,  200); 

DrawStringCeach  time  for  each  parameter;  one  for  position,  the  other  justifi*  '); 
MoveTo(10,  215); 

Drawstring  ('cation  to  draw.'); 

MoveTo(40.  260); 

DrawStringC  CLICK  THE  MOUSE  ONCE  TO  RETURN  THE  MAIN  MENU  "'); 
repeat  until  Button; 

HideWindow(Window[HelpWind].P); 

SetPort(SavePort): 

end; 

procedure  InfoPrint; 

[=  Bring  up  the  Help  information  of  print  out.  =  } 

var 

SavePort  :  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  :  Integer;  {  Window  offset  from  top  } 

begin 

GetPort(SavePort);  {  save  the  current  grafport } 

Select Wind(Help Wind,  FALSE); 

with  Window[HelpWind],P''.portRect  do 

begin 

InLeft  >  (XMaxGlb  -  (Right-Left))  div  2;  {  Calculate  window  offsets  } 

InTop  :«  (YMaxGlb  -  (Bottom-Top))  div  2; 
end; 

MoveWindow(Window(HelpWind].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
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SetVisibility(HelpWind,  TRUE); 

TextFont(SystemFont): 

TextSize(14): 

TextStyle([l): 

MoveTo(160.  15); 

DrawStringCPRINT  OUT'); 

TextStyle([l); 

TextSlze(12); 

MoveTo(20,  35); 

DrawStringCIn  order  to  print  out  your  work,  there  are  a  few  ways  available.'); 
MoveTo(20.  50); 

DrawStringC-.Select  the  print  command  in  the  File  menu.  This  allows  the  '); 
MoveTo(10,65); 

DrawStringCuser  to  get  a  hard  copy  of  any  plot  displayed.'); 

MoveTo(20,  80); 

DrawStringC-.  Hold  the  command  key  and  then  type  the  number  4  to  print  ’); 
MoveTo(10,  95); 

DrawStringCthe  contents  of  the  active  window  immediately.'); 

MoveTo(20,  110); 

DrawStringC*.  Press  the  command  and  shift  key  and  type  the  number  3  to  '); 
MoveTo(10.  125); 

Drawstring  ('create  a  MacPaint  document.'); 

MoveTo(40.  260); 

DrawStringC"  CLICK  THE  MOUSE  ONCE  TO  RETURN  THE  MAIN  MENU  "'); 
repeat  until  Button; 

HideWindow(Window[HelpWind].P); 

SetPort(SavePort); 

end; 

begin 

end.{unit  message} 
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unit  My  Dialog  (7000); 


This  unit  supports  four  dialog  boxes  for  input  data. 

{  Set  the  bundle  bit } 

{$R  MacRootLocus.Rsrc}  {  Identify  resource  file  for  menu  and  icon  info  } 
{$T  APPLFFTD}  {  Set  the  application  type  and  creator } 

{$S+}  { generate  segmented  code} 

{$!-}  {  Turn  off  I/O  error  checking  } 

{$U  SpecVar} 

{$U  RootsFinder) 

{$U  MakeRoot} 

{$U  TurboGraph} 

{$U  Message} 

interface 

uses 

MemTypes,  QuickDraw,  OSIntf,  ToolIntf.PackIntf, 
PasPrinter.SANE.MacPrint,  RootsFinder,  SpecVar, 

{$S  Second  Segment} 

TurboGraph, 

Message, 

{$S  Third  Segment} 

MakeRoot; 


procedure  MakeRegend; 
procedure  GetEQParameter; 
procedure  GetCoeff; 
procedure  PlotOneParameter; 
procedure  PLotTwoParameter; 

implementation 

procedure  MakeRegend; 

This  procedure  provides  the  informaton  box  with 
the  capabilety  of  the  text  editor 

var 

SavePort  ;  GrafPtr; 

InLeft  :  integer;  {  Window  offset  from  left } 

InTop  :  integer;  {  Window  offset  from  top  } 
begin 

GetPort(SavePort);  {  save  the  current  gra^rt  } 
SelectWindO,  FALSE); 
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Textsizei^); 

with  Window[3].P''.portRect  do 
begin 

InLeft  >  Round((XMaxGlb  -  (Right-Left))  /  1.38);  {  Calculate  window  offsets  } 
InTop  Round((YMaxGlb  -  (Bottom-Top))  /1 .05); 
end; 

MoveWindow(Window[3].P,  InLeft,  InTop,  TRUE);  {  Center  the  window  } 
SetVisibility(3,  TRUE); 
txRect  thePort'^.portRect; 
textH  :«  TENew(txRect,  txRect); 

TEIdle(textH); 

TextInputEnabled  >  True; 
end; 


procedure  GetEQParameter; 

{=  This  procedure  provides  *^e  dialog  box  for  the  = } 

{=  equation  parameter  of  the  characteristic  equation.  >) 


const 

CancelBtn  =  1 ; 

OKBtn  =  2; 

InitDegreeText  =  3; 

InitGuessText  « 4; 

MaxIterText  =  5; 

ToleranceText  =  6; 
var 

tempdegree,templter ;  Longint; 
tempGuess :  TNcomplex; 
tempTolerance  :  Extended; 

BoxPlot :  Boolean; 
begin 

theDialog  GetNewDialog(256,  nil,  pointer(-l)); 
GetDltem(theDialog,  3,  theType,  h3,  r); 
GetDltem(theDialog,  4,  theType,  h4,  r); 
GetDltem(theDialog,  5,  theType,  h5,  r); 
GetDltem(theDialog,  6,  theType,  h6,  r); 

Done :«  False; 
tempGuess.Re  >  1.0; 
tempGuess.lm  :«  0.0; 
tempiter  100; 
tempTolerance  >  IE-7; 

BoxPlot  ;*  False; 
repeat 

ModalDialog(nil,  itemHit); 
case  itemHit  of 
CancelBtn  :  done  :«  True; 
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OKBtn 


:  begin 

InitOegtee  :« tempdegree; 

InitGuess.Re  >  tempGuess.Re; 

InitGuess.im  >  tempGuess.lm; 

Maxiter  :>  tempiter; 

Tolerance :» tempTolerance; 
done  :=  True; 
end; 

InitDegreeText  :  begin 

GetlText(h3,  s); 

StringToNum(s,  tempdegree); 

InitDegreeStatus  >  True; 
if  (tempdegree  <1  )  or  (tempdegree  >10)  then 
begin 

sysbeep(lO); 

end; 

end; 

InitGuessText  :  begin 

GetlText(h4,  s); 
tempGuess.Re  :=  Str2Num(s); 
if  (tempGuess.Re  <0 )  then 
begin 

sysbeep(IO); 

end; 

end; 

MaxIterText  :  begin 

GetlText(h5,  s); 

StringToNum(s,  tempiter); 
if  (tempiter  <0  )  then 
begin 

sysbeep(10}; 

end; 

end; 

ToleranceText  '.begin 

GetlText(h6,  s); 
tempTolerance  :=  Str2Num(s); 

If  (tempTolerance  <=  0 )  then 
begin 

sysbeep(IO); 

end; 

end; 

end; 

until  done; 

DisposDialog(theDialog); 

end; 

procedure  GetCoeff; 
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[=  This  procedure  provides  the  dialog  tx>x  for  the  -} 
[=>  coefficients  of  the  characteristic  equation.  -} 

type  char_set  -  set  of  char; 
const 

CancelBtn  - 1 ; 

OKBtn  =  2: 

SOText  -  3; 

SIText  -  4; 

S2Text  =  5; 

S3Text  =  6; 

S4Text  =  7; 

SSText  ■  8; 

S6Text  =  9; 

STText  =  10; 

SSText  =11; 

SSText  =12; 

SIOText  =13; 


ldno,tempitemhit,sssindex,ssslength 
:  Integer; 

h  :  Arrayfl.. 13]  of  Handle: 
set_valid  :char_set; 
op_set  :  char_set; 
sss:  str255; 
sssc:  char; 

noerrorcheck.leftparen:  boolean; 

Procedure  skipblank  ; 
begin 

while  ((sss[sssindex]  =  ' ')  and  (sssindex  <=  ssslength))  do 
sssindex  ;=  sssindex  +  1 ; 

end; 

Procedure  Checknext  (  setofnext1,setofnext2;char_set); 

var  opchar;  char; 

begin 

opchar  :*  sssc; 
sssindex  :=  sssindex  +1 ; 
sssc  ;=  sss[sssindex]; 
if  (sssc  =  ' ')  then  begin 
skipblank; 

sssc sss[sssindex]; 
if  sssindex  <=  ssslength  then  begin 
if  not  (sssc  in  setofnextl)  then 
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noerrorcheck:-  False; 
end  else 

if  (opchar  in  (op_set  +  [’(']))  then 
noerrorcheck false 

end  else 

if  not  (sssc  in  setofnext2)  then 
noerrorcheck:-  false 

end; 


Procedure  CheckExpression(matchparen:boolean); 
begin 

leftparen:-  true; 

while  ((noerrorcheck  and  leftparen)  and  (sssindex  <  ssslength))  do 
case  (sssc)  of 

Checknext(op_set+['.',Tl.op_set+ro*..‘9’;.'/)*l); 

'a'.'b':  Checknext(op_set+  (')'].op_set+['a’,’b’, ')’]); 

Checknext([*0*..'9’,'a’.*b',*('],r0‘..’9*, ‘a’, ‘b’, •(■]); 
Checknext(('0’..’9'],[’0’..‘9']); 

'(';begin 

Checknext(C0’..'9' .'a', •b';-\Tl.t'0*./9'.'a'.'b', •-',•(■]); 
CheckExpression  (false); 
if  sssc  =  ')'  ttien 

if  sssindex  <  ssslength  then  begin 
sssc sss[sssindexl; 

if  sss(sssindex  1]  -  '  *  then  skipblank; 
if  sssindex  <  ssslength  then 

Checknext(op_set+[')'I,op_set+(')’]) 
else  if  not  matchparen  then  noerrorched^ false 
end  else  begin 

if  not  matchparen  then  noerrorcheck  false 
end 

else  noerrorcheck  false; 
end; 

')':begin 

if  matchparen  then  noerrorcheck;-  false; 
leftparen  false; 
end; 

'  ':begin 

skipblank ; 

if  sssindex  -  ssslength  then  begin 

if  not  (sss[sssindex]  in  ['0'..'9','a','bl)  then 
noerrorcheck:-  false 

end  else  if  sssindex  <  ssslength  then  sssc  :«  sss[sssindex]; 
end 
end; 


if  not  noerrorcheck  then  sysbeep  (30); 
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end; 

{begin 

if  not  InitDegreeStatus  then 
begin 

SysBeep(30); 

AlertBoxt ; 
end 
else} 
begin 

set_valid  :«  r0’..'9'.'a',’b*.’+'.'**.7*.*-*; 
op_set  := 

if  not  InitDegreeStatus  then 
begin 

SysBeep(30); 

AlertBoxt; 

end 

else 


begin 

Idno  300  +  InitDegree; 

theDialog  :=  GetNewDialog(ldno,  nil.  Pointer  (-1)); 
for  n  ;=  3  to  InitDegree  +  3  do 
GetDltem(theDialog,  n,  theType,  h[nl,  r); 
done  :=  False; 

ModalDialog(nil,  itemHit); 


repeat 

tempitemhit  :=  itemHit; 
case  itemHit  of 
CancelBtn  :  done  :=  True; 
OKBtn  :  done  :=  True; 


3..12  :  begin 

while  tempitemhit  =  itemHit  do 
ModalDialog(nil,  tempitemhit); 

GetIText  (h{itemHit].sss); 
noerrorcheck  :>  true; 
sssindexta  1 ; 
ssslength Length  (sss); 
sssc sss[11; 

if  not  (sss(11  in  ['0'..'9’, ’a’, 'b', then 
noerrorcheck  :>  False; 

CheckExpression  ( True); 

GetlText(h[itemHit],lnfixArray(lnitDegree  +  3  -  itemHit]); 
GetCoeffStatus True; 
itemhit tempitemhit; 
if  itemhit  >  OKBtn  then  done  >  true 
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end; 


end; 

until  done; 

DisposDialog(theDialog) ; 
end; 
end; 

{end;} 

procedure  PlotOneParameter; 

{=  This  procedure  provides  the  dialog  box  for  the  «} 

{=  plot  data  of  the  one  parameter  root  locus  method.  -} 

const 

CancelBtn  =  1 ; 

OKBtn  =  2; 

UnearBtn  « 3; 

LogarithmicBtn  =  4; 

AutoBtn  s  5; 

ManualBtn  =  6; 

AMinGainText  =  7; 

AMaxQainText  « 8; 

XMinText  =  9; 

XMaxText  « 10; 

YMinText  *11; 

YMaxText  =  12; 

PlotPointsText  =  13; 


saveSoundVol :  Integer; 

radButton  :  array  [3..6]  of  ControlHandle; 

h  :  Array(7..131  of  Handle; 

tempAMinGain  :  Extended; 

tempAMaxGain  :  Extended; 

tempXMin  :  Extended; 

tempXMax  :  Extended; 

tempYMin  :  Extended; 

tempYMax  ;  Extended; 

tempPoints  :  Longint; 

OkPlot  iboolean; 
begin 

if  not  InitDegreeStatus  or  not  GetCoeffStatus  then 
begin 

SysBeep(30}; 

AlertBox2; 

end 

else 
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begin 

GetSoundVol  (saveSoundVol); 

SetSoundVol  (1): 

theDialog  >  GetNewDialoQ(257,  nil,  Pointer  (-1)): 
for  n  >  UnearBtn  to  ManualBtn  do 

GetDltem(theDialog,n.theType,Handle(radButton(n]).r): 

SetCtIValue  (radButton[UnearBtn],1 ); 
linear  :*  True: 

SetCtIValue  (radButton[AutoBtn],1 ); 

AutoScale True; 
for  n  7  to  13  do 

GetDlten){theDialog,  n,  theType,  hin],  r); 
tempAMinGain  :■  Str2Num('0.1'); 
tempAMaxGain  :■  Str2Num('10000'); 
tempXMin  :=  Str2Num(’-10'); 
tempXMax  :=  Str2Num('5'): 
tempYMin  :=  Slr2Num('-10'): 
tempYMax  :=  Str2Num('10'): 

StringToNum('50',  tempPoints); 
done  :=  False; 

OkPlot :»  False; 
repeat 

ModalDialog{nil,itemHit); 
case  item  Hit  of 
CancelBtn  :  done  :»  True: 

OKBtn  :  begin 

done  :=  True; 

OkPlot  ;=  True; 
end; 

UnearBtn  ;  begin 

SetSoundVol  (1): 

for  n  :«  UnearBtn  to  LogarithmicBtn  do 

SetCtlValue(radButton[n].Ord{n  *=  itemHit)); 
linear  :=  True; 
end; 

LogarithmicBtn  :  begin 

SetSoundVol  (1); 

for  n  UnearBtn  to  LogarithmicBtn  do 

SetCtlValue(radButton[nl,Ord(n  -  itemHit)): 
linear False; 
end; 

AutoBtn  :  begin 

SetSoundVol  (1); 

for  n  >  AutoBtn  to  ManualBtn  do 

SetCtlValue(radButton(nl,Ord{n  -  itemHit)): 
AutoScale True; 


109 


end; 


ManualBtn  :  begin 

SetSoundVol  (1); 

for  n  >  AutoBtn  to  ManualBtn' do 

SetCtlValue(radButton[n],Ord(n  -  itemHit)); 
AutoScale  >  False; 
end; 

AMinGainText  :  begin 

GetiText(h[itemHit] .  s); 
tempAMinGain  :«  Str2Num(s); 
if  (tempAMinGain  <  0  )  or  (tempAMinGain  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

AMaxGainText  :  begin 

GetlText(h[ilemHit] ,  s); 
tempAMaxGain  :=  Str2Num(s); 
if  (tempAMaxGain  <  0  )  or  (tempAMaxGain  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

XMinText  -.begin 

GetlText(h(itemHit] ,  s); 
tempXMin  >  Str2Num(s); 
if  (tempXMin  <-1e7  )  or  (tempXMin  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

XMaxText  :  begin 

GetiText(h[itemHit] ,  s); 
tempXMax  >  Str2Num(s); 
if  (tempXMax  <-1e7  )  or  (tempXMax  >  1e7)  then 
begin 

sysbeep(10); 

end; 

end; 

YMinText  :  begin 

GetlText(h[itemHit] ,  s); 
tempYMin  >  Str2Num(s); 
if  (tempYMin  <-1e7  )  or  (tempYMin  >  1e7)  then 
begin 

sysbeep(10); 

end; 
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YMaxText 


end; 

:  begin 

GetlText(h[itemHit] ,  s); 
tempYMax  >  Str2Num(s); 
if  (tempYMax  <-1e7  )  or  (tempYMax  >  1e7)  then 
begin 

sysbeep(10); 
end; 
end; 

PiotPointsText  :  begin 

GetlText(h[itemHit] ,  s); 

StringToNum(s,  tempPoints); 
if  (tempPoints  <-  0)  or  (tempPoints  >  150)  then 
begin 

sysbeep(IO); 

end; 

end; 

end;{case  end} 
untii  done; 

DisposDiaiog(theDialog) ; 

if  OkPiot  then 

begin 

AMinGain  :=  tempAMinGain; 

AMaxGain  >  tempAMaxGain; 

XMn  :=  tempXMin; 

XMx  tempXMax; 

YMn  :=  tempYMin; 

YMx  :=  tempYMax; 

Points  >  tempPoints; 

HideCursor; 

Rewrite(OutFiie,‘RootLocus1.data'); 

GetRootl ; 

Close(OutFile); 

ShowCursor; 

MakeRegend; 

end;{if} 

end;{case} 

e  nd ;{ PlotOneParameter) 


procedure  PLotTwoParameter; 

{s  This  procedure  provides  the  dialog  box  for  the  =} 

{s  plot  data  of  the  two  parameter  root  loois  method. 

const 

CancelBtn  -  1 ; 
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OKBtn  « 2; 

LinearBtn  « 3; 

LogarithmicBtn  «  4; 

AMarkStBtn  -  5; 

AMarkEdBtn  « 6; 

AJustifyRBtn  =  7; 

AJustifyLBtn  « 8; 

BMarkStBtn  « 9; 

BMarkEdBtn  » 10; 

BJustifyRBtn  =11; 

BJustifyLBtn  =12; 

AMinGainText  =13; 

AMaxGainText  =14; 

BMinGainText  =15; 

BManGainText  =16; 

LociText  =17; 

PlotPointsText  =18; 

XMinText  =  19; 

XMaxText  =  20; 

YMinText  =  21 ; 

YMaxTexl  =  22; 

var 

saveSoundVol.i ;  Integer; 

radButton  :  array  (3. .121  of  ControlHandle; 

h  :  Array[13..22]  of  Handle; 

tempAMinGain  :  Extended; 

tempAMaxGain  :  Extended; 

tempBMinGain  :  Extended; 

tempBMaxGain  :  Extended; 

tempXMin  :  Extended; 

tempXMax  :  Extended; 

tempYMin  ;  Extended; 

tempYMax  ;  Extended; 

tempStep,  tempPoints  :  Longint; 

OkPiot ;  Boolean; 
begin 

if  not  InitDegreeStatus  or  not  GetCoeffStatus  then 
begin 

SysBeep(30); 

Alert  Box2; 
end 
else 
begin 

GetSoundVol  (saveSoundVol); 

SetSoundVol  (1); 

theDialog  GetNewDialog(258,  nil,  Pointer  (-1)); 
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for  n  :«  LinearBtn  to  BJustifyLBtn  do 
GetDltem(theOialog,n.theType.Handle(radButton[n]),r); 
SetCtI Value  (radButton[LJnearBtn],1 ); 
linear  True; 

SetCtIValue  (radButton[AMarkStBtn],1 ); 

AMarkStatus  False; 

SetCtIValue  (radButton[AJustifyRBtn].1 ); 

ARJustification  >  True; 

SetCtIValue  (radButton[BMarkStBtn],1); 

BMarkStatus  >  False; 

SetCtIValue  (radButton[BJustifyRBtn].1 ); 

BR Justification  >  True; 
for  n  >  13  to  22  do 

GetDltem(theDialog,  n.  theType,  h[nl,  r); 
tempAMinGain  :=  Str2Num('0.1'); 
tempAMaxGain  ;=  Str2Num(’10000'); 
tempBMinGain  >  Str2Num('0.1'); 
tempBMaxGain  Str2Num(‘10000'); 
tempXMin  >  Str2Num('-10'); 
tempXMax  :=  Str2Num(’5'); 
tempYMin  >  Str2Num('-10'); 
tempYMax  :=  Str2Num('10'); 

StringToNum('4',  tempStep); 

StringToNum('50',  tempPoints); 
done  >  False; 

OkPlot False; 
repeat 

ModalOialog(nil,ltemHit); 
case  itemHit  of 
CancelBtn  :  done  :=  True; 

OKBtn  :  begin 

done  :=  True; 

OkPlot  :=  True; 
end; 

LinearBtn  :  begin 

SetSoundVol  (1 ); 

for  n  >  LinearBtn  to  LogarithmicBtn  do 
SetCtlValue(radButton[n],Ord(n  -  ItemHit)); 
linear :«  True; 
end; 

LogarithmicBtn  :  begin 

SetSoundVol  (1); 

for  n  LinearBtn  to  LogarithmicBtn  do 
SetCtlValue(radButton[n],Ord(n  -  itemHit)); 
linear  >  False; 
end; 
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AMarkStBtn  :  begin 

SetSoundVol  (1); 

for  n  >  AMarkStBtn  to  AMarkEdBtn  do 

SetCtlValue(radButton[n].Ord(n  itemHit)); 
AMarkStatus  >  False; 
end; 

AMarkEdBtn  :  begin 

SetSoundVol  (1); 

for  n  >  AMarkStBtn  to  AMarkEdBtn  do 

SetCtlValue(radButtonIn],Ord(n  « itemHit)); 
AMarkStatus  >  True; 
end; 

AJustifyRBtn  ;  begin 

SetSoundVol  (1); 

for  n  ;=  AJustifyRBtn  to  AJusfrfyLBtn  do 
SetCtlValue(radButton[n],Ord(n  =  itemHit)); 
ARJustification  ;=  True; 
end; 

AJustifyLBtn  :  begin 

SetSoundVol  (1); 

for  n  :=  AJustifyRBtn  to  AJustifyLBtn  do 
SetCtlValue(radButton[n],Ord(n  =  itemHit)); 
ARJustification  ;»  False; 
end; 

BMarkStBtn  :  begin 

SetSoundVol  (1); 

for  n  :=  BMarkStBtn  to  BMarkEdBtn  do 
SetCtlValue(radButton[n],Ord(n  =  itemHit)); 
BMarkStatus False; 
end; 

BMarkEdBtn  :  begin 

SetSoundVol  (1 ); 

for  n  BMarkStBtn  to  BMarkEdBtn  do 

SetCtlValue(radButton[n],Ord(n  -  itemHit)); 
BMarkStatus :«  True; 
end; 

BJustifyRBtn  :  begin 

SetSoundVol  (1); 

for  n  BJustifyRBtn  to  BJustifyLBtn  do 
SetCtlValue(radButton[n],Ord(n  >  itemHit)); 
BR Justification True; 
end; 

BJustifyLBtn  :  begin 
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SetSoundVol  (1); 

for  n  :«  BJustifyRBtn  to  BJustifyLBtn  do 
SetCtlValue(radButton[n].Ord{n  -  itemHit)); 
BRJustification  >  False; 
end; 

AMinGainText  :  begin 

GetlText(h[itemHit],s) ; 
tempAMinGain  >  Str2Num(s); 
if  (tempAMinGain  <  0  )  or  (tempAMinGain  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

AMaxGainText  :  begin 

GetlText(h[itemHitJ,s); 
tempAMaxGain  :=  Str2Num(s); 
if  (tempAMaxGain  <  0  )  or  (tempAMaxGain  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

BMinGainText  :  begin 

GetlText(h(itemHit],s); 
tempBMinGain  ;=  Str2Num(s); 
if  (tempBMinGain  <  0  )  or  (tempBMinGain  >  1e7)  then 
begin 

sysbeep(lO); 

end; 

end; 

BManGainlext  :  begin 

GetlText(h(itemHit],s); 
tempBMaxGain  ;=  Str2Num(s); 
if  (tempBMaxGain  <  0  )  or  (tempBMaxGain  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

XMinText  :  begin 

GetlText(h[itemHit],s) ; 
tempXMin  >  Str2Num(s); 
if  (tempXMin  <-1e7  )  or  (tempXMin  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

XMaxText  :  begin 
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GetIT  ext(h[itemHitl,s) ; 
tempXMax  >  Str2Num(s); 
if  (tempXMax  <-1e7  )  or  (tempXMax  >  leT)  then 
begin 

sysbeep(IO); 

end; 

end; 

YMinText  :  begin 

GetlText(h[itemHltl,s); 
tempYMin  :=  Str2Num(s); 
if  (tempYMin  <-1e7  )  or  (tempYMin  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

YMaxText  ;  begin 

GetlText(h[itemHit],s); 
tempYMax  :=  Str2Num(s); 
if  (tempYMax  <-1e7  )  or  (tempYMax  >  1e7)  then 
begin 

sysbeep(IO); 

end; 

end; 

LociText  ;  begin 

GetlText(h(itemHjti,s); 

StringToNum(s,  tempStep); 
if  (tempStep  <=  0)  or  (tempStep  >  15)  then 
begin 

sysbeep(IO); 

end; 

end; 

PIotPointsText :  begin 

GetIT  ext(h[itemHit],s) ; 

StringToNum(s,  tempPoints); 
if  (tempPoints  <=  0)  or  (tempPoints  >  150)  then 
begin 

sysbeep(IO); 

end; 

end; 

end;{case  end} 
until  done; 

OisposDialog(theDialog) ; 

if  OkPlot  then 

begin 

AMinGain  tempAMinGain; 

AMaxGain  :>  tempAMaxGain; 

BMinGain  tempBMinGain; 
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BMaxGain  :=  tempBMaxGain: 

XMn  >  tempXMin; 

XMx  :=  tempXMax; 

YMn  :=  tempYMin; 

YMx  tempYMax; 

Step  :*  tempStep: 

Points  >  tempPoints; 

HideCursor; 

Rewrite{OutFile,'RootLocus2.data'): 

GetRoot2; 

Close(OutFile); 

ShowCursor; 

MakeRegend: 

end:{if} 

end:{case} 

end;{PlotTwoParameter} 

begin 

end.{unit  MyDialog} 


unit  TurboGraph(1 1000); 


{= 

{=  Turbo  Pascal  Numerical  Methods  Toolbox 
{=  Copyright  (C)  1987  Borland  International 

{= 

{=  This  unit  provides  routines  for  displaying  graphics. 

{$U  RootsFinder} 

{$U  SpecVar} 

interface 

uses 

MemTypes,  QuickDraw,  OSIntf,  Toolintf,  PackIntf.PasPrinter,  SANE,  MacPrint  , 
RootsFinder,  SpecVar; 

const 

MaxWorldsGIb  =10;  {  The  maximum  number  of  worlds  that  can  be  defined  } 
MaxWindowsGIb  =  1 0;  {  The  maximum  number  of  windows  that  can  be  defined  } 
{MaxPlotGIb  =  1500;}{  The  maximum  number  of  points  in  a  plot  array  } 

XI  Offset  =40;  {  The  upper  left  corner  of  the  axis  } 

Y1  Offset  =10; 

X20ffset  =  30;  {  The  lower  right  corner  of  the  axis  } 

Y20ffset  =  90; 

type 

WrkString  =  Str255;  {  A  general  string  type  } 

WorldType  =  record  {  Used  to  store  world  coordinates  } 

XI,  Y1,X2,  Y2  :real; 
end; 

WindowType  =  record  {  Used  to  store  window  information  } 

XI,  Y1,  X2,  Y2  :  integer;  {  The  windows  screen  coordinates  } 

Header :  WrkString;  {  Header  for  a  window } 

W  :  WindowReoord;  {  Mac  window  record } 

P  :  WindowRr;  {  Mac  window  pointer } 

H  :  PicHandle;  { A  handle  to  a  picture } 
end; 

Worlds  «  array[1.. MaxWorldsGIb]  of  WorldType;  {  Holds  world  info  } 
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Winciows  n  array[1  ..MaxWindowsGIb]  of  WindowType;{  Holds  window  info  } 
{  PlotArray  -  array[1  ..MaxPlotGIb,  1..21  of  Extended;} 

var 

{  Coordinates  of  the  currently  active  window } 

XMinGIb,  YMinGIb,  XMaxGlb,  YMaxGlb  :  integer; 

{  Coordinates  of  the  currently  selected  world  } 

XIWIdGIb,  X2WldGlb.  YIWIdGIb,  Y2WldGlb  :  real; 

{  World  coordinate  scaling  factors  } 

AxGlb,  AyGIb,  BxGlb,  ByGIb  :  real; 

{  Coordinates  of  an  axis  if  one  is  defined } 

XIRefGIb,  X2RefGlb,  YIRefGIb,  Y2RefGlb  :  integer; 

{  The  maximum  world  and  window  number  defined  } 

MaxWorldGIb,  MaxWindowGIb  :  integer; 

{  The  currently  selected  world  and  window  } 

WorldNdxGlb,  WindowNdxGlb  :  integer; 

{  Aspect  ratio  for  a  true  circle  } 

AspectFactor :  real; 

AspectGIb  :  real; 

{  Flags  if  an  axis  is  defined  } 

AxisGib :  boolean; 

{  Flags  if  hatching  is  turned  on  } 

HatchGIb :  boolean; 

{  Flags  if  clipping  should  be  performed  } 

ClippingGIb  ;  boolean; 

{  The  currently  selected  line  style  } 

LineStyleGIb  :  integer; 

{  The  current  foreground  color } 

ForeColorGIb  :  integer; 

{  Holds  the  worlds  defined  by  the  user } 

World  :  Worlds; 

{  Holds  the  windows  defined  by  the  user } 

Window  :  Windows; 
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procedure  Error(ProcName  :  WrkString); 

{  Report  that  an  error  occurred  in  the  procedure  named  "ProcName* } 

procedure  SetForegroundColor(Color  :  integer); 

{  Set  the  foreground  drawing  color } 

procedure  SetBackgroundCoior(Color  ;  integer); 

{  Set  the  background  color } 

procedure  DP(X,  Y  ;  integer); 

{  Plot  a  pixel  at  position  (X,  Y) } 

function  PD(X,  Y  :  integer)  :  boolean; 

{  Return  true  if  the  color  of  the  pixel  at  (X,  Y)  matches  ForeColorGIb  } 

procedure  DrawStraight(X1 ,  X2,  Y  :  integer); 

{  Draw  a  horizontal  line  from  XI  ,Y  to  X2,Y  } 

procedure  lnvertWindow(WinNum  :  integer); 

{  Invert  the  window  referenced  by  WinNum  } 

function  RealToStr(R  :  real)  :  WrkString; 

{  Returns  the  string  representation  of  the  real  R.  } 

function  lntToStr{l  :  Integer)  :  WrkString; 

{  Returns  the  string  representation  of  the  integer  I.  } 

procedure  SetLineStyle(Style  :  integer); 

{  Select  the  current  line  style  } 

function  GetLineStyle  :  integer; 

{  Returns  the  current  line  style  } 

procedure  DefineWorld(WorldNum  :  integer;  X_1,  Y_1,  X_2,  Y_2  :  real); 

{  Defines  a  world  coordinate  system  with  a  specific  number } 

procedure  SelectWorid(WorldNum  :  integer); 

{  Select  the  world  associated  with  WorldNum  } 

procedure  DefineWindow(WinNum,  XLo,  YLo,  XHi,  YHi,  WindowType  :  integer); 
{  Defines  a  window  with  a  specific  window  number  and  window  type  } 

procedure  SetVisibility(WinNum  :  integer;  Visibte  :  boolean); 

{  Sets  the  visibility  of  a  window  to  either  TRUE  or  FALSE  } 

procedure  SelectWind( WinNum  :  integer;  Visible  :  boolean); 

{  Selects  a  window  as  visible  or  invisible  } 
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procedure  DefineHeader(WinNum  :  integer;  Hdr  :  WrkString); 

{  Define  a  header  for  a  window  } 

procedure  RemoveHeader(WinNum  :  integer); 

{  Clears  a  header  from  a  window } 

procedure  ReDefineWindow(WinNum,  XI,  Y1.  X2,  Y2,  WindowType  :  integer); 
{  Redefines  the  coordinates  of  an  existing  window } 

function  WindowX(X  :  real)  :  integer; 

{  Converts  an  X  world  coordinate  into  a  screen  coordinate  } 

function  WindowY(Y  :  real)  :  integer; 

{  Converts  a  Y  world  coordinate  into  a  screen  coordinate  } 

procedure  ResetWorlds; 

{  Resets  all  worlds  to  the  maximum  dimensions  of  the  screen  } 

procedure  InitGraphic; 

{  Initializes  the  graphics  system  } 

function  Clip(var  XI,  Y1,  X2,  Y2  :  integer)  :  txx>tean; 

{  Clips  a  line  to  the  current  window  and  returns  TRUE  if  any  part } 

{  of  the  line  is  still  in  the  window  after  the  clip  operation  } 

procedure  DrawPoint(Xr,  Yr  :  real); 

{  Draw  a  point  in  world  coordinates  } 

function  PointDrawn(Xr,  Yr  :  real)  :  tx>olean; 

{  Returns  TRUE  if  the  point  at  (Xr,  Yr)  is  drawn  } 

procedure  DrawLine(X1,  Y1,  X2,  Y2  ;  real); 

{  Draw  a  line  in  world  coordinates  } 

procedure  DrawLineClipped  (XI,  Y1,  X2,  Y2  :  integer); 

{  Draw  a  line  clipped  in  screen  coordinates  } 

procedure  DrawCrossDiag(X,  Y,  Scale  :  integer); 

{  Draw  a  cross  at  screen  coordinate  (X,  Y)  with  a  scaling  factor } 

procedure  DrawWye(X,  Y,  Scale  :  integer); 

{  Draw  a  Y  at  screen  coordinate  (X,  Y)  with  a  scaling  factor } 

procedure  DrawDiamond(X,  Y,  Scale  :  integer); 

{  Draw  a  cross  at  screen  coordinate  (X,  Y)  with  a  scaling  factor } 
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procedure  DrawCircleDirect(Xr.  Yr.  R  :  integer): 

{  Draw  a  circle  at  saeen  coordinate  (Xr,  Yr)  with  a  radius  R  } 

procedure  DrawCircle(X_R,  Y_R,  Xradius  :  real); 

{  Draw  a  circle  at  world  coordinate  (X_R,  Y_R)  with  a  radius  R  } 

procedure  DrawCross(X1 ,  Y1.  Scale  :  integer); 

{  Draw  a  cross  at  screen  coordinate  (X,  Y)  with  a  scaling  factor } 

procedure  DrawStar(X,  Y,  Scale  :  integer); 

{  Draw  a  star  at  screen  coordinate  (X,  Y)  with  a  scaling  factor } 

procedure  DrawSquareC(X1 ,  Y1,  X2.  Y2  :  integer;  Fill  :  boolean); 

{  Draw  a  square  in  screen  coordinates  with  optional  filling  } 

procedure  DrawSquare(X1 ,  Y1,  X2,  Y2  :  real;  Fill  :  boolean); 

{  Draw  a  square  in  world  coordinates  with  optional  filling  } 

procedure  DrawAscii(X,  Y  :  integer;  Size.  CharByte  :  byte); 

{  Draw  a  character  with  ASCII  code  CharByte  } 

procedure  DrawText(X,  Y,  Scale  ;  integer;  Txt  :  Str255); 

{  Draw  a  string  at  screen  coordinate  (X.  Y)  with  a  scaling  factor } 

procedure  DrawTextW{X,  Y  :  real;  Scale  :  Integer;  Txt  :  WrkString); 

{  Draw  a  string  at  world  coordinate  (X,  Y)  with  a  scaling  factor } 

procedure  TextStyle(Face  :  Style); 

{  Face  B  (bold,  italic,  underline,  outline,  shadow,  condense,  extend) } 
procedure  HardCopy(TopWin  :  boolean); 

{  Do  a  screen  dump  of  the  currently  selected  window  (TopWin  =  TRUE) } 
{  or  the  whole  screen  (TopWin  -  FALSE) } 

procedure  OpenPic(WinNum  :  integer;  ShowPic  :  boolean); 

{  Open  a  picture  for  a  specific  window  and  only  show  the  drawing  } 

{ if  ShowPIc  is  set  to  TRUE  ) 

procedure  DrawPic(WinNum  :  integer); 

{  Draw  the  picture  associated  with  a  window ) 

procedure  ErasePic(WinNum  ;  integer); 

{  Erase  the  picture  associated  with  a  window } 

procedure  ClearWindow(WinNum  :  integer); 

{  Clear  the  content  portion  of  a  window  } 
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function  WhereX  :  integer; 

{  Returns  the  X  cursor  position  } 

function  WhereY  :  integer; 

{  Returns  the  Y  cursor  position  } 

procedure  SetWindow(X1 ,  Y1 ,  X2,  Y2  :  integer); 

{  Creates  an  invisible  window  inset  from  the  currently  selected  one  } 

procedure  FindWorld(l  :  integer;  A  :  PlotArray;  NPoints  :  integer); 

{  Finds  a  world  to  fit  a  polygon  defined  in  A } 

procedure  FindWorld1(l  :  integer;  XMn.YMn.XMx.YMxrExTended  ); 

{  Finds  a  world  to  fit  a  polygon  defined  in  A } 

procedure  DrawAxis(Footer1 ,  Footer2  ;  WrkString;  Arrows  :  boolean); 
{  Draws  an  axis  with  Footers  and  optional  arrows  on  the  axis  } 

procedure  ResetAxis; 

{  Sets  AxisGib  to  TRUE } 

procedure  DrawPoiygon(A  ;  PlotArray;  First,  NPoints,  Line,  Scale, 
Lines  :  integer;  CrossHairs  :  boolean); 

{  Draws  a  polygon  defined  In  A  with  "NPoints"  points,  line  style  "Line" } 
{  and  optional  Lines  from  the  axis  to  the  points  and  cross  hairs.  } 

procedure  Hatch(X_1,  Y_1,  X_2,  Y_2,  Delta  :  real); 

{  Hatch  a  bar  in  a  histogram  } 

procedure  DrawHistogram(A  :PlotArray;  NPoints  :  integer; 

Hatching  ;  boolean;  HatchStyle  :  integer); 

{  Draw  a  histogram  defined  in  A  with  "NPoints"  points  and  optional } 

{  hatching  with  a  given  hatch  style  } 

implementation 

procedure  Error{{ProcName  :  WrkString)}; 

{  Report  that  an  error  occurred  in  the  procedure  named  "ProcName" } 
begin 

DrawStringCERROR  in  '); 

DrawString(  ProcName) ; 

DrawStringC  Press  the  Button  to  exit.'); 
repeat  until  Button; 

Halt; 

end;  {  Error  } 
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procedure  SetPoregroundColor{(Color  :  integer}}; 

{  Set  the  foreground  drawing  color } 
begin 

case  Color  of 
0 ;  begin 

ForeColor(BlackColor) ; 

ForeColorGIb 0; 
end; 

1  ;  begin 

ForeColor(  WhiteColor)  ; 

ForeColorGIb  :=  1  ; 
end; 
end; 

end;  { SetForegrourtdColor } 

procedure  SetBackgroundColor{(Color  ;  integer)}; 

{  Set  the  background  color } 
begin 

case  Color  of 

0  ;  BackColor(BlackColor); 

1  :  BackColor( WhiteColor); 
end; 

end;  { SetBackgroundColor } 

procedure  DP{(X,  Y  ;  integer)}; 

{  Plot  a  pixel  at  position  (X,  Y) } 
begin 

MoveTo{X,  Y); 

UneTo(X,  Y); 
end;  {  DP } 

function  PD{(X,  Y  :  integer)  ;  boolean); 

{  Return  true  if  the  color  of  the  pixel  at  (X,  Y)  matches  ForeColorGIb  } 
var 

BlackPixel  ;  boolean; 

PixelColor  :  integer; 
begin 

BlackPixel  GetPixel(X,  Y); 
if  BlackPixel  then 
PixelColor 0 
else 

PixelColor  1; 

PD  >  PixelColor  ■  ForeColorGIb; 
end;  {  PD } 

procedure  DrawStraight{(X1 ,  X2,  Y  ;  integer)}; 

{  Draw  a  horizontal  line  from  X1  ,Y  to  X2,Y  } 
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begin 

MoveTo(X1,  Y); 

UneTo{X2.  Y); 
end;  {  DrawStraight  } 

procedure  lnvertWindow{(WinNum  :  integer)}; 

{  Invert  the  window  referenced  by  WinNum  } 
begin 

if  WinNum  in  [1..MaxWindowGlb]  then 

lnvertRect(Window[WinNum].W.Port.PortRect) 

else 

Error(’lnvert  Window'); 
end;  {  InvertWindow  } 

function  RealToStr{(R  :  real)  :  WrkString); 

{  Returns  the  string  representation  of  the  real  R. } 
var 

Int,  Frac  :  Longlnt; 

S1 ,  S2  :  Str255: 

Negitive  :  boolean; 
begin 

51  ;=  "; 

52  :=  "; 

if  R  <  0.0  then 
Negitive  :■  TRUE 
else 

Negitive  :*  FALSE; 

R  :=  ABS(R); 

Int  ;=  Trunc(R); 

Frac  :=  Round(100.0  *  (R  -  Int)); 
NumToString(lnt,  S1); 

NumToString(Frac,  S2); 
if  Length{S2)  =  1  then 
S2  ;=  S2  +  'O'; 

S2  :=  S1  +  +  S2; 

if  Negitive  then 
RealToStr '-'  +  S2 
else 

RealToStr ;«  "  +  S2; 
end;  {  RealToStr } 

function  lntToStr{(l  :  integer)  :  WrkString); 

{  Returns  the  string  representation  of  the  integer  t.  } 
var 

Form  :  DecForm; 

Str  :  DecStr; 
begin 
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Form.Styte  :«  FixedOecimal; 

Form.Digits  >  0; 

Num2Str(Form,  I,  Str); 

IntToStr  >  Str; 
end;  {  IntToStr } 

procedure  SetLineStyle{(Style  :  integer)}; 
var 

LineStyle  :  Pattern; 
begin 

case  Style  of 
0  ;  LineStyle  Black; 

1  :  LineStyle  >  White; 

2  :  LineStyle  >  Gray; 

3  :  LineStyle  :=  LtGray; 

4  :  LineStyle  :=  DkGray; 
otherwise 

Style  :=  0; 

LineStyle  :=  Black; 
end; 

LineStyleGIb  Style; 

PenPat(LlneStyle); 
end;  {  SetLineStyle  } 

function  GetLlneSty(e{  ;  Integer); 
begin 

GetUnestyle  LineStyleGIb; 
end;  { GetLineStyle } 

procedure  DefineWorld{(WorldNum  :  integer;  X_1,  Y_1,  X_2,  Y_2  :  real)}; 
begin 

if  ((X_1  <>  X_2)  and  (Y_1  <>  Y_2))  and 
(WorldNum  in  [1  ..MaxWorldsGIbj)  then 
with  World[WorldNum]  do 
begin 

XI  :=X_1; 

Y1  ;*  Y_1; 

X2  X_2; 

Y2  >  Y_2; 

if  WorldNum  >  MaxWortdGIb  then 
MaxWorkJGIb  WorldNum; 
end 

else  if  WorldNum  in  [1..MaxWorldsGlb]  then 
Error('Define World  #1') 
else 

Error('DeflneWorld  #2’); 
end;  {  DefineWorld  } 
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procedure  Select World{(WoridN urn  :  integer)}; 
begin 

if  (WorldNum  in  [1  ..MaxWorldGib])  then 
with  World[WortdNum]  do 
begin 

WoridNdxGib  >  WorldNum; 

XtWIdGIb:-  XI; 

YtWIdGIb  :=  Yt; 

X2WldGlb  >  X2; 

Y2WldGlb  :*  Y2: 
end 
else 

Error('Select  World'); 
end;  {  SelectWorld } 

procedure  DefineWindow{(WinNum,  XLo,  YLo,  XHi,  YHi,  WindowType  :  integer;  Visible 

boolean)}; 

var 

BoundsRect :  Rect; 

Title  :  Str255; 

RefCon  :  Longint; 

Visible  :  boolean; 

GoAwayRag ;  boolean; 

begin 

if  WinNum  in  [1  ..MaxWindowsGIb]  then 
begin 

if  WinNum  >  MaxWindowGIb  then 
MaxWindowGIb  :=  WinNum; 
with  BoundsRect  do 
begin 

Left XLo; 

Top  :=  YLo; 

Right XHi; 

Bottom  YHi; 
end; 

Window(WinNumJ.X1  :■  XLo; 

Window{WinNum].Y1  YLo; 

Window[WinNum].X2  :*  XHi; 

Window(WinNumi.Y2  :*  YHi; 

Title  "; 

GoAwayFlag  :■  TRUE; 

Visible FALSE; 

RefCon WinNum; 
if  WindowType  ■>  documentProc  then 
WindowType  16  *  WindowType  +  8;  {  Add  Zoom  box  } 
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Window[WinNum].P  :=  NewWindow(@Window[WinNum].W,  BoundsRect, 
Title,  Visible,  Wir^wType, 

POINTER(-I),  GoAwayFlag,  RefCon); 
end 
else 

Error('DefineWindow’): 
end;  {  DefineWindow } 

procedure  SetVisibility{(WinNum  :  integer;  Visible  :  boolean)}; 
begin 

ShowH ide (Window[Win N urn] . P ,  Visible) ; 
end;  {  Set  Visibility  } 

procedure  SelectWind{(WinNum  ;  integer;  Visible  :  boolean)}; 
begin 

if  (WinNum  in  [1  ..MaxWindowGIbj)  then 
begin 

SelectWindow(Window[WinNum].P); 
if  Visible  then 

ShowWindow(Window[WinNum).P); 
SetPort(@Window[WinNum].W.Port); 
with  Window[WinNum]  do 
begin 

WindowNdxGlb  :=  WinNum; 

XfRefGIb  :=  W.Port.PortRect.Left; 

Y1  RefGIb  :=  W.Port.PortRect.Top; 

X2RefGlb  ;=  W.Port.PortRect.Right; 

Y2RefGlb  :=  W.Port.PortRect.Bottom; 

BxGlb  (X2RefGlb  -  XfRefGIb  -  16)  /  (X2WldGlb  •  XIWIdGIb); 
ByGIb  :=  (Y2RefGlb  -  Y1  RefGIb  -  16)  /  (Y2WldGlb  -  YIWIdGIb); 
AxGlb  :»  XfRefGIb  -  XIWIdGIb  *  BxGlb; 

AyGIb  ;=  Y1  RefGIb  -  YIWIdGIb  *  ByGIb; 

AxisGIb  ;=  FALSE; 
end; 
end 
else 

Error('SelectWind'); 
end;  {  SelectWind } 

procedure  DefineHeader{( WinNum  :  integer;  Hdr  :  WrkString)}; 
begin 

WindowfWinNumj.Header  Hdr; 

SetWTitle( Window[WinNumJ ,P,  Hdr) ; 
end;  {  DefineHeader } 

procedure  RemoveHeader{(WinNum  :  integer)}; 
begin 
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DefineHeader(WinNum,  "); 
end;  { RemoveHeader } 

procedure  ReDefineWindow{(WinNum,  XI,  Y1,  X2,  Y2,  WindowType  ;  integer)}; 
begin 

if  (WinNum  in  [1  ..MaxWindowsGIb])  then 
b^in 

CioseWindow(Window[WinNum].P); 

WindowIWinNumJ.P  >  NiL; 

DefineWindow(WinNum,  X1,  Y1,  X2,  Y2,  WindowType); 
SelectWind(WinNum.  FALSE); 

DefineHeader( WinNum ,  Window[WinNum] .Header) ; 
end 
eise 

Error('ReDefineWindow'); 
end;  {  ReDefineWindow } 

function  WindowX{(X  :  real)  :  integer}; 
var 

Temp  :  real; 
begin 

Temp  :=  AxGlb  +  BxGlb  *  X; 
if  Temp  >  Maxint  then 
WindowX  Maxint 
else  if  Temp  <  -32767  then 
WindowX  :=  -32767 
else 

WindowX  :=  trunc(Temp); 
end;  { WindowX } 

function  WindowY{(Y  :  real)  :  integer}; 
var 

Temp  :  real; 
begin 

Temp  ;=  AyGIb  +  ByGIb  *  Y; 
if  Temp  >  Maxint  then 
WindowY  :=  Maxint 
else  if  Temp  <  -32767  then 
WindowY  ;=  -32767 
else 

WindowY  :=  trunc(Temp); 
end;  {  WindowY } 

procedure  ResetWorlds; 
var 

I  :  integer; 
begin 
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for  I  >  1  to  MaxWorldsGIb  do 
DefineWorld(l,  XMinGIb,  YMinGIb,  XMaxGlb,  YMaxGlb); 
SelectWorld(l); 
end:  {  ResetWorlds } 

procedure  InitGraphic; 


Index  ;  integer; 

BoundsRect :  Rect; 

Title  :  Str25S; 

Visible  :  boolean; 

RefCon  :  Longint; 

GoAwayRag :  boolean; 

WindowType  ;  integer; 
begin 

XMinGIb  screenBits.bounds.Left; 

YMinGIb  screenBits.bounds.Top; 

XMaxGlb  >  screenBits.bounds.Right; 

YMaxGlb  screenBits.bounds.Bottom; 
for  Index  1  to  MaxWindowsGIb  do 
begin 

Window[lndex].P  :=  NIL; 

Window[lndexl.H  ;=  NIL; 
end; 

ResetWorlds, 

MaxWorldGib  0; 

MaxWindowGIb  ;=  0; 

WindowNdxGlb  :=  0; 

WorldNdxGlb  :=  0; 

AspectFactor  :=  0.44; 

AspectGIb  :»  ABS(AspectFactor)  *  AspectFactor; 
AxisGIb  :=  false; 

HatchGIb  >  false; 

ClippingGIb  :=  true; 

SetLineStyle(O);  {  Solid  Black  lines  } 
end;  {  InitGraphic  } 

function  Clip{(var  XI,  Y1,  X2,  Y2  :  integer)  :  boolean}; 
var 

1x1,  lyl,  1x2,  Iy2,  Dummy,  XILoc,  X2Loc  :  integer; 
ClipLoc ;  boolean; 

Temp  :  real; 

function  lnside(X,  Xx1,  Xx2  ;  integer)  ;  integer; 
begin 

Inside  >  0; 
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if  X  <  Xx1  then 
Inside  >  -1 
else  if  X  >  Xx2  then 
Inside  1 ; 
end;  { Inside } 

begin  {  Clip } 

Clip  true; 

ClipLoc  true; 
if  ClippingGIb  then 
begin 

X1Loc:-  XIRefGIb; 

X2LOC  X2RefGlb; 

1x1  >  lnside(X1,  XILoc,  X2Loc); 
lyl  :=  lnside(Y1,  YlRefGIb,  Y2RefGlb); 

1x2  :=  lnside(X2,  XILoc.  X2Loc); 

Iy2  ;=  lnside(Y2,  YlRefGIb,  Y2RefGlb); 
if  (1x1  or  1x2  or  lyl  or  Iy2)  <>  0  then 
begin 

if  XI  <>  X2  then 
begin 

if  1x1  <>0  then 
begin 

if  Ixl  <  0  then 
Dummy  :=  XI  Loc 
else 

Dummy  X2Loc; 
if  Y2  <>  Y1  then 
begin 

Temp  :=  (Y2  -  Y1)  /  (X2  -  X1)  *  (Dummy  -  XI); 
if  Temp  >  Maxint  then 
Temp  :=  Maxint 
else  if  Temp  <  -32767  then 
Temp  ;=  -32767; 

Y1  :*  Y1  +  trunc(Temp); 
end; 

XI  :=  Dummy; 
end; 

if  (1x2  <>  0)  and  (XI  <>  X2)  then 
begin 

if  1x2  <  0  then 
Dummy :«  XILoc 
else 

Dummy ;»  X2Loc; 
if  Y2  <>  Y1  then 
begin 

Temp  >  (Y2  -  Y1)  /  (X2  -  XI)  *  (Dummy  -  XI); 
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if  Temp  >  Maxint  then 
Temp Maxint 
else  if  Temp  <  -32767  then 
Temp  >  -32767; 

Y2  :=  Y1  +  truncfTemp): 
end; 

X2  Dummy; 
end; 

lyl  >  lnside(Y1,  YIRefGIb,  Y2RefGlb); 

Iy2  lnside(Y2.  YIRefGIb,  Y2RefGlb); 
end; 

if  Y1  <>  Y2  then 
begin 

if  lyl  <>  0  then 
begin 

if  Iy1  <  0  then 
Dummy  :=  Y1  RefGIb 
else 

Dummy  :=  Y2RefGlb; 

If  XI  <>  X2  then 
begin 

Temp  ;=  (X2  -  X1)  /  (Y2  -  Y1)  *  (Dummy  -  Yl); 
if  Temp  >  Maxint  then 
Temp  ;=  Maxint 
else  if  Temp  <  -32767  then 
Temp  :a  -32767; 

X1  :=  X1  +  trunc(Temp); 
end; 

Yl  :=  Dummy: 
end; 

if  Iy2  <>  0  then 
begin 

if  Iy2  <  0  then 
Dummy  :=  Y1  RefGIb 
else 

Dummy  :=  Y2RefGlb; 
if  X1  <>  X2  then 
begin 

Temp  >  (X2  -  X1)  /  (Y2  -  Yl)  *  (Dummy  -  Yl); 
if  Temp  >  Maxint  then 
Temp  >  Maxint 
else  if  Temp  <  -32767  then 
Temp  >  -32767; 

X2  XI  +  trunc(Temp); 
end; 

Y2  :»  Dummy: 
end; 
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end; 

lyl  >  lnside(Y1.  YIRefQlb.  Y2RefGib); 

Iy2  >  lnside(Y2.  YIRefGIb,  Y2RefGlb): 
if  (Iy1  <>  0)  or  {Iy2  <>  0)  then 
ClipLoc  >  false; 
if  CiipLoc  then 
begin 

1x1  >  lnside(X1,  XlLoc,  X2Ljoc); 

1x2  lnside(X2,  XlLoc,  X2IjOC); 
if  (1x2  <>  0)  or  (1x1  <>  0)  then 
ClipLoc false; 
end; 

Clip  :«  ClipLoc; 
end; 
end; 

end;  {  Clip } 

procedure  DrawPoinl{(Xr,  Yr  :  real)}; 
var 

X,  Y  ;  integer; 
begin 

X  :=  WindowX(Xr); 

Y  :=  WindowY(Yr); 

DP(X.  Y); 
end;  {  DrawPoint  } 

function  PointDrawn{(Xr,  Yr  :  real)  :  boolean); 
begin 

PointDrawn  :=  PD(WindowX(Xr),  WindowY(Yr)); 
end;  {  PointDrawn  } 

procedure  DrawLine{(X1,  Y1,  X2,  Y2  :  real)}; 
begin 

MoveTo(WindowX(X  1 ),  WindowY( Y 1 )) ; 
LineTo(WindowX(X2),  WindowY(Y2)); 
end;  {  DrawLine  } 

procedure  DrawLineClipped{(X1 ,  Y1,  X2,  Y2  :  integer)}; 
begin 

if  Clip(X1,  Y1,  X2,  Y2)  then 
begin 

MoveTo(X1,  Y1); 

LineTo(X2,  Y2); 

{ end 
else 
begin 

Moveto(X1  -  3.Y1  -  3); 
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TextSize(9); 

DrawString(Ds); 

LabelSet  :*  False;} 
end; 

end;  {  DrawLineClipped  } 

procedure  DrawCrossDiag{(X,  Y,  Scale  :  integer)}; 
begin 

DrawLineClipped(X  -  Scale,  Y  +  Scale,  X  +  Scale  +  1 
DrawLineCiipped{X  -  Scale,  Y  -  Scale,  X  -i-  Scale  1 , 
end;  {  DrawCrossDiag } 

procedure  DrawWye{(X,  Y,  Scale  :  integer}}; 
begin 

DrawUneClipped(X  -  Scale,  Y  -  Scale,  X,  Y); 
DrawLineClipped(X  +  Scale,  Y  -  Scale,  X,  Y); 
DrawUneClipped(X,  Y,  X,  Y  +  Scale); 
end;  {  DrawWye  } 

procedure  DrawDiamond{(X,  Y,  Scale  :  integer)}; 
begin 

DrawLineClipped{X  -  Scale,  Y,  X,  Y  -  Scale  -  1); 
DrawLineClipped(X,  Y  -  Scale  +  1,  X  +  Scale,  Y  +  1); 
DrawlJneClipped(X  +  Scale,  Y  +  1 ,  X,  Y  +  Scale); 
OrawLineClippedix,  Y  +  Scale,  X  -  Scale,  Y); 
end;  {  OrawOiamond } 

procedure  DrawCircleDirect{(Xr,  Yr,  R  :  integer)}; 
type 

Circ  =  array[1..14]  of  integer; 
var 

Xkl ,  Xk2,  Yk1 ,  Yk2,  Xp1 ,  Yp1 ,  Xp2,  Yp2  :  integer; 
Xfact,  Yfact :  real; 

I  :  integer; 

X  :  Circ; 

procedure  InitX; 
begin 
X(11  :=  0; 

X(21  121; 

X(3)  :=  239; 

X[41  355; 

X15}  465; 

X[6]  568; 

X[7]  >  663; 


,  Y  -  Scale  -  1); 
Y  +  Scale  +  1); 
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X(8]  >  749: 

X[9]  >  823; 

X[10]  :>  885: 

X[111  935: 

X(12]  971: 

X[13]  >  993; 

XII 4]  1000: 

end;  {  InitX  } 

begin  {  DrawCircleDirect  } 

InitX; 

Xfact  :=  abs(R); 

Yfact  :=  Xfact  *  AspectGIb; 
if  Xfact  >  0.0  then 
begin 

Xkl  :=  trunc(X[11  *  Xfact  +  0.5): 

Ykl  :=  trunc(X(14]  *  Yfact  +  0.5); 
for  I  :=  2  to  14  do 
begin 

Xk2  :=  trunc(X[l]  *  Xfact  +  0.5); 

Yk2  :=  trunc(X[14  -  I  +  1]  *  Yfact  +  0.5): 

Xpl  :=  Xr  -  Xkl; 

Ypl  :=  Yr  +  Ykl ; 

Xp2  :=  Xr  -  Xk2: 

Yp2  :«  Yr  +  Yk2; 

DrawLine(Xp1 ,  Ypl,  Xp2,  Yp2): 

Xpl  ;=  Xr  +  Xkl ; 

Xp2  :=  Xr  +  Xk2; 

DrawLine(Xp1 ,  Ypl,  Xp2,  Yp2): 

Ypl  :=  Yr  -  Ykl; 

Yp2  ;=  Yr  -  Yk2: 

DrawLine(Xp1,  Ypl  +1,  Xp2,  Yp2  +  1); 

Xpl  ;-Xr-Xk1; 

Xp2  Xr  -  Xk2: 

DrawLine(Xp1 ,  Ypl  +  1,  Xp2,  Yp2  +  1); 

Xkl  ;>  Xk2; 

Ykl  ;=  Yk2: 
end; 
end 
else 

DP{Xr,  Yr); 

end;  {  DrawCircleDirect  } 

procedure  DrawCircle{(X_R,  Y_R,  Xradius  :  real)}; 
begin 

DrawCircleDirect(WindowX(X_R),  WindowY{Y_R),  trunc(Xradius)): 
end;  {  DrawCircle  } 
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procedure  DrawCross{(X1 .  Y1.  Scale  :  integer)}; 
begin 

DrawLineCiipped(X1  -  Scaie,  Y1,  X1  -i-  Scale  +  2,  Y1); 
DrawLineClip^(X1 ,  Y1  -  Scale,  XI.  Y1  Scale  +  1); 
end;  {  DrawCross } 

procedure  DrawStar{(X.  Y,  Scale  :  integer)}; 
begin 

DrawLineClipped(X  -  Scale,  Y  -t-  Scale,  X  +  Scale  +  1,  Y  -  Scale  -  1); 
DrawUneClipped(X  -  Scale,  Y  -  Scale.  X  Scale  1,  Y  +  Scale  +  1); 
DrawLineClipped(X  -  Scale  -  2,  Y.  X  +  Scale  +  4.  Y); 
end;  {  DrawStar } 

procedure  DrawSquareC{(X1 ,  Y1,  X2.  Y2  :  integer;  Fill  :  boolean)}; 
var 

I  :  integer; 

procedure  DSC(X1,  X2,  Y  ;  integer); 
begin 

DrawStraight(X1 ,  X2,  Y); 
end;  { DSC } 

begin  {  DrawSquareC } 

If  not  Fill  then 
begin 

DrawLineClipped(X1 ,  Y1,  X2,  Y1); 

DrawLineClipped(X2,  Y1,  X2,  Y2); 

DrawLineClipped(X1,  Y2.  X2,  Y2); 

DrawLineClipped(Xl ,  Y2,  X1,  Y1); 
end 
else 

for  I  :=  Y2  to  Y1  do 
DSC{X1,  X2, 1); 
end;  { DrawSquareC } 

procedure  DrawSquare{(X1 ,  Y1,  X2.  Y2  :  real;  Fill  :  boolean)}; 
var 

I.  XILoc,  YIIjOC,  X2LOC,  Y2Loc  :  integer; 

DirectModeLoc :  boolean; 

procedure  DS(X1,  X2,  Y  :  integer); 
begin 

if  LineStyleGIb  -  0  then 
DrawStraight(X1,  X2,  Y) 
else 

DrawLine(X1.  Y,  X2,  Y); 


end;  { DS } 

procedure  DSC(X1,  X2,  Y  :  integer): 
begin 

DS(X1.  X2.  Y): 
end:  { DSC } 

procedure  DrawSqr(X1,  Y1,  X2,  Y2  :  integer;  Fill  :  boolean): 
var 

I  :  integer; 
begin 

if  not  Fill  then 
begin 

DS(X1.X2,  Y1): 

DrawLine(X2.  Y1,  X2,  Y2): 

DS(X1.  X2.  Y2): 

DrawLine(X1,  Y2,  XI.  Y1); 
end 
else 

for  I  ;=  Y1  to  Y2  do 
DS(X1,  X2.  I); 
end;  {  DrawSqr } 

begin  {  DrawSquare  } 

XILoc  :»=  WindowX(XI); 

YILOC  :*  WindowY(YI): 

X2LOC  :=  WindowX{X2); 

Y2Loc  :=  WindowY(Y2); 
if  not  Fill  then 
begin 

DSC(X1  Loc,  X21jOC,  Y1  Loc); 

DrawLineClipped(X2Loc,  YIL0C,  X2Loc,  Y2Loc): 
DSC(X1Loc,  X2L0C,  Y2Loc): 

DrawLineClipped(X1Loc,  Y2Loc,  XILoc,  Y1Loc); 
end 
else 

for  I  >  Y1  Loc  to  Y2Loc  do 
DSC{X1Loc,  X2LOC.  I): 
end;  {  DrawSquare } 

procedure  DrawAscii{(X,  Y  :  integer;  Size,  CharByte  ;  byte)); 
begin 

MoveTo(X.  Y): 

TextSize(Size  *  12); 

DrawChar(Chr(CharByte)): 
end;  {  DrawAscii  } 
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procedure  DrawText{(X,  Y.  Scale  :  integer;  Txt  :  WrkString)}; 
var 

Index  :  integer; 

EscStr :  boolean; 

StringLen  :  integer; 

AsciiValue  :  integer; 

SymbolScaie  :  integer; 

SymboiCode :  integer; 
begin 
Index  >  1 ; 

EscSir :«  FALSE; 

StringLen  Length(Txt); 

while  (Index  <«  StringLen)  and  (not  EscStr)  do 

begin 

if  Txt[lndex]  =  #27  then 
EscStr  :=  TRUE; 

Index  ;=  Index  +  1 ; 
end; 

if  not  EscStr  then 
begin 

MoveTo(X.  Y); 

TextSize(Scale  *  12); 

DrawString(Txt): 
end 
else 
begin 
Index  ;=  1 ; 

while  Index  <=  StringLen  do 
begin 

AsciiValue  :=  Ord(TxtIlndex]); 
if  AsciiValue  -  27  then 
begin 

SymbolScaie  ;=  Scale; 

Index  Index  +  1 ; 
if  Index  <«  StringLen  then 
begin 

SymboiCode  Ord(T xt[lndex})  -  48; 

If  (Index  +  2  <=  StringLen)  and  (Ord(TxtIlndex  +  1J)  =  64)  then 
begin 

SymboiCode  Ord(Txt[lndexl)  -  48; 

Index  Index  +  2; 
end; 

case  SymboiCode  of 

1  :  DrawCross(X  +  SymbolScaie,  Y  +  Scale,  SymbolScaie); 

2  :  DrawCrossDiag(X  +  SymbolScaie,  Y  +  Scale,  Symbol^le); 
3,4  :  DrawSquareC(X,  Y  +  (SymbolScaie  shl  1)  -  1, 

X  +  (SymbolScaie  shl  1),  Y  - 1 ,  (SymboiCode  -  4)); 
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5  :  begin 

DrawDiamond{X  +  trunc{1.5  *  SymbolScale), 

Y  +  SymbolScale  - 1,  SymbolScale  +  1): 

X  X  SymbolScale; 
end; 

6  :  DrawWye(X  +  SymbolScale.  Y  +  SymbolScale  - 1 ,  SymbolScale): 

7  :  begin 

DrawStar(X  +  SymbolScale  shl  1.  Y  +  SymbolScale  - 1,  SymbolScale): 
X  X  +  SymbolScale  shl  1; 
end; 

8  :  DrawCircleDirect(X  +  SymbolScale.  Y  +  (SymbolScale  shr  1), 
SymbolS^le  + 1); 

end; 

X  X  +  3  *  SymbolScale: 

SymbolScale  :=  Scale; 
end; 
end 
else 

DrawAscii(X,  Y.  Scale,  AsciiValue); 

Index  :=  Index  +  1 ; 
end; 
end; 

end;  {  DrawText } 

procedure  DrawTextW{(X,  Y  :  real;  Scale  :  integer;  Txt  :  WrkString)}: 
begin 

DrawTexl{WindowX(X),  WindowY(Y).  Scale.  Txt); 
end;  {  DrawTextW  } 

procedure  TexlStyle{{Face  :  Style)}; 

{  Face  =  (bold,  italic,  underline,  outline,  shadow,  condense,  extend) } 
begin 

TextFace(Face); 
end:  (  TextStyle  } 

procedure  HardCopy{(TopWin  :  boolean)}; 
begin 

PrDrvrOpen; 
if  TopWin  then 
{  Print  the  top  folder. } 

PrCtlCall(iPrEvtCtl,  LPrEvtTop,  0,  LScreenBits) 
else 

{  Print  the  whole  screen. } 

PrCtlCall(iPrEvtCtl,  LPrEvtAil,  0,  LScreenBits): 

PrDrvrClose; 
end;  { HardCopy } 
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procedure  OpenPic{(WinNum  :  integer;  ShowPic  :  boolean)}; 
begin 

if  WindowIWinNumJ.H  <>  NIL  then 
begin 

KillPicture(WindowIWinNum].H); 

Window[WinNum].H  >  NIL; 
end; 

RectRgn(Window[WinNum].W.Port.cfipRgn,  ScreenBits.bounds); 
Window[WinNum].H  :=  OpenPicture(Wlndow(WinNum].W.Port.PortRect); 
if  ShowPic  then 
ShowPen 
end;  { OpenPic } 

procedure  DrawPic{(WinNum  :  integer)}; 
var 

PictRect :  Rect; 
begin 

PictRect  :=  Window[WinNum}.W.Port.PortRect; 
with  PictRect  do 
begin 

if  {(Bottom  -  Top)  <  200)  OR  ((Right  -  Left)  <  200)  then 
begin 

Right  :=  Right  -  16;  {  so  we  don't  overwrite  the  grow  region  } 

Bottom  ;=  Bottom  - 16:{  on  a  small  window.  } 

end; 
end; 

DrawPicture(Window(WinNum).H,  PictRect); 
end;  {  DrawPic  } 

procedure  ErasePic{(WinNum  :  integer)}; 
begin 

if  Window[WinNum].H  <>  NIL  then 
begin 

KiiiPicture(Window[WinNum].H); 

Window[WinNum].H  ;=  NIL; 
end; 

end;  {  ErasePic } 

procedure  ClearWindow{(WinNum  :  integer)}; 
begin 

EraseRect{Window  (WinNum]  .W.Port.PortRect) ; 
end;  {  ClearWindow } 

function  WhereX{  :  integer}; 
var 

Pt  :  Point: 
begin 
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GetPen(R); 

WhereX  >  Pt.H; 
end;  { WhereX } 

function  WhereY{  :  integer}; 
var 

R  :  Point; 
begin 

GetPen(R); 

WhereY  :=  R.V; 
end;  { WhereY } 

procedure  SetWindow{(X1 ,  Y1,  X2,  Y2  :  integer)}; 
begin 

XtRefGIb  :=  X1; 

YIRefGIb  ;=  Y1; 

X2RefGlb  :=  X2; 

Y2RefGlb  :=  Y2; 

BxGlb  :=  (X2  •  X1)  /  (X2WldGlb  -  XtWIdGib); 

ByGIb  :=  (Y2  -  Y1)  /  (Y2WldGlb  -  YIWIdGIb); 

AxGlb  :=  XI  -  XlWIdGIb  *  BxGlb; 

AyGIb  Y1  -  YIWIdGIb  *  ByGIb; 

AxisGIb  :=  FALSE; 
end;  { SetWindow } 

procedure  FindWorld{(l  ;  integer;  A  :  PlotArray;  NPoints  :  integer)}; 
var 

J  :  integer; 

Xmax,  Ymax,  Xmin,  Ymin,  Xmid,  Ymid,  Xdiff,  Ydiff :  real; 
begin 

NPoints  :=  abs(NPoints); 
if  NPoints  >  2  then 
if  I  in  [1..MaxWorldsGlb]  then 
begin 

Xmax  :=  A[1,  1]; 

Ymax  Ajl,  2]; 

Xmin  Xmax; 

Ymin  Ymax; 

for  J  :«  2  to  NPoints  do 

begin 

if  A[J,  1]  >  Xmax  then 
Xmax  >  A[J,  1] 
else 

if  A[J,  1]  <  Xmin  then 
Xmin  :=  A[J,  1]; 
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if  AIJ,  2]  >  Ymax  then 
Ymax  >  A[J,  2] 
else 

if  A[J,  2]  <  Ymin  then 
Ymin  :=  A[J,  2]; 

end; 

Xmin  >  Round(Xmin); 

Ymin  >  Round(Ymin)  -  0.5; 

Xmax  >  Round(Xmax); 

Ymax  >  Round(Ymax)  +  0.5; 

DefineWorld(l,  Xmin,  Ymin,  Xmax,  Ymax); 

SelectWorld(l); 

end 

else 

Error('FindWorld  #1') 
else 

Error('FindWor!d  #  2'); 
end;  {  FindWorld  } 

procedure  FindWorld1{(l  :  integer;  A  :  PlotArray;  NPoints  :  integer)}; 
var 

J  :  integer; 

Xmax,  Ymax,  Xmin,  Ymin,  Xmid,  Ymid,  Xdiff,  Ydiff :  real; 
begin 

Xmin  :=  XMn; 

Ymin  :=  YMn{Round(YMn)}  ; 

Xmax  :=  XMx{Round(XMx)}; 

Ymax  :=  YMx{Round(YMx)}  ; 

DefineWorld{l,  Xmin,  Ymin,  Xmax,  Ymax); 

SelectWorld(l); 
end;  {  FindWorld  } 

procedure  DrawAxis{(Footer1 ,  Footer2  :  WrkString;  Arrows  ;  boolean)}; 
var 

LineStyleljOC,  XkO,  YkO,  Xk1,  Ykl,  Xk2,  Yk2, 

MaxExponentX,  MaxExponentY,  TIckPoInt :  Integer; 

TickSmail,  TickLarge,  Max,  Min,  Tick,  Offset,  Diff  :  real; 

function  Log(X  :  real)  :  real; 

{  Base  10  logarithm  of  X. } 
begin 

Log  >  Ln(X)  /  Ln(IO.O); 
end;  {  Log } 

function  ALog(X  ;  real)  ;  real; 
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{Ten  raised  to  the  X  power. } 
begin 

ALog  Exp(X  *  Ln(IO.O)): 
end;  { Alog } 

function  Frac{R  :  real)  :  real; 

{  Return  the  fractional  part  of  the  real  number  R. } 
begin 

Frac  :«  R  -  Int(R): 
end:  {  Frac } 

procedure  Ticks(NTicks  :  integer;  Max,  Min  :  real; 
var  TickSmall,  TickLarge  :  real); 

{  NTicks  ;  The  approximate  number  of  tick  marks  in  the  interval.  ) 
{  Max,  Min  :  World  coordinates  of  the  axis  extremes.  } 

{  TickSmall,  TickLarge  :Tick  mark  intervals,  in  world  coordinates  .) 
var 

TickLog  :  array [1.. 4]  of  real; 

I,  J,  ChA  :  integer: 

Delta.  XTicks,  LogTicks,  Mant,  MinDiff,  Diff  :  real; 
begin 

TickLogff]  :=  0.0; 

TickLog[2]  :=  Log(2.0): 

TickLog[3]  :=  Log(5.0): 

TickLogI4]  1.0; 

XTicks  ;=  NTicks; 

Delta  Max  -  Min; 

XTicks  :=  Delta  /  XTicks; 

LogTicks  :*  Log(XTicks): 

ChA  :a  trunc(LogTicks): 
if  LogTicks  <  0.0  then 
ChA  :=  ChA  - 1  ; 

MinDiff  :=  1.0; 

Mant  :=  LogTicks  -  ChA;  {  Fractional  part  of  logarithm  } 

for  I  :=  1  to  4  do 

begin 

Diff  :=  Abs(Mant  -  TickLog(l)); 
if  (Diff  <  MinDiff)  then 
begin 

MinDiff  Diff; 

J  :=l: 
end; 
end; 

LogTicks  >  ChA  +  TickLog(J];  {  Logarithm  of  tick  mark  } 
TickUrge  :=  A LOg( LogTicks):  {  The  tick  mark  } 

{  Find  the  small  tick  marks,  that  are  two  tick  scales  smaller } 
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J  :=  J  -  2; 

LogTicks  :=  ChA  +  Tickljog[Jl: 
end 
else 
begin 
J  :=  J  +  1; 

LogTicks  :=  ChA  +  TickLog[Jl  - 1 : 
end; 

TickSmall  >  ALog( LogTicks); 
end;  {  Ticks  } 

function  StringNumber(X1  :  real;  MaxExponent  :  integer)  :  WrkString; 
begin 

StringNumber  :=  RealToStr{X1  *  Exp(-MaxExponent  *  Ln(10.0))); 
end;  {  StringNumber  } 

function  GetExponent(X1  :  reai)  :  integer; 
begin 

GetExponent  ;=  0; 
if  X1  <>  0.0  then 
if  ABS(X1)  >=  1.0  then 

GetExponent  :=  trunc(Ln(ABS(X1))  /  Ln(IO.O)) 
gIsg 

GetExponent  :=  -trunc(ABS(Ln{ABS(X1)))  /  Ln(IO.O)  +  1.0); 
end;  { GetExponent } 

procedure  DrawNum{X1,  Y1,  MaxExponent  :  integer;  Number  :  real); 
var 

StrNumber  :  WrkString; 
begin 

TextSize(9); 

StrNumber  :=  StringNumber(Number.  MaxExponent); 

Y1  :=  Y1  -  3; 

MoveTo(Xl,  Y1); 

DrawString(StrNumber) ; 

TextSize{12); 
end;  {  DrawNum  } 

procedure  DrawExponent(X1 ,  Y1,  MaxExponent  :  integer); 
var 

NumStr  :  WrkString; 
begin 

MoveTo(X1,  Y1); 

TextSize(9); 

DrawChar('x'); 


DrawCharC  '); 

DrawChar(’l'): 

DrawChar('O'): 

XI  WhereX  +  1; 

Y1  WhereY  -  3; 

MoveTo(X1,  Y1); 

NumStr  >  IntToStr(MaxExponent); 

TextSize(7): 

DrawString(NumStr); 

TextSi2e(12): 
end;  {  DrawExponent } 

begin  {  DrawAxis  } 

UneStyleLoc  :=  LinestyleGIb; 

SetLineStyle(O);  { Black } 

XkO  :=  XIRefGIb  +  X1  Offset: 

YkO  :=  Y2RefGlb  -  Y20ffset: 

Xk1  :=  XkO; 

Yk1  :=  YIRefGIb  +  Y1  Offset; 

Xk2  ;=  X2RefGlb  -  X20ffset: 

Yk2  :=  YkO; 

MoveTo(XK0,  YKO);  {  Draw  the  Y  axis  with  optional  Arrows  } 
LineTo(XK1.  YK1); 
if  Arrows  then 
begin 

MoveTo(XkO,  Yk1): 

LineTo(XkO  -  4,  Yk1  +  4); 

MoveTo(XkO,  Ykt); 

UneTo(XkO  +  4,  Yk1  +  4); 

DP(XkO,  Yk1  -  1): 
end; 

MoveTo(XkO,  YkO);  {  Draw  the  X  axis  with  optional  Arrows  } 
LineTo(Xk2  +  1,  Yk2): 
if  Arrows  then 
begin 

MoveTo(Xk2,  Yk2): 

UneTo(Xk2  -  4,  Yk2  -  4); 

MoveTo(Xk2,  Yk2); 

LineTo(Xk2  -  4.  Yk2  +  4); 
end; 

if  Footerl  <>  ”  then  {  Draw  the  1  st  footer  below  the  X  axis  } 
begin 

MoveTo(XkO,  YkO  +  45); 
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TextSi2e(9); 

DrawString(Footer1 ) ; 
end; 

if  Footer2  <>  "  then  {  Draw  the  2nd  footer  befow  the  X  axis  } 
begin 

MoveTo(XkO,  YkO  +  65); 

TextSi2e(9); 

DrawString(Footer2) ; 
end; 

if  (ABS(YkO  -  Ykl)  >=  35)  and  (ABS{Xk2  -  Xkl)  >-  150)  then 
begin 

if  ABS(Y2WldGlb)  >  ABS(YIWIdGlb)  then 
MaxExponentY  :=  GetExponent(Y2WidGlb) 
else 

MaxExponentY  :=  GetExponent(YIWIdGlb); 

if  MaxExponentY  <?.  0  then  {  Draw  the  power  of  ten  on  top  of  Y  axis  } 
DrawExponent{Xk1  -  30,  Ykl  +  2,  MaxExponentY); 

TickPoint  :=  YkO; 
if  Y1  WIdGIb  >  Y2WldGlb  then 
begin 

Max  Y1  WIdGIb; 

Min  :=  Y2WldG!b; 
end 
else 
begin 

Max  :=  Y2WldGlb; 

Min  :=  Y1  WIdGIb; 
end; 

{  Using  the  Max  and  Min  values,  this  procedure  call  calculates  } 

{ large  and  small  Tick  Marks  in  world  coordinates.  } 

Tlcks(5,  Max,  Min,  TickSmall,  TickLarge); 

Offset  ;=  Min  /  TIcklarge; 

Offset  Offset  -  Frac{Offset); 

Tick  Offset  *  TickLarge; 

If  Tick  <  Min  then 
Tick  >  Tick  +  TickLarge; 

{  Tick  is  the  world  coordinate  at  which  the  tick  mark  is  to  be  drawn  } 
Diff  >  Max  •  Min; 

{  Plot  large  tick  marks  and  Numeric  labels  } 
while  Tick  <-:  Max  do 
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begin 

TickPoint  ;=  YkO  -  Trunc((YkO  -  Yk1)  *  {Tick  -  Min)  /  Diff); 
MoveTo(XkO,  TickPoint); 

LineTo(XkO  -  4,  TickPoint); 

DrawNum(X1  Offset  -  30,  TickPoint  +  7,  MaxExponentY,  Tick); 
Tick  ;=  Tick  +  TickLarge; 
end; 

{  The  same  repeated  for  the  small  tick  marks, } 

{  only  without  axis  numbering.  } 

Offset  >  Min  /  TickSmall; 

Offset  >  Offset  -  Frac(Offset); 

Tick  :=  Offset  *  TickSmall; 
if  Tick  <  Min  then 
Tick  :=  Tick  +  TickSmall; 
while  (Tick  +  0.01)  <  Max  do 
begin 

TickPoint  :=  YkO  -  Trunc((YkO  -  Ykl)  *  (Tick  -  Min)  /  Diff); 
MoveTo(XkO,  TickPoint); 

LineTo(XkO  -  2,  TickPoint); 

Tick  :=  Tick  +  TickSmall; 
end; 

if  ABS(X2WldGlb)  >  ABS(XIWIdGlb)  then 
MaxExponentX :»  GetExponent{X2WldGlb) 
else 

MaxExponentX  :>=  GetExponent(XIWIdGlb); 
if  MaxExponentX  <>  0  then  {  Draw  power  of  ten  label  on  X  axis  } 
DrawExponent(Xk2  -  25,  YkO  +  28,  MaxExponentX); 

{  This  is  the  same  as  for  the  Y  axis,  but  the  window } 

{  and  world  are  appropriate  for  the  X  axis.  } 

TickPoint  :=  XkO; 
if  XI  WIdGIb  >  X2WldGlb  (hen 
begin 

Max  :=  XI  WIdGIb; 

Min  :=  X2WldGlb; 
end 
else 
begin 

Max  X2WldGlb; 

Min  >  XI  WIdGIb; 
end; 

Ticks(5,  Max,  Min,  TickSmall,  TickLarge); 

Offset  Min  /  TickLarge; 

Offset  >  Offset  -  Frac(Offset); 

Tick  :«  Offset  *  TickLarge; 
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if  Tick  <  Min  then 
Tick  >  Tick  +  TickLarge; 

Diff  >  Max  -  Min; 
while  Tick  <«  Max  do 
begin 

TickPoint  XkO  +  Trunc((Xk2  -  XkO)  *  (Tick  -  Min)  /  Diff); 
MoveTo(TickPoint,  YkO); 

LjneTo(TickPoint,  YkO  +  4); 

DrawNum(TickPoint  -  14,  YkO  +  20,  MaxExponentX,  Tick); 
Tick  :=  Tick  +  TickLarge; 
end: 

Offset  >  Min  /  TickSmall; 

Offset  :=  Offset  -  Frac(Offset): 

Tick  :*  Offset  *  TickSmall; 
if  Tick  <  Min  then 
Tick  ;=  Tick  +  TickSmall; 
while  (Tick  +  0.01)  <  Max  do 
begin 

TickPoint  .=  XkO  +  Trunc((Xk2  -  XkO)  *  (Tick  -  Min)  /  Diff); 
MoveTo(TickPoint,  YkO); 

UneTo(TickPoint.  YkO  +  2); 

Tick  ;=  Tick  +  TickSmall; 
end; 
end; 

Set  LineSty  ie  ( Lin  eStyle  Loc) : 

AxisGIb  :=  TRUE; 
end;  {  DrawAxis  } 

procedure  ResetAxis; 
begin 

AxisGIb  true; 
end;  {  ResetAxis } 

procedure  DrawPolygon{(A  :  PlotArray;  First,  NPoints,  Line,  Scale, 
Lines  :  integer;  CrossHairs  :  boolean)}; 
var 

I,  XI,  X2,  Y1,  Y2,  XOffset,  YOffset, 

XI  RefLoc,  YIRefLoc,  X2RefLoc,  Y2RefLoc, 

DeitaY,  XOsI ,  XOs2,  YOsI ,  YOs2  :  integer; 

AutoClip,  DirectModeLoc,  PlotUne,  PlotSymbol,  Flipped  :  boolean; 
XI  Loc,  Y1  Loc,  X2LOC,  Y2Loc  :  integer; 

Temp  :  real; 

LineStyleLoc2  :  integer; 

DrawR  :  Boolean; 

procedure  DrawPointClipped(X,  Y  :  integer); 
begin 
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if  (XI  >  XIRefGIb)  and  (X2  <  X2RefGlb)  then 
if  (Y1  >  YIRefGIb)  and  {Y2  <  Y2RefGlb)  then 
DP(X,  Y): 

end;  {  DrawPointClipped  } 

procedure  Drawltem{X,  Y  ;  integer); 
var 

LineStyleLoc  :  integer; 
begin 

LineStyleLoc  :=  LineStyleGIb; 

SetLineStyie(O);  {  Biack  } 
case  Line  of 

2  :  DrawCrossDiag(X,  Y,  Scale); 

3,  4  :  DrawSquareC(X  -  Scale,  Y  +  Scale.  X  +  Scale.  Y  -  Scale,  (Une  =  4)) 

5  ;  DrawDiamond(X,  Y,  Scale  +  1); 

6  :  DrawWye(X,  Y,  Scale  +1); 

1  :  DrawCross(X,  Y,  Scale); 

8  :  DrawCircleDirect(X.  Y,  Scale  +  1); 

9  :  begin 

PlotLine  :=  false; 
if  AutoClip  then 
DrawPointClipped(X,  Y) 
else 

DP(X,  Y); 
end; 

7  ;  DrawPoint(X,  Y){,  Scale)}; 

end; 

SetLineStyle(  LineStyleLoc); 
end;  {  Drawitem  } 

begin  {  DrawPolygon  } 
if  AxisGIb  then 
Flipped  :=  FALSE 
else 
begin 

Flipped  :=  TRUE; 

Temp  >  World[WorldNdxGlbl.Y1 ; 

World{WorldNdxGlbl.Y1  >  World[WorldNdxGlbI.Y2; 
World[WorldNdxGlb].Y2  ;=  Temp; 

Select  World(WorldNdxGlb) ; 

SelectWind(WindowNdxGlb,  TRUE); 
end; 

if  abs(NPoints  -  First)  >=  2  then 
begin 

AutoClip  :«  (N  Points  <  0): 

NPoints  >  abs(NPoints); 

XOsI  :=  1 ; 
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X0s2  1  : 

YOsI  >  6: 

YOs2  6; 
if  AxisGIb  then 
begin 

XOs1  XI  Offset; 

XOs2 X20ffset; 

YOs1  :-Y1  Offset: 

YOs2  >  Y20ffset: 

if  ((X2RefGlb  -  X0s2  -  XIRefGIb  +  XOsl)  >  (XOsI  +  XOs2))  and 
((Y2RefGlb  -  YOs2  -  YIRefGIb  +  YOs1)  >  (YOsl  +  YOs2))  then 
begin 

XtRefLoc XIRefGIb; 

X1  X1RefGlb  +  XOs1: 

YtRefLoc  :=  YIRefGIb; 

Y1  :=  YIRefGIb +  YOs1: 

X2RefLoc  >  X2RefGlb; 

X2  ;=  X2RefGlb  -  XOs2: 

Y2RefLjOC  :=  Y2RefGlb: 

Y2  :=  Y2RefGlb  -  YOs2: 

SetWindow(X1,  Y1,  X2.  Y2): 

AxisGIb  :=  TRUE; 
end; 
end; 

PlotLine  ;=  (Line  >=  0); 

PlotSymbol  :=  (Line  <>  0); 

Line  abs(Line); 

Scale  ;=  abs(Scale); 
if  Lines  <  0  then 

DeltaY  ;=  Trunc(1.0  /  (abs(YIWIdGlb)  +  abs(Y2WldGlb))  * 

abs(YIWIdGlb)  *  abs(Y2RefGlb  -  YIRefGIb))  +  1 
else 

DeltaY  :=  0; 
if  (NPoints  <  2)  then 
Error('DrawPolygon  #1') 
else 
begin 

if  CrossHairs  then 
begin 

LineStyleljoc2  LineStyleGIb; 

SetLineStyle(3);  {  Light  Gray  } 

MoveTo(XlRefGlb,  YIRefGIb  +  Y2RefGlb  -  WindowY(O.O)): 
LineTo(X2RefGlb,  YIRefGIb  +  Y2RefGlb  -  WindowY(O.O)); 
MoveTo(WindowX(0.0),  YIRefGIb); 

UneTo(WindowX(0.0),  Y2RefGlb); 
SetLineStyle(UneStyleLoc2); 
end; 
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XI  >  WindowX{A[First,  1]); 

Y1  YIRefGIb  +  Y2RefGlb  -  WindowY(A[Flrst,  2])  ; 

Drawltem(X1,  Yl): 
if  Abs(Lines)  -  1  then 
if  AutoClip  then 

DrawUneClipped(X1,  Y2RefGlb  -  DeltaY,  XI,  Y1) 
eise 
begin 

MoveTo(X1.  Y2RefGlb  -  DeitaY); 

LineTo(X1,  Yl): 
end; 

DrawR  :■  True: 
for  i>  First  1  to  NPoints  do 
begin 

X2  :=  WindowX(A[l,  1]); 

Y2  :=  Y2RefGlb  +  YIRefGIb  -  WindowY(A[l,  2)): 

Drawltem(X2,  Y2): 
if  StepA  then  begin 
if  not  AMarkStatus  then 
begin 

if  (A[l.  2]  >  0)  and  DrawPt  then  {Clip(X1.Y1,X2,Y2)} 
begin 

if  ARJustification  then 
begin 

MoveTo(X2+15,  Y2+5): 

TextSize(9): 

DrawString(Os); 

end 

else 

begin 

MoveTo{X2-65,  Y2): 

TextSize(9): 

DrawString(Ds); 

end: 

DrawPt  :=  False; 
end; 
end 
else 
begin 

if  (I  >  (NPoints  -  InitDegree))  and  (All,  2]  >  0)  and  DrawPt  then 
begin 

if  ARJustification  then 
begin 

MoveTo(X2+15,  Y2+5): 

TextSize(9): 

DrawString(Ds): 

end 
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else 

begin 

MoveTo(X2-65,  Y2): 

TextSize(9); 

DrawString(Os); 

end; 

DrawR  :=  False: 
end; 
end; 

end;{StepA} 
if  not  StepA  then  begin 
if  not  BMarkStatus  then 
begin 

if  (A[l,  2]  >  0)  and  DrawR  then  {Clip{X1  .Y1  .X2.Y2)} 
begin 

if  BRJustification  then 
begin 

MoveTo(X2+15.  Y2+5): 

TextSi2e(9); 

DrawString(Os); 

end 

eise 

begin 

MoveTo(X2-65,  Y2): 

TextSize(9); 

DrawString(Ds); 

end: 

DrawPt  :=  False; 
end; 
end 
else 
begin 

if  (I  >  (NPoints  -  InitDegree))  and  (Ali,  2J  >  0)  and  DrawPt  then 
begin 

if  BRJustification  then 
begin 

MoveTo(X2+15.  Y2+5); 

TextSize(9); 

DrawString(Ds); 

end 

else 

begin 

MoveTo(X2-65.  Y2): 

TextSize(9); 

DrawString(Ds); 

end; 

DrawR  :=  False; 
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end; 

end; 

end;{not  StepA} 

if  Abs(LJnes)  «  1  then 
if  AutoClip  then 

DrawUneClipped(X2,  Y2RefGlb  -  DeltaY,  X2,  Y2) 
else 
begin 

MoveTo(X2,  Y2RefGlb  -  DeltaY); 

UneTo(X2.  Y2); 
end; 

if  PlotLine  then 
if  AutoClip  then 

DrawLineClipped(X1 ,  Y1,  X2,  Y2) 
else 
begin 

MoveTo(X1,  Y1): 

UneTo{X2,  Y2); 
end; 

XI  :=  X2; 

Y1  :=  Y2; 
end; 
end; 

if  AxisGIb  then 
begin 

SetWindow(X1  RefLoc,  YtRefLoc,  X2RefLoc,  Y2RefLoc); 
AxisGIb  :=  false; 
end; 
end 
else 

Error(’DrawPolygon  #  1’); 
if  Flipped  then 
begin 

Temp  :=  World[WorldNdxGlb].Y1; 

WorkJIWorldNdxGlbJ.YI  :=  WorldlWorldNdxGlbl.Y2; 
World[WorldNdxGlb}.Y2  :=  Temp; 

SelectWorld(  WorldNdxGlb) ; 

SelectWind(WindowNdxGlb,  TRUE); 
end; 

end;  {  DrawPolygon  } 

procedure  Hatch{(X_1,  Y_1,  X_2,  Y_2,  Delta  :  real)}; 
var 

XI,  Y1,  X2,  Y2  :  integer; 

procedure  HatchDirect(X1,  Y1,  X2,  Y2,  Delta  :  integer); 
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var 

I,  Yst,  Yen,  Count ;  integer; 

XI  RefLoc,  X2Refljoc,  Y1  RefLoc,  Y2RefLoc  :  integer; 
ClippingLoc  :  boolean; 

X1D,  Y1D,  X2D,  Y2D  :  integer; 

begin  {  Hatch  Direct } 
if  Delta  <>  0  then 
begin 

HatchGb  >  true; 

ClippingLoc  >  ClippingGIb; 

ClippingGIb  ;=  true; 

XI  RefLoc XIRefGIb; 

XIRefGIb  >  XI; 

X2RefLoc  :=  X2RefGlb; 

X2RefGlb  :=  X2: 

Y1  RefLoc  :=  YtRefGIb; 

Y1  RefGIb  :=  Y1 ; 

Y2RefLoc  :=  Y2RefGlb: 

Y2RefGlb  :=  Y2; 

Yst  :=  Y1  +  Delta; 

Yen  :=  Y1  -  X2  +  X1  +  Delta; 

if  Delta  <  0  then 

begin 

Delta  ;=  -Delta; 

I  :=  Yst; 

Yst  :=  Yen; 

Yen  :=  I; 
end; 

Count  :=  (Y2  -  Y1  +  X2  -  X1  +  X2  -  X1)  div  Delta; 
for  I  :=  1  to  Count- 1  do 


begin 

X1D 

:=  X1; 

Y1D 

:»  Yst; 

X2D 

:=  X2; 

Y2D  : 

Yen; 

if  Clip(X1D,  Y1D,  X2D,  Y2D)  then 
begin 

MoveTo(X1D,  Y1D); 

LineTo(X2D,  Y2D); 
end; 

Yst :«  Yst  +  Delta; 

Yen  :«  Yen  +  Delta; 
end; 

ClippingGIb  ;>  ClippingLoc; 
HatchGIb  >  false; 

XIRefGIb X1  RefLoc; 
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X2RefC:lb X2RefLoc: 

YIRefGib YIRefLoc; 

Y2RefGlb Y2Reflj0c: 
end; 

end;  {  HatchOirect } 
begin  {  Hatch } 

HatchDirect(trunc(X_1),  trunc(Y_1),  trunc(X_2),  trunc(Y_2),  trunc(Delta)) 
end;  {  Hatch } 

procedure  DrawHistogram{{A  rPlotArray;  NPoints  :  integer; 

Hatching  :  boolean;  HatchStyle  :  integer)}; 


var 

XI,  X2,  Y2.  NPixels,  Delta,  NDiff,  YRef,  LineStyleLoc,  I  :  integer; 
Fract,  S,  Y  :  real; 

DirectModeLoc,  Negative ;  boolean; 

X1  Loc,  Y1  Loc,  X2LOC,  Y2Loc  :  integer; 

X1  RefLx>c,  Y1  RefLoc,  X2RefLoc,  Y2RefLoc,  YAxis  :  integer; 

Temp  :  real; 

begin  {  DrawHistogram  } 
if  ABS(NPoints)  >=  2  then 
begin 

LineStyleLoc  ;=  LinestyleGIb; 

SetLineStyle(O);  {  Black } 
if  AxisGIb  then 
begin 

XI  RefLoc  :=  Window[WindowNdxGlb].X1 ; 

YIRefLoc  :=  Window[WindowNdxGlbi.Yl ; 

X2RefLoc  :=  Window(WindowNdxGlbj.X2; 

Y2RefLoc  ;=  Window[WindowNdxGlbj.Y2; 

SetWindow(X1  RefGIb  +  X1  Offset,  YIRefGib  +  Y1  Offset, 
X2RefGlb  -  X20ffset,  Y2RefGlb  -  Y20ffset); 

AxisGIb  :=  TRUE; 
end; 

Negative  >  NPoints  <  0; 

NPoints  ABS(NPoints); 

NPixels  X2RefGlb  -  X1  RefGIb; 

Delta  >  NPixels  div  NPoints; 

NDiff  >  NPixels  -  Delta  *  NPoints; 

Fract  :■  NDiff  /  NPoints; 

S  :=  -Fract; 

X1  :=  XI  RefGIb; 

Temp  >  Y2RefGlb  +  YIRefGib  -  AyGIb; 
if  Temp  >  Maxint  then 
Temp :«  Maxint 
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else 

if  Temp  <  -32767  then 
Temp  :=  -32767; 

YRef  :=  trunc(Temp): 
if  Negative  then 

DrawStraight(X1 ,  X2RefGib,  YRef); 

YAxis  >  Y1  RefGib; 
if  BYGIb  >  0  then 
YAxis  Y2RefGlb; 
for  I  :=  1  to  NPoints  do 
begin 

X2  XI  +  Delta; 

Y  :=  A{l.  2]; 

if  not  Negative  then 

Y  :=  ABS(Y); 

Temp  AyGIb  +  ByGib  *  Y; 
if  Temp  >  Maxint  then 
Temp  ;=  Maxint 
else 

if  Temp  <  -32767  then 
Temp  :=  -32767; 

Y2  :=  Y2RefGlb  +  Y1  RefGib  -  trunc(Temp); 
if  not  Negative  then 
begin 

MoveTo(X1,  YAxis); 
yneTo{X1,  Y2); 

MoveTo{X1 ,  Y2); 

UneTo(X2,  Y2); 

MoveTo(X2,  Y2); 

LineTo(X2,  YAxis); 
if  Hatching  then 
if  Odd(l)  then 

Hatch(X1 ,  Y2,  X2,  YAxis.  HatchStyle) 
else 

Hatch(X1,  Y2,  X2,  YAxis.  -HatchStyle); 
end 
else 
begin 

MoveTo(X1.  YRef); 

UneTo(X1.  Y2); 

MoveTo(X1.  Y2); 

UneTo(X2.  Y2); 

MoveTo(X2,  Y2); 

UneTo{X2.  YRef); 
if  Hatching  then 
if  YRef  -  Y2  <  0  then 
if  Odd(l)  then 


156 


Hatch(X1,  YRef,  X2.  Y2,  HatchStyle) 
else 

Hatch(X1.  YRef,  X2,  Y2,  -HatchStyle) 
1?0dd(l)then 

Hatch(X1 ,  Y2.  X2,YRef.  HatchStyle) 
else 

Hatch(X1,  Y2,  X2,  YRef,  -HatchStyle); 

end; 

X1  :=  X2: 
end: 

if  AxisGIb  then 
begin 

SetWindow(X1RefLoc,  Y1  RefLoc,  X2RefLoc,  Y2RefLoc): 
AxisGIb  :=  FALSE; 
end; 

SetLineStyle(IJneStyleLoc); 

end 

else 

Error('DrawHistogram'); 
end;  {  DrawHistogram  } 

begin 

end.  {TurboGraph} 
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*  This  is  the  resource  file  that  defines  the  menus  and  icons  for  MacRootLocus. 
MacRootLocus.Rsrc 
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EditTextItem  Enabled 
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EditTextItem  Enabled 
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EditTextItem  Enabled 
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TYPEDUOG 


,305  (36) 

Characteristic  Equation  Coefficient  Data 
75  56  285  456 
Visible  NoGoAway 
1 
0 

305 

TYPE  DITL 
,305 
15 

Btnitem  Enabled 
15  320  35  380 
CK 

Btnitem  Enabled 
50  320  70  380 
Cancel 

EditTextItem  Enabled 
110  20130  130 

EditTextItem  Enabled 
110  145130  255 

EditTextItem  Enabled 
110  270  130  380 

EditTextItem  Enabled 
165  20  185  130 

EditTextItem  Enabled 
165145185  255 

EditTextItem  Enabled 
165  270185  380 

StateTextItem  Enabled 
2015  40  305 

Characteristic  Equation  Coefficient  Data 

StateTextItem  Enabled 
85  23  105  83 


169 


S**5 

StateTextItem  Enabled 
85148105  208 
S**4 

StateTextItem  Enabled 
85  273  105  333 
S**3 

StateTextItem  Enabled 
140  23160  83 
S**2 

StateTextItem  Enabled 
140148160  208 
S**1 

StateTextItem  Enabled 
140  273  160  333 
S**0 

TYPEDUOG 
,306  (36) 

Characteristic  Equation  Coefficient  Data 
45  56  310  456 
Visible  NoGoAway 
1 
0 

306 

TYPE  DITL 
.306 
17 

Btnitem  Enabled 
15  320  35  380 
CK 

Btnitem  Enabled 
50  320  70  380 
Cancel 

EditTextItem  Enabled 
110  270130  380 


170 


EditTextItem  Enabled 
165  20185130 

EditTextItem  Enabled 
165  145  185  255 

EditTextItem  Enabled 
165  270  185  380 

EditTextItem  Enabled 
220  20  240130 

EditTextItem  Enabled 
220  145  240  255 

EditTextItem  Enabled 
220  270  240  380 

StateTextItem  Enabled 
2015  40  305 

Characteristic  Equation  Coefficient  Data 

StateTextItem  Enabled 
85  273105  333 
S**6 


StateTextItem  Enabled 
140  23  160  83 
S**5 


StateTextItem  Enabled 
140148160  208 
S**4 


StateTextItem  Enabled 
140  273160  333 
S**3 

StateTextItem  Enabled 
195  23  215  83 
S**2 


StateTextItem  Enabled 
195  148  215  208 
S**1 


StateTextItem  Enabled 


171 


195  273  215  333 
S'*0 

TYPEDUDG 
,307  (36) 

Characteristic  Equation  Coefficient  Data 
45  56  310  456 
Visible  NoGoAway 
1 
0 

307 

TYPE  DITL 
.307 
19 

Btnitem  Enabled 
15  320  35  380 
CK 

Btnitem  Enabled 
50  320  70  380 
Cancel 

EditTextItem  Enabled 
110  145  130  255 

EditTextItem  Enabled 
110  270  130  380 

EditTextItem  Enabled 
165  20  185  130 

EditTextItem  Enabled 
165  145  185  255 

EditTextItem  Enabled 
165  270  185  380 

EciiTextItem  Enabled 
220  20  240  130 

EditTextItem  Enabled 
220  145  240  255 

EditTextItem  Enabled 


172 


220  270  240  380 


StateTextItem  Enabled 
2015  40  305 

Characteristic  Equation  Coefficient  Data 

StateTextItem  Enabled 
85148105  208 
S**7 


StateTextItem  Enabled 
85  273  105  333 
S**6 


StateTextItem  Enabled 
140  23160  83 
S**5 


StateTextItem  Enabled 
140148  160  208 
S**4 


StateTextItem  Enabled 
140  273160  333 
S**3 


StateTextItem  Enabled 
195  23  215  83 
S**2 


StateTextItem  Enabled 
195  148  215  208 
S**1 


StateTextItem  Enabled 
195  273  215  333 
S**0 

TYPEDLDG 
,308  (36) 

Characteristic  Equation  Coefficient  Data 
45  56  310  456 
Visible  NoGoAway 
1 
0 

308 


173 


TYPE  DITL 
.308 
21 

Btnitem  Enabled 
15  320  35  380 
CK 

Btnitem  Enabled 
50  320  70  380 
Cancel 

EditTextItem  Enabled 
110  20130130 

EditTextItem  Enabled 
110  145130  255 

EditTextItem  Enabled 
110  270  130  380 

EditTextItem  Enabled 
165  20  185  130 

EditTextItem  Enabled 
165  145  185  255 

EditTextItem  Enabled 
165  270  185  380 

EditTextItem  Enabled 
220  20  240  130 

EditTextItem  Enabled 
220  145  240  255 

EditTextItem  Enabled 
220  270  240  380 

StateTextItem  Enabled 
20  15  40  305 

Characteristic  Equation  Coefficient  Data 

StateTextItem  Enabled 
85  23  105  83 
S**8 


174 


StateTextItem  Ensiled 
85148105  208 


StateTextItem  Enabled 
85  273  105  333 
S**6 


StateTextItem  Enabled 
140  23160  83 
S**5 


StateTextItem  Enabled 
140148  160  208 
S**4 


StateTextItem  Enabled 
140  273  160  333 
S**3 


StateTextItem  Enabled 
195  23  215  83 
S**2 

StateTextItem  Enabled 
195148  215  208 
S**1 


StateTextItem  Enabled 
195  273  215  333 
S**0 


TYPEDUX3 
,309  (36) 

Characteristic  Equation  Coefficient  Data 
40  56  315  456 
Visible  NoGoAway 
1 
0 

309 

TYPE  DITL 
,309 


175 


Btnitem  Enabled 
15  320  35  380 
CK 

Btnitem  Enabled 
50  320  70  380 
Cancel 

EditTextItem  Enabled 
75  20  95  130 

EditTextItem  Enabled 
130  20150130 

EditTextItem  Enabled 
130  145150  255 

EditTextItem  Enabled 
130  270  150  380 

EditTextItem  Enabled 
185  20  205  130 

EditTextItem  Enabled 
185145  205  255 

EditTextItem  Enabled 
185  270  205  380 

EditTextItem  Enabled 
240  20  260  130 

EditTextItem  Enabled 
240  145  260  255 

EditTextItem  Enabled 
240  270  260  380 

StateTextItem  Enabled 
2015  40  305 

Characteristic  Equation  Coefficient  Data 

StateTextItem  Enabled 

50  23  70  83 

S**9 


StateTextItem  Enabled 


176 


105  23125  83 
S**8 


StateTextItem  Enabled 
105148125  208 
S**7 

StateTextItem  Enabled 
105  273125  333 
S**6 

StateTextItem  Enabled 
160  23180  83 
S**5 

StateTextItem  Enabled 
160148180  208 
S**4 

StateTextItem  Enabled 
160  273  180  333 
S**3 


StateTextItem  Enabled 
215  23  235  83 
S**2 


StateTextItem  Enabled 
215  148  235  208 
S**1 


StateTextItem  Enabled 
215  273  235  333 
S**0 


TYPEDUDG 

,310(36) 

Characteristic  Equation  Coefficient  Data 
40  56  315  456 
Visible  NoQoAway 
1 
0 

310 

TYPE  DITL 
,310 


177 


25 


Btnitem  Enabled 
15  320  35  380 

a< 

Btnitem  Enabled 
50  320  70  380 
Cancel 

EditTextItem  Enabled 
75  20  95  130 

EditTextItem  Enabled 
75145  95  255 

EditTextItem  Enabled 
130  20  150130 

EditTextItem  Enabled 
130145  150  255 

EditTextItem  Enabled 
130  270  150  380 

EditTextItem  Enabled 
185  20  205  130 

EditTextItem  Enabled 
185  145  205  255 

EditTextItem  Enabled 
185  270  205  380 

EditTextItem  Enabled 
240  20  260130 

EditTextItem  Enabled 
240  145  260  255 

EditTextItem  Enabled 
240  270  260  380 

StateTextItem  Enabled 
2015  40  305 

Characteristic  Equation  Coefficient  Data 


178 


StateTextItem  Enabled 

50  23  70  83 

S**10 


StateTextItem  Enabled 
50  148  70  208 
S**9 

StateTextItem  Enabled 
105  23125  83 
S**8 

StateTextItem  Enabled 
105148125  208 
S**7 

StateTextItem  Enabled 
105  273  125  333 
S**6 


StateTextItem  Enabled 
160  23180  83 
S**5 


StateTextItem  Enabled 
160148180  208 
S**4 


StateTextItem  Enabled 
160  273  180  333 
S“3 


StateTextItem  Enabled 
215  23  235  83 
S**2 


StateTextItem  Enabled 
215  148  235  208 
S**1 


StateTextItem  Enabled 
215  273  235  333 
S**0 


TYPE  MENU 

,1000 


179 


\14 

About  MacRootLx)Cus ... 
(- 

,1001 

File 

EQ  Parameter/E 
Get  Coefl/G 
Print  Screen  /S 
Print  Window  /R 
Quit/0 

.1002 

Edit 

Undo  /U 
(- 

Cut  /X 
Copy  /C 
Paste  N 
Clear 

,1003 

Plot 

One  Parameter/0 
Two  Parameter/! 

.1004 

Help 

EQ  Parameter/A 
Get  Coeff/F 
One  Parameter/N 
Two  Parameter/W 
Print  Out/P 


180 


LIST  OF  REFERENCES 


1.  Apple  Computer,  Inc., Macintosh  System  Software  User’s  Guide,  pp. 
123-154,1988. 

2.  Graham,  N.,  Introduction  to  Computer  Science:  A  Structured  Approach,  2nd 
ed.,  pp.  239  —  249,  West  Publishing  Company,  1982 

3.  Ralston,  A.,  A  First  Course  in  Numerical  Analysis,  pp.  368  -  371,  McGraw- 
Hill  Book  Co.,  1978. 

4.  Borland  International,  Inc.,  Turbo  Pascal  Tutor  (a  Self-study  Guide  to 
Turbo  Pascal  on  the  Macintosh),  1987. 

5.  Borland  International,  Inc.,  Turbo  Pascal  for  the  Mac  (User’s  Guide  and 
Reference  Manual),  1986. 


181 


INITIAL  DISTRIBUTION  LIST 


NO. Copies 

1.  Defence  Technical  Information  Center  2 

Cameron  Station 

Arlington,  Virgini  22304-6145 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  CA  93943—5002 

3.  Chairman,  Code  62  1 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey  CA  93943—5000 

4.  Professor  George  J.  Thaler,  Code  62Tr  5 

Naval  Postgraduate  School 

Monterey  CA  93943—5000 

5.  Professor  Harold  A.  Titus,  Code  67Ts  1 

Naval  Postgraduate  School 

Monterey  CA  93943—5000 

6.  Professor  SeHung  Kwak,  Code  52Kw  1 

Naval  Postgraduate  School 

Monterey  CA  93943—5000 

7.  Professor  James  W.  Bamham  1 

Assoc.  Head.  Dept  of  Mech  Engineering  k  Mechanics 

College  of  Engineering  Drexel  Univ. 

Philadelphia  PA,  19104 

8.  Major  Ko,  Sung  Hoon  5 

286-38  SooYu  4  Dong  DoBong  Gu  Seoul  .^32 

Seoul  Korea 

9.  Park,  Seung  Chin  1 

3401  N.  Columbus  #16j 

Tucson  AZ.  85712 

10.  Kim  Yoo  Chang  1 

460  W.  Forest  ave  #903 

Detroit  Michigan  48201 


182 


1 


1 1 .  Lcdr  Kenneth  Mac  Donald 
1524  Dolphin  Court 
Orange  Park  Florida  32673 


12.  Lcdr  Roy  Wood  1 

1522  Oak  Knoll  Road 

Virginia  Beach,  Virginia  23414 

13.  Lcdr  Hwang,  Jung  Sub  1 

SMC1209  Naval  Postgraduate  School 

Monterey,  CA  93943 

14.  Seo  Yong  Seok  1 

SMC1448  Naval  Postgraduate  School 

Monterey,  CA  93943 

15.  Hong-on  Kim  1 

SMC2665  Naval  Postgraduate  School 

Monterey,  CA  93943 

1C.  Kang,  Mung  Hung  1 

SMC1375  Naval  Postgraduate  School 
Monterey,  CA  93943 

17.  Kwon,  Hui  Man  1 

SMC  1375  Naval  Postgraduate  School 

Monterey,  CA  93943 

18.  Yang,  Young  Hyee  1 

SMC2440  Naval  Postgraduate  School 

Monterey,  CA  93943 


183 


